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Die Biegeschwingungen einseitig eingespannter gekriimmter 


Stabe und Rohre. 
Von W. Wuest. 


1. Einfiihrung. Kreisformig gekrimmte Rohre von ovalem Querschnitt werden in der Druck- 
meBtechnik in Gestalt der. sogenannten Bourdonrohre! angewandt. Das Schwingungsverhalten. 
solcher Rohrenfedern interessiert deswegen, weil nicht selten beim Einbau der Druckmesser 
an Stellen, die starken Erschitterungen ausgesetzt sind (z. B. bei Pumpen, Kompressoren, 
Dieselmaschinen usw.) Resonanzschwingungen auftreten, die sehr unerwiinscht sind. Ferner ist 
es wertvoll, die Eigenschwingungszahl von Réhrenfedern in denjenigen Fallen zu kennen, in 
denen diese bei Registriergeraten zur Aufnahme von Druckschwingungen als MeBorgane dienen. 
In der Technik treten Schwingungsaufgaben auf, die eine gewisse Ahnlichkeit mit dem hier 
behandelten Problem haben?. Inshesondere ist auf die Schwingungen der Rahmen von Wechsel- 
strommaschinen sowie auf die Schwingungen von Kreisbogentragern im Briickenbau hinzu- 
weisen. Bei diesen Aufgaben ist der Kreisbogen beiderseitig fest eingespannt oder gelenkig 
befestigt, wahrend in unserem Fall die Réhrenfeder einseitig fest eingespannt ist und das andere 
Ende frei schwingen kann, soweit keine Behinderung durch Ubertragungsorgane stattfindet. 

Die Ermittlung der Eigenschwingungszahlen kreisbogenformig gekrimmter Stabe wurde 
erstmalig von H. Lamb? in Angriff genommen und unter der Annahme geringer Anfangskrimmung 
und dehnungsloser Stabachse genahert durchgefiihrt. Wegen der rechnerischen Schwierigkeiten 
bei beliebig groBem Offnungswinkel wendet I. P. Den Hartog+ das Naherungsverfahren von 
Rayleigh- -Ritz auf den Kreisbogen bei festgehaltenen Enden an. K. Federhofer® beriicksichtigt 
noch eine radiale Belastung des Kreisbogens. Dabei taucht aber die Frage auf, ob nicht gerade 
die vernachlassigten Achsdehnungen von entscheidendem Einflu8 auf das Schwingungsverhalten 
sind. Diese Frage wird von F. W. Waltking® untersucht, der fiir den Zweigelenkbogen die Lésungen 
der vollstandigen Differentialgleichung unter Beriicksichtigung der Achsdehnung entwickelt. 
Bei den gegensymmetrischen Schwingungen (Schwingungsknoten in der Mitte) ist die Achs- 
dehnung ohne EinfluB, bei den symmetrischen Schwingungen macht sich dagegen ein wesentlicher 
Einflu8 geltend. Wenn der kreisbogenformig gekrimmte Stab nur am einen Ende fest ein- 
gespannt ist, darf man sicher annehmen, daB die Achsdehnungen keinen merklichen EinfluB 
auf das Schwingungsverhalten haben. Aus diesem Grunde beschrankt sich die poly poe Unter- 
suchung auf den dehnungslosen Fall. 


re Ableitung der Frequenzgleichung. Der normale Be- 
wegungsablauf der Réhrenfeder als Druckmefglied spielt 
sich ausschlieBlich in der Ringebene ab. Aus diesem Grunde 
beschranken wir uns im folgenden auf Biegeschwingungen 
in der Ringebene, zumal diese wegen der geringeren Steifig- 
keit der Réhrenfeder in dieser Richtung leichter angeregt 
werden. Allgemeinere Ansatze sind von K. F ederhofer’ 
angegeben worden. 

‘Im folgenden sei mit u die N. ormalverschiebung, mit v. 
die Tangentialverschiebung eines Ringelementes bezeich- 
net. Ferner sei die Masse in der Langeneinheit ° des 7+d7 : 
Umfanges, und r der Krimmungshalbmesser der neutralen M 
Faser. sBezeiclinet man- suBerdem eee none Abb. Oba CE Krifte und Momente an einem Federelement, 
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In der Tat kann man bei derartigen Roh- 
renfedern Lissajoussche Figuren dieses | 
Frequenzverhaltnisses beobachten. We- 
gen des nicht genau ganzzahligen Verhalt- 
nisses werden die verschiedenen Figuren 
rasch durchlaufen. Bei gréBeren Schwin- 
gungsamplituden bemerkt man, da die 
Hauptschwingung auf einer nach aufen — 
(also entgegengesetzt zur Krimmung der 
Réhrenfeder) gebogenen Bahn verlauft. 
Die schnellere Nebenschwingung wird . | 
rasch gedampft und tritt daher weniger 
in Erscheinung. 


; 4. Die gleichmaBig mit Masse belegte 
Windungswinke! > Feder, Wir wenden uns jetzt der allge- 
@ esha a HE mais meinen Untersuchung von Gleichung (9) 
zu, die eine Gleichung dritten Grades 
fiir n2 darstellt. Indem wir n?=x setzen, 
tragen wir zundchst die Kurve y=x (1 
+ x)? auf (Abb. 3). Dann legen wir durch jf 

den Punkt y=0, x=1 ein Geradenbiischel y=k(x—1) mit k als Parameter. Die Schnittpunkte’ jj 
faa der Geraden mit der Kurve dritten Grades ergeben die Lésungen x,, x9, x, Wie man aus Abb. 3 jf 
CS ae erkennt, stellen sich Grenzfalle dadurch ein, daB die Gerade die Kurve beriihrt. Durch diese 
Grenzlagen wird das Geradenbischel in drei 
Bereiche getrennt, wobei im mittleren Be- 
reich (II) nur ein reeller Schnittpunkt vor- 
handen ist. 
Analytisch sind die Grenzfalle gegeben durch 


Freguenzverhaltnis 


Abb. 2. Verhaltnis der Frequenzen der Neben- und Haupt- 
schwingung einer masselosen nur am Ende belasteten Feder. 


7 
(8 Tk+2)=0 (16) 
mit den Lésungen 
Bae 
k eta RG — 0,11340, 
os py 175636 6 


Kine der Wurzeln von n?= x ist stets reell 
negativ, wahrend die anderen beiden ein ver- 
schiedenartiges Verhalten zeigen, wie aus 
Tabelle 2 ersichtlich ist. In dieser Tabelle — 
sind zur Kennzeichnung der verschiedenen — 
Lésungsbereiche und Grenzlagen die gleichen 
Symbole wie in Abb. 3 verwendet (1. Spalte). 
Im Lésungsbereich II sind die Wurzeln n? 
und n; konjugiert komplex. Uber die Be- 
deutung der ebenfalls in Tabelle 2 angegebe- 
nen GréBen A, u, ¥, UL, V, &, f, siehe nachste 

~ Seite. ; 

Am fest eingespannten Federende p=g* lauten offenbar die Randbedingungen wie im vor- 
hergehenden Abschnitt . 


__Abb. 3. Graphische Lésung der Gleichung (9) 
mit Einzeichnung der Grenzlagen. 


DSU = 0, SSO oad eae . AOA RY 
Die Striche bezeichnen hierbei wiederum die Ableitungen nach y. Am freien Ende y~=0 muB 


andererseits das Moment, die Normal- und Tangentialkraft verschwinden. Aus (5). (6 
und (12) folgt dann fir y=0 the tip a : Se 


y ; vo”! Bye ie 


a hl ie | ; 
“. ‘ 8 
vty =ky’ . 1) 


= Z - : r : 
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Tabelle 2 


ocak der Gleichung (9). 
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| | 


k ni 4 A i v 
Bee: Loy aE ae | 1 Vee | eevee Vk 
0,025 —1,2135 —0,7593 | —0,027131 1,1016 0,8714 0,16471 
0,050 —1,29756 —0,64.246 —0,05999 |  1,13910 0,80153 0,244.92 
I 0,075 —1,36072 —0,53656 | —0,10274 1,16650 0,73250 0,32055 
_| 0,100 —1,41334 —0,41714 —0,16964 | 1,18884 0,64587 0,41287 
0,1125 —1,43691 —0,31319 | —0,25000 1,19867 0,55964 0,50000 
b 0,11340 —1,43849 | —0,28076 —0,28076 1,1994 0,5298 0,5298 
k | tas A 4 B 
0,15 -—1,50135 —0,24932 +0,1938 i 1,2253 0,18256 0,5320 
0,30 —1,6909 —0,15455 § +0,3918 i 1,3003 0,3655 0,5365 
Il 2/3 =o + VIP i 1,41421 0,53729 0,53729 
3 33 0,5 + [3/4 i 1,73205 0,86603 0,50000 
36/5 =A, 1,0 + V4 i 2,00000 1,08195 0,41369 
40/3 5 bs + 5/12 i 2,23607 1,25160 0,25786 
c | 17,6366 —5,5617 17808 67. 2,3583 1,3345<~ 0 
k nt n3 n3 A Le y 
. 150/7 — 6,0000 | 1,34535 2,65465 2.4495 1,1600 1,3393 
100 —11,4291 1,04360. 8,3853 83,3807 1,0218 2,8957 
IIT |~ 211, 765. | —16,0000 1,01964 12,9804 4,0000 1,00981 3,60283 
1000 —33,1070 1,00404 30,0927 5,7539 1,00202 5,4857 
3000 56,2577 1,00133 53,2558 7,5004 1,000668 71,2976 
Si fate ey ae hee 
d k—-00 GVE | 1 | Vk Vk 1 Ve 


¥ 


Die Lésungen haben in den drei in Abb. 3 gekennzeichneten Lésungsbereichen eine perros scene 
Form. Wir behandeln daher die drei Bereiche. getrennt: 


a) Lésungsbereich I: 0Sk0,11340. Alle sechs Wurzeln der Gleichung (9) sind rein 
prcear: eo zwar paarweise mit entgegengésetztem Vorzeichen gleich. Wir setzen daher 


=41A, 


Ng = +ip, 


Ny = 41y.- 


E Die Werte A, yw, » finden sich in den drei rechten Spalten der Tabelle 2. Die Lésung hat dann 
im Bereich I folgende Form: 


3 
= 24 (Am COS Om rt Bot sin Om P) cos(wt+e), 


‘eas 


=A : 


Og=/4 , 


6 


6.=V- 


4 


(19) 


. “ist. Die Rednibedinguigen (17). und (18) liefern ein System von sechs Gleichungen. fir die Am 
und Bm. Das Verschwinden der Determinante der Koeffizienten ergibt die ,,Frequenzgleichung“, 


_wobei i im folgenden wieder statt y* ertecher, ue eras wird: 


ae eagle Ap sin Ay "eos Lp sin UY cos vp sin yy 
eee —Asin 2G a cosAp —psin LP pcos up —Vvsin vp y Cos Vp | - 
| -Aeos Ap, —A?sindp —p? cos wp = —p? sin up —v cos yp —v® sin yp 
eae (L—A4) 2 2 y hone pw) Buea) y? (1 —y?) 0 =0 (20), 
(ee es y(L—p?) 0 » (1-92) : 
eee ee ALA), 0 3 (1 —?) 0 v8 (1 —y?) ‘\ 
leaner tees A pat eae Seren oa Ty 1+ y2 
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Die Auflésung dieser Determinante fihrt auf folgende Gleichung: 
a, cos Ay-+-a, cos up + ay cos vp + b cos Ay cos UP Cos VP + C4 COS Ags sin wy sin ag = | (21) 
+ cy, cos wg sin vp sin Ay-+ cy cos yp sin Ap sin up=0. 
Hierbei ist : 


is: y 1—?? wo lw 

a a, =2 A, Ay, ca = —Ay Ay (7) eae au lo)? 

ae A yA 

ae a, =2 Ay Ax, |: b= AP+AR+ A}, | Gu=—4y a (| Ga + ope) 5) 

= Pigs Wary yee GN ve < 3 
ay = 2 A, Ay, cy = —Ay Ay (5 1— 72 ay ary, ar 


_Dabei haben wir zur Abkiirzung gesetzt 
Ay = (0° m2) 42 (12), | 
Au = (22-22) w® (Lp), (23) 
Ay = (u?—2?) v7 (1—v?). 
b) Losungsbereich II: 0,11340 Sk = 17,6366. Im zweiten Lésungsbereich sind nz und 
n2 konjugiert komplex. Daher werden auch ~ und » konjugiert komplex, und wir setzen 
w2=A+iB, p=at+ifp, 
y= A=-iB, v=a-ip. 
_ Die A liencnacnas (21) geht damit iber in . 
ae a, cos ap Gof By + a, sin ag Gin By + 
+ by cos Am + b, cos Ay cos? ay + by cos Ap of? Bye + 
+c, sin Ag cos ap sin ay + ¢, sin Ay Gof Py Gin Bo = 0, 
wobei 


ead ee 4A, Ay, 
a A2 3 2 2 2 2 2 ames 
ibys Ate doe, By = | AP +2(Aj—A; 424 + Ai) k= 


: Pi ae pe it 
a= 3 tn (A, -14,) (2 + 48-4 A-iB) w(1—p) 


®t und s kennzeichnen hier Real- und Imaginarteil des Klammerausdruckes. Ferner ist © | 
CABS ot pero I Bee Ge ae oA 

aa is = eee [4? + B?—A-+ 42 (1 say), «4 
ts = (1—A) (A? + B2~ A?) — B22, 
. N = (a? + 6?) (1 —A)P + B*). sg 
; c) Lasubetberdat Ill: 17,6366 Sk ce, co. Im Ebranerbereicl Il mid ny und Ls te | 
é _ Positiv, oe wir setzen daher niet ho Bl 
i 3 psi, vie, Pi, 
i Die F paciienadlolahiily (21) geht in leicht ersichtlicher Weise uber in i 


Se : 


ohm 


@ cos hp + a, Cof “P+ ay Got ?p +5 cos Ay Col Zp Cof » Dp ee Ng 
Dae €1 COs ua Gin Bp aie Vp + Ze sin Ap nal ue Gin Vp +e cy sin Age Sin wep Gof 7p J 


aid eee teeny ano Ay ss 
ey De en ie wy 
wy re aA 72 A2 A: 4 =. Fe iv A i; ad v 147? * 1 
vn, | b= Att Ai + 4S, GUM Ocal ah 
wih { 4 mw 1+? A 1-2 
A4ys : i, G=—-AA | 1-P ees aE ey 
gh hy Cc tae ay 0-2), Be ia et en ee 
Be Tale Sa = — (22-99) fi? (Vay Wai Lg ety OCTET 
fate uae Lk ate ae eae ye 7) 
si \ : 4 . ; \ : a 
i i Fis 3 ‘, 
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Tabelle 3. Wurzeln der Frequenzgleichung. 


te Po | Bo Pi. | By Po Be 
0,025 551° 3,83 
0,05 301,7° 2,490 | 474° 3,91 564° 4,66 
0,075 247,9° 2,265 390° 3,57 520° 4,75 
0,10 223,3° 2,184 
0,1125 213,0° 2,156 348° 3,52 477° 4,83 
0,15 193,14° 2,098 323° 3,51 454° 4,94 
0,30 156,28° 2,020 268,5° 3,47 401,22 5,19 
0,667 124,17° 1,956 
3 SoU be L018 163,1° 3,75 266,2° 6,12 
7.2 | _ 66,45° 1,901 136,5° 3,91 226,0° 6,46 
21,43 50,50° 1,890 104,0° 4,09 181,6° 6,83 
211,77 28,23° 1,881 66,00° 4,390 110,5° 1,35 
1000 19,15° 1,879 46,25° 4,545 17,14° 7,58 
3000 14,54° 1,876 35,66° 6,610 59,23° 7,66 
oo 0 1,875 0 4,694 0 7,855” 
k | Ps Bs Ps | Ba | Ps Bs 
0,05 1073° 8,86 
0,075 816° 7,46 a 
0,10 674 SS 6,62 870° 8,55 987° 9,70 
0.1125 |: 652° 6,60 826° 8,37 958° 9,69 
0,15 612° 6,25 759°. 8,24 906° he One 
0,30 539° 6,97 677° _ 8,76 816° 10,55 
3 370,0° 8,50 474° 10,89 578° 13,27 
7,2 315,7° 9,03 405,7° 11,60 496° 14,18 
21,43 255,0° 9,59 328,5° 12,35 402° 15,11 
211,77 155,5° 10,34 200,50) >| 13,35 245,5° 16,34 
1000 108,4° 10,65 139,7° 13,72 171,0° 16,80 
3000 83,23° 10,17 107,2 13,86 131,2° 16,95 
00 0 10,996 0 14,137 0 11,279 


In den Grenzfallen k—>0 und k—>oo nimmt die Frequenzgleichung eine besonders einfache 
Form an, namlich fir 


k—>0 : cos Vko [1 + cos 2k @] =O, (26) 
; 4 4 
k— oo: 1+ cos Vk p Goji ky a= OG (27) 


Die letzte Formel stimmt aber iiberein, mit der F requenzgleichung fir den einseitig eingespannten 
geraden Stab. 

Fir die praktische Durchfihrung der Rechnung hat es sich als zweckmaBig erwiesen, von 
gegebenen Werten von k auszugehen. Damit sind auch die Wurzeln A, uu, » usw. gegeben. Je nach 
dem Lésungsbereich werden dann aus der Frequenzgleichung die zugehérigen Windungswinkel 
durch schrittweise Annaherung errechnet. Tabelle 3 enthalt die Rechenergebnisse fiir die Grund- 
frequenz und die ersten finf Oberschwingungen. Es empfiehlt sich, die Eigenfrequenz durch 
folgende Formel anzugeben: . 


By ES f 


Hierbei ist l=rq die abgewickelte Feéderlange, gFl die Gesamtmasse der Feder. Zwischen fj 
und k besteht die Beziehung 


bee Cw 


In Tabelle 3 sind daher auch die GréBen fj angegeben. In Abb. 4 ist fj tber dem Windungs- 
winkel aufgetragen. Die Werte titber p=3060° sind praktisch nicht zu verwirklichen, wenn man 
sie nicht naherungsweise gewendelten oder spiralfoérmigen Federn zuordnen will. Die Grenzen 
der Lésungsbereiche erscheinen wegen (29) mit k=konst. als unter 45° geneigte gerade Linien. 


10° 20° uy? 80° 60° CD” 64° —- 1000". 


Abb. 4. Die Eigenfrequenzen einer kreisférmig gekriimmten Feder in logarithmischer Auftragung. 
) Die gestrichelten Linien grenzen die drei Liésungsbereiche ab. - 


(Eingegangen am 7. August 1948.) 
Anschrift des Verfassers: Dr. Walter Wuest, (20b) Géttingen, ZeppelinstraBe 7. 
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Zur Berechnung des Profilwiderstands. 


Von H. B. Helmbold. 


1. Grundlagen. Den Profilwiderstand umstrémter Kérper kann man auf Grund der Er- 
gebnisse der Grenzschichtforschung einigermafen zutreffend vorausberechnen, falls es sich um 
schlanke zylindrische Koérper in der stationdren, ebenen (oder auch zylindrischen), inkompres- 
siblen Strémung bei hohen Reynoldszahlen handelt, wobei die Strémung iberall an der Korper- 
oberflache entlang strémen (.,anliegen“‘) soll, wenn nur die Geschwindigkeitsverteilung langs der 
Kontur und die Umschlagstellen der Grenzschicht bekannt sind!. Im allgemeinen geniigt es, 
der Berechnung die theoretische Geschwindigkeitsverteilung der Potentialstromung um den 
gegebenen Korper zugrunde zu legen; der wirkliche Geschwindigkeitsverlauf weicht zwar vom 
theoretischen am hinteren Kérperende dadurch merklich ab, da® sich zwischen Kérper und 
Potentialstrémung die durch Reibung gebremste Grenzschicht. zwischenschiebt; wir werden 
aber die Berechnung so einrichten, daB wir die Abweichung von der Theorie nicht weiter zu 
beriicksichtigen brauchen. ; 

Im wesentlichen besteht die Abweichung darin, daB sich die Potentialstrémung am Kérper- 
heek nicht mehr.schlieBt, so daB ihre Geschwindigkeit dort nicht auf Null fallt, sondern einen 
endlichen Betrag Uy behalt. Dem entsprechen’nach der Bernoullischen Gleichung geringere 
statische Driicke am Ké6rperheck als in der reinen Potentialstromung,-fir welche seit 
D’Alembert bekannt ist, daB der Widerstand eines umstrémten Kérpers Null ist; wir erhalten 
vielmehr als sekundadre Wirkung der Flissigkeitsreibung einen Druckwiderstand neben der 
Oberflachenreibung, welche von den durch Reibung unmittelbar verursachten Schubspannungen 


~*~. 


os el * » s ° - - . . . ; oe ¢ 
herrithrt. Fur die Berechnung ist es jedoch unnétig, den Profilwiderstand in Oberflachen- : cage 
reibung und Druckwiderstand zu zerlegen, wenn wir nur den Strémungszustand in groBer Ent- lect 
fernung vom Kérper betrachten. Dieser ist dadurch bestimmt, da sich die Potentialstrémung auch af ae 


hinter dem Kérper nicht mehr schlieSt, weil die von der Grenzschicht erfaBten Flissigkeitsteile 
als Nachlauf (Kielwasser, Wirbelschleppe) dem Kérper nacheilen. Das aus dem umstrémten ~ 
KGrper, der Grenzschicht und dem Nachlauf bestehende Ganze kann hinsichtlich seiner Fern- 
wirkung als ein sog. Halbkérper (halbunendlich langer Kérper) oder auch als ein System von 
Quellen angesehen werden. Die Dicke des Halbkérpers im Unendlichen wird als die Verdrangungs- 
dicke des Nachlaufs bezeichnet und ergibt sich aus der Kontinuitatsbedingung als. 


+) 1 
a if (U5 —to) dy 


_ (U. ungestérte Geschwindigkeit im Unendlichen, u..(y) Geschwindigkeit im Nachlauf. weit 
-hinter dem Koérper, y' Koordinate senkrecht zu U,.). 
Um den Widerstand nach dem Impulssatz zu berechnen, umgeben wir den Kérper mit einer 
 Kontrollflache, die tberall in sehr groBer Entfernung vom Koérper verlauft; ihre Seitenflachen 
ollen aus Stromlinien gebildet werden, und die AbschluSflachen sollen stromaufwarts und 
‘stromabwarts gleich weit vom Korper entfernt sein und senkrecht zu den Stromlinien verlaufen. 
Die AbschluBflachen unterscheiden sich nur um die endliche Flache 03 2, wenn z die Lange des 
Kérpers senkrecht zur Stromungsebene ist; ebenso gro ist die Projektion der Seitenflachen 
in Richtung von U,. Die mit dem Halbkorper gleichwertige Quelle verursacht auf diesen end- 
lichen Flachen Druckunterschiede, die mit zunehmender Entfernung r der Kontrollflache vom 
Kérper mindestens wie 1/r verschwinden. Die ubrigen Teile der AbschluBflachen ergeben schon 
aus Symmetriegriinden die Druckresultierende Null. Demnach ist die in Richtung von U,, fallende 
Komponente des Druckintegrals auf der gesamten Kontrollflache gleich Null (L. Prandtl). Es 
 pleibt als einziger Beitrag zum Widerstand der Impulsstrom des Nachlaufs durch die Kontroll- 
- flache F hindurch zu bericksichtigen, und somit ist der Profilwiderstand des Kérpers 


We = f eus(Ua'wa) d P= 004 22 9a). ayy 


Durch die zweite Schreibweise der rechten Seite haben wir eine weitere charakteristische 
 Dicke des Nachlaufs, seine ,,[mpulsdicke“ definiert und 0, geschrieben, weil diese GréBe 


1 J, Pretzsch, Jb. 1938 der Dt. Luftf.-Forschg., S. 160. 
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sich aus den ,;[mpulsverlustdicken“ der Grenzschichten von Ober- und Unterseite des um- 
stromten Zylinders entwickelt (X Summenzeichen). Wir sehen uns also fur die Widerstands- 


berechnung auf die Berechnung der Impulsverlustdicken der Grenzschicht verwiesen. Diese © 


beginnt am Kérpervorderteil als laminare (glatt strémende) Grenzschicht, schlagt an der Um- 
schlagstelle in den turbulenten Zustand (Mischbewegung) um und setzt sich hinter dem Kérper 
als Nachlauf fort. Die Entwicklung der Impulsverlustdicke mu von der Rechnung durch diese 
drei physikalisch verschieden bedingten Zustande hindurch verfolgt werden. 

Hierfiir wahlen wir folgende Bezeichnungen: g fir die Dichte der Flissigkeit, mw fir ihre 
Zahigkeit, y= u/o far ihre kinematische Zahigkeit, x fir die Wegkoordinate oder Entfernung 
eines Konturpunktes vom Spaltungspunkt (vorderer Staupunkt) der Strémung, langs der Kontur 
gemessen, y fur den Wandabstand eines Punktes, von der Kontur senkrecht nach aufen ge- 
messen, z fiir die Lange des zylindrischen Kérpers senkrecht zur Strémungsebene und L fir 
die gréBte Lange im Profil des umstrémten Zylinders; dabei bleibt y in der Grenzschicht immer 
klein gegen L. Ferner bezeichnen wir mit U(x) die Geschwindigkeit an der Grenze zwischen Grenz- 
aoe schicht und Potentialstrémung, mit u(x,y)< U(x) die Geschwindigkeit in der Grenzschicht. 
Fir y=6 wird u=U, also ist 6 die Grenzschichtdicke. Die Grenzschicht beansprucht fur die in 
ihr stroémende, durch Reibung verlangsamte Flissigkeit einen breiteren Querschnitt und drangt 


die Potentialstrémung um einen entsprechenden Betrag, die Verdrangungsdicke 


h 
6*(x)= 77 [(U-u)dy, 


von der Wand ab; dabei ist h <6 eine konstante Lange, aber noch klein gegen L. 
Das fir uns wichtigste DickenmaB aber ergibt sich, wenn der Impulsverlust der Grenzschicht 
gegentiber der reinen Potentialstromung 


h 
Jz=z f[ ou(U—u)dy =zeU?6 
0 


gesetzt wird: die Impulsverlustdicke 
h 


(x) =a fu(U—u)ay. 


0 


“Wird noch die Geschwindigkeit in Richtung der 4uBeren Wandnormale y mit v und die 
Schubspannung mit 7 bezeichnet, dann lauten die Prandtlschen Grenzschichtgleichungen des 
stationdren ebenen Problems bei hinreichend schwacher und stetiger Wandkrimmung 

du du dU OT 


OAS yee a me aay” 


du 092 V5 
see 


Die erste von ihnenenthalt links zuerst die Tragheitsglieder, dann das Druckglied 
dp/dx =—@U dU/dx, welches der Grenzschicht von der Potentialstrémung aufgepragt wird, 
und rechts das Reibungsglied; die zweite ist die Kontinuitatsbedingung. Durch Integration 
_ dieses Gleichungssystems von y=0 bis y=h erhalt man eine der Karmanschen Impulsgleichung 
m gleichbedeutende Integralbedingung! y 


LEM rey : aie 
af Be erp en UR a ree al 7 7 (2) 


, worin %)(x) die Wandschubspannung bezeichnet. Diese Gleichung ist die allen folgenden Be- 
; rechnungen von Impulsverlustdicken gemeinsame Grundlage. Sabor 
_ Benotigt wird schlieBlich noch ein Formparameter, der die Gestalt der vorkommenden 
: Grenzschichtprofile (Geschwindigkeitsverteilungen in der Grenzschicht langs Querschnitten 
_ x=const) kennzeichnen soll; der meist verwendete Formparameter ist die Verhaltniszahl 
H= 0"/. Die Gestalt des Grenzschichtprofils selbst wird durch die Funktion u/U=f(y/9) 
-beschrieben. Hk Tasty ae Dike 3 e 
ohh Nunmehr kénnen wir zu den Einzelberechnungen iibergehen. | 


LY, 


1 E, Gruschwitz, Ing.-Arch. 2 (1931) S. 321. 


Wt 
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_ 2. Die laminare Grenzschicht. s der laminaren Grenzschicht ist die Schubspannung 
i udujdy. An der Wand (y=0) verschwinden u und » wegen der Haftbedingung, und die 


erste Grenzschichtgleichung vereinfacht sich dort zu 


a2 
eu Ft u(S4) =0 


ay 
oder in dimensionsloser Schreibweise ~ 
Oo? dU. 072 u/U 
“yp dx (ay), ee 
Diese Form der Wandbedingung zeigt, daB se : 
ian ee 
= yp dx vy U dx 


ein Hordiparuideter des laminaren Crenseohichtwobls ist. Wenn wir nun annehmen darion: daB 
das laminare Grenzschichtprofil bereits durch einen einzigen Formparameter eindeutig gekenn- 
zeichnet ist, dann besteht auch zwischen H und 4 sowie allen tibrigen Formparametern ein ein- 
deutiger Zusammenhang. Nun kennt man fiir gewisse oe Geschwindigkeitsverteilungen 
langs der Kontur, namlich U(x)=a+bx und U(x)=mx", die exakten Lésungen der Prandtl-. 
schen Grenzschichtgleichungen (Howarth und Hartree), und der Vergleich dieser Lésungen zeigt 
in der Tat, daB die erwahnte Annahme bei beschleunigter Grenzschicht (A> 0) sehr genau zu- 
rift. Unter diesen giinstigen Umstanden erscheint es dann auch zulassig, im Einzelfalle auf die | 
exakte ‘Lésung der Grenzschichtgleichungen zu verzichten und die der Gleichung-(2) entspre- 
chenden Integralaussagen der vorliegenden speziellen Lésungen zu verallgemeinern. 
-Setzen wir demgema®B fiir die Wandschubspannung der laminaren Grenzschicht t,=4(0u/0y), 
in die Integralbedingung (2) ein, so erhalten wir nach enfachen Umformungen die Differentials 
gleichung- der Jaminaren Grenzschicht? 

eed. U6 \ Oo , du 
a dx ae ay Cy )o 


Te ee 
vy dx 


CHS 


. 


8 dU au/U 
faa ist auch aan y18 0 


aber die _Einpaametig Mex laminaren Cpe aaa ist also die rechte 


ye -ein solcher. “Auf Grund Aer bee 


or, wo die Guisohlagstelle a i aoe aah des Ac ischyidiskeitane anita Xm Tee. 

y re , Grenzschicht vom Staupunkt bis zur Umschlagstelle durchweg beschleunigt 
dann ist die zugrundeliegende Annahme erfallt und die exakte Form der rechten 
) kann nach einer Untersuchung von A. Walz) in guter ‘Naherung ein fir allemal 
F unktion «—fA mit den Konstanten | a=0,441 und B ~ 4,16 ersetzt werden. — 

ee 441 ist aus a - Blasiusschen Grenzschichtberechnung fir ‘die ebene 

+ Poter ge area ange (A=0) bekannt, den Konstantenwert 8 hat * 

sind dabei die letzten Dezimalen bereits Bes Man hat © 
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3. Die turbulente Grenzschicht. Die turbulente Grenzschicht, welche nach einer kurzen, : 
in der Rechnung zu vernachlassigenden. Ubergangsstrecke an die laminare Grenzschicht an- || 
schlieBt, wird ebenfalls mittels (2) berechnet; doch sind die Zusammenhange zwischen Wand- |} 
schubspannung, Druckanstieg (Verzégerung) und Grenzschichtprofil noch nicht genau genug || 
erforscht, so daB sich nur grob orientierende Angaben machen lassen.*) Wir fuihren deshalb die |} 
Wandschubspannung 1,(x) =e 0 U?/2 mit einer mittleren Reibungszahl cy, die wir anschlieBend |}j 
noch abschatzen wollen, und den Formparameter H(x) mit einem entsprechenden Mittelwert H 
in (2) ein und erhalten dann als vereinfachte Differentialgleichung der turbulenten Grenzschicht 

dd = 6 dU & 5 

de po sea eeoL ae | 

und daraus durch Integration von der Umschlagstelle bis zur Hinterkante des umstrémten | 

Korpers die Impulsverlustdicke am Kérperheck? | | 
Ly 


U, \2+H Sere fe UT eee 

ons (Ge) (+H f(a) ae). (1) 
Hier ist Uy die Potentialgeschwindigkeit an der Hinterkante, welche, wie eingangs bereits be- | 
merkt, in der wirklichen Stromung einen endlichen positiven Wert besitzt. Wenn man trotzdem | 
bei der Berechnung des Integralausdrucks die theoretische Geschwindigkeitsverteilung einer rei- | 
bungslosen Kérperumstrémung zugrunde legt, so begeht man keinen schwerwiegenden Fehler, 
da sich der Unterschied zwischen dem wirklichen und dem angenommenen Verlauf von U(x) | 
auf einen kleinen Bereich nahe der Hinterkante beschrankt. Den Wert von @, itbernimmt man jf 
aus (5). im 

Der Mittelwert H des Formparameters der turbulenten Grenzschicht am umstrémten zy- 
lindrischen Kérper hangt tber cy indirekt noch von der Reynoldszahl] ab; nach den bisherigen 
Erfahrungen? ergaben sich die Wertepaare c—0,0040; H=1,5 und cpA=0,0034; H=1,4. 

Um nun den noch fehlenden Mittelwert cg der Reibungszahl im Turbulenzteil der Grenz- 
schicht am umstrémten zylindrischen Kérper, der sich aus den bereits genannten_ 
_ Griinden nicht exakt bestimmen 1]a4Bt, in moglichst guter Annaherung abzuschatzen, ersetzen 
wir ihn durch den Mittelwert der Reibungszahl im Turbulenzteil der Grenzschicht an der 
glatten ebenen Platte, die in ihrer Ebene angestrémt wird; die Reynoldszahl und die Lage 
der- Umschlagstelle an der Platte sollen dabei dieselben Werte wie am umstrémten Zylinder 
haben. 

Die Potentialgeschwindigkeit ist langs der ganzen Platte U=konst.=U,,, und R= U, L|y 
ist die Reynoldszahl der Platte von der Gesamtlange L. Bekannt sind die Formeln fir die jf 


mittlere Reibungszah] der Laminarstrecke von. der Lange x,=é,:L und der Reynoldszahl 
R,= Un %u/v = uh f - 


a , 1,328 z 
Cfu = 2, Ve = TR. nach Blasius 
und fir die mittlere Reibungszahl bei vollturbulenter Grenzschicht (Turbulenz von der Vorder- 

bis zur Hinterkante) ; . 


a 


0,427 

% GD) eg 0,407)264 

Um-hieraus die benétigte mittlere Reibungszahl im Turbulenzteil der Grenzschicht hinter 

einer laminaren Anfangsstrecke zu berechnen, denken wir uns das Plattenvorderteil bis zur _ 

Umschlagstelle durch ein kirzeres Plattenvorderteil mit turbulenter Grenzschicht ersetzt, das _ 

denselben Widerstand hat und deshalb tiberall dahinter den gleichen Beitrag zur Impulsverlust- 

_ dicke ergibt®. Mit 9%, bezeichnen wir die Reynoldszahl des turbulenten Ersatzstiickes allein © 

und mit 3% die Reynoldszah] der Platte mit dem turbulenten Ersatzstiick. Dann ist wegen der 
verlangten Widerstandsgleichheit von Ersatzstiick und Laminarteil 


Rto Cft (Reo) = Ru Cfu » 


nach , Schultz-Grunow. 


*) (Zusatz bei der Korrektur:) Einen neuen Fortschritt hierin bringt die Unt h udwi 
und W. Tillmann, Ing.-Arch. 17 (1949) S. 288. — ie She. dren ae ney eye 


1 A. Buri, Diss. Eidg. Techn. Hochsch. Ziirich 1931. 
. E. Gruschwitz, siehe FuBnote von S, 274 sowie H. Peters, Journ. of Aer. Sciences 1935, S. 7. 
L. Prandtl in W. F. Durand, Aerodyn. Theory, Bd. 3. Im Ergebnis etwas abweichend. 


f 
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woraus tj zu bestimmen ist; damit wird 


Ne = RR + Rto = #(1 —é,) ae Mto« 


Weiter folgt aus der Widerstandsgleichheit der Platte mit laminarer Anfangsstrecke und 
der Platte mit Ersatzstiick die Beziehung 


RK Cf = Re Cft (Re) ’ 


welche zur Ermittlung der mittleren Reibungszahl c¢ der Platte mit laminarer Anfangsstrecke 
dient. Der gesuchte Mittelwert ¢f des Turbulenzteils der Grenzschicht an der Platte mit lami- : 
narer Anfangsstrecke folgt dann aus rare 


ef = Su cfu + (1—Eu) eyes 
er ist in Zahlentafel 1 tabuliert. 


5 Zahlentafel 1. ae : 
CC OO OO oS ————ESSSSSSeSeSese mith 
K 108 Cft i i 
ea) 0,1 0,2 0,3 0,4 0,5 0,6 Oi blo 0.8 0,9 — oo 


1-108 4,55 4,31 4,31 4,34 4,38 A Ad. 4,52 4,62 4,74 5,10 j 
2 3,95 3,80 3,81 3,85 3,89 3,96 4,03 4,13 4,25 4,50 ae 
5 3,34 3,26 3,27 3,31 3,37 | 3,43 | 3,50 3,57 3,66 3,80 “A 
1-107 2,96 2,91 2,92 2,96 3,01 3,07 3,15 8.455. 1 A538 3,41 
2 2,61 2,59 2,61 2,64 2,69 2,74 2,82 2,90 3,00 3,11 
5 2,24 | 2.24 2,26 2,29 2,33 2,38 2,44 253 S128 63 SONOS 
PPOs > 2,02 2,02 2,04 2,07 2,11 2,16 9,22 |» 2:30 2,39 2,50 


~ 


4. Der Nachlauf hinter dem Kérper. Die Grenzschicht setzt sich hinter dem umstrémten 
Kérper als Nachlauf fort. Auch hier ist die Gleichung (2) zugrunde zu legen. Mit dem Fortfall 
der begrenzenden Wand verschwindet die Wandschubspannung (t,=0) und (2) geht in die 
Differentialgleichung des Nachlaufs iber: 3 Me 


i ¢ 


do ben OU 
<< =-(2+H) 4. (8) 


Fir den Formparameter H des Nachlaufprofils wurde experimentell1 die GesetzmaBigkeit 
H=1+ konst. In 3° | ui (9) 


gefunden. Diese GesetzmaBigkeit scheint die friher verwendete Hypothese?, welche die Um- 
formung des Nachlaufs in Kérpernahe als iberwiegend von Druckgefalle beherrscht annimmt, 
wahrend die Vermischungsvorgange dagegen zuriicktreten, im wesentlichen zu bestatigen; nur 
die unmittelbare Umgebung der Hinterkante dirfte hiervon auszunehmen sein. Weit hinter 
dem Kérper wird die Impulsdicke des Nachlaufs konstant, weil dort die Druckunterschiede 
mindestens wie 1/r auf Null abklingen. Wird nun (9) in (7) eingesetzt, so bekommt man durch 
Integration von der Hinterkante (Uz) bis ins Unendliche (U..) die asymptotische Impulsdicke 
des Nachlaufs aus der Grenzschicht einer Kérperseite  _ 


’ aI, 


2 ot 5+HH 


oe (PA (10) 


~ Die gesamte Impulsdicke des Nachlaufs setzt sich aus den Beitragen der Ober- und Unter- 
seite des zylindrischen Kérpers zusammen, Die Potentialgeschwindigkeit Ug ist den beiden am 
Korperheck zusammenflieBenden Grenzschichten gemeinsam. Der Formparameter an der Hinter- 
kante betragt etwa Hy=1,8 bis 2,0, wenn die Grenzschicht dort turbulent ist. Die Werte von 
Ox sind fir beide Korperseiten nach (7) zu berechnen. 
aan Reites i: . 


+o i] 
re 


Bo Sauce Young. ‘ex. Res. Comm.-Rep. and Mem. 1838. 
“a PPh: siehe FuBnote von S. 273; vgl. auch H. B. Helmbold, Ing.-Arch. 2 (1931) S. 275. 


ae 


rf ! pene Platte ohne Anstellwinkel, bei welcher U= konst. =U, ist, geht (13a) in 
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5. Der Profilwiderstand. Bis hierher sind die verschiedenen Abschnitte einer Wider- || 
standsberechnung aus dlteren Arbeiten im wesentlichen bereits bekannt. Das Ergebnis,. das 


wir gewinnen, indem wir die Teilergebnisse (5), (7), (10) umd (1) der Reihe nach in die je- | 
weils nachste Gleichung dieser Folge einsetzen und die ubliche dimensionslose Schreibweise || 


lip = Wald Re zL (Widerstandszahl) anwenden, lautet nun 


oder 
Ux ote U,,\2+H8 Poi Tey 6 B Ch Sa, U A ae wl 
Cup = 22 (TG -) Un) ta oD coon (ze) (11) | 


u 


worin € = x/L ist. Abgesehen davon, daB dies Ergebnis nicht sehr ibersichtlich ist, weit es den: | 
Nachteil auf, daB noch die Potentialgeschwindigkeit an der Hinterkante Ug darin vorkommt; | 
es ist aber ziemlich schwierig, fir diese Gré8e begriindete Annahmen zu machen. Deshalb ist 


es ein sehr erwiinschter Umstand, daf die Exponenten 4 (5+ Hy) und 2+H zahlenmaBig immer 
so nahe beieinander liegen, daB man sie naherungsweise auch gleichsetzen kann. Dann fallt Un 
aus dem Ergebnis iberhaupt ganzlich heraus (daB das Ergebnis einer Widerstandsberechnung 
von der angenommenen Gréfe der Potentialgeschwindigkeit Ug nahezu unabhangig ist, hatte | 


A. Walz, Dt. Luftf.-Forschg., Unt. u. Mitt. 3060, bereits gelegentlich bemerkt). Mit der Annahme 
1. =e y 
as 3 (9+ Ha) =2+H a 


und mit dem bekannten Ausdruck fir die mittlere Reibungszahl einer glatten ebenen Platte. 2) 
von der Gesamtlange L bei vollaminarer Grenzschicht 


| a2 lara a nach Blasius Fe 
ergibt sich also einfach : | 
ye U U eae oe 
ce Th) Eto S (ae ) ae). (12) jf 


Dadurch, daB nun auch noch die Exponenten 1+2 H und B sia abi wenig voneinander ~ 
verschieden sind, ist eine weitere Vereinfachung des Ergebnisses méglich. Mit der Annahme | 
ey 1+2H=f re (aq | 
erhalten wir das Endergebnis ; ae | 

U 245 ae 7am 
o =z (ag Pag ae af (a ) ae) ; A. eee a) 
Dic Annahmen (1) und (II) lassen Wins durch ace Wahl der an sich nicht sehr’ genau | a 


-festliegenden Konstanten und Parameter miteinander in Einklang bringen, indem beispielsweise _ 
folgende rechnerisch bequemen Werte gewahlt werden: 


pb =4,0; Rie Fae) eg: Hu = 2,0. , K: Side ay 


z Dies nie einer mittleren Reibungszahl i im Turbulenzteil ep=0, 0040 und folglich nach Zahlen- 
_ tafel 1 Reynoldszahlen ®t von der Sere ae, 2/108 ens Byechon: Mit dieser Konstanten- 
_wahl lautet das nae ; 


Ras us 


Cw 


sie Der Wert Eq= ey |L ist bei schlanken Kérpern ones chilighy: gréoBer als’ Eins: Fir eine elanbe ‘ 


Cup = [en Ves cp (Su — -é,)] =2 [ene fut ep (1— é.)] 


ee tiber, - was mit dep letzten Gleichung des. Abselinitens 3 véllig gleichbedeutend ist. Der Vergleick § | 
von (13a) mit der Plattenformel zeigt, wie im olgcrneten, Fall des umstrémten zylindrischen a) 
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Kérpers die Widerstandszahl auBer von den Reibungszahlen cy und Umschlagstellen &, auch 
noch von der Geschwindigkeitsverteilung U(£)/U. langs der Kontur und der Geschwindigkeit 
U./U.. an der Umschlagstelle abhangt. 

Bei mittleren Reynoldszahlen ®% von der GréBenordnung 10° und darunter 
liegt die Umschlagstelle erfahrungsgemaB in der Nahe der laminaren Ablésestelle, die durch 
das Kriterium 1 ~ —0,084 festgelegt ist; in diesen Fallen fihren also die Bedingungen fir die 
Ablésung der laminaren Grenzschicht von der. Wand nicht zur Ablésung, sondern zum Um- 
schlag in den turbulenten Strémungszustand, in welchem die von der AuBenstrémung auf die 
Grenzschicht durch turbulente Mischbewegung (Impulsaustausch) ausgettbten Scheinschubspan- 
nungen stark genug sind, um die energiearme Grenzschicht gegen den Druckanstieg mitzu- 
schleppen und auf diese Weise die Ablésung zu verhindern. Die laminare Grenzschicht durch- 
lauft hier hinter dem Geschwindigkeitsmaximum (xm) bis zur Umschlagstelle (x) noch einen 
Verzégerungsbereich (Druckanstieg) mit A<0, in welchem zwar die in Abschnitt 2 benutzte 
Annahme iiber die Einparametrigkeit der Grenzschichtprofile nicht mehr genau zutrifft, die 
bestehenden Zusammenhange aber doch noch wenigstens im Groben angenahert werden kénnen, 
indem die zweite Konstante der Gleichung (3a) mit B=5,0 angenommen wird. Man bekommt 
dann an Stelle von (4) fir x> xm 

2 *m a 
eas =n (Umax f U*%dx+ f U5 dx), (4a) 


0 “8 
m 


v 


woraus insbesondere fiir x=xy>%m als Ablése- bzw. Umschlagkriterium zur Bestimmung der 


Umschlagstelle &, 


BO tates 1272 Ue, 40 f U \ 5.0 
da = a Cae) ee (oe) der ti Ga) ae) ~ — 0,084 (14) 
und an Stelle von (6) 
Ou reer Su U ,4,0 g Su U ,5,0 y 
fa Vat (Se Sh) dE+ f(y) 48) (64) 


folgt. Dies wird unter Beriicksichtigung der Konstantenwahl (III) in die Vorform der Gleichung (11) 


ap eS 
2 Dice ry Uy, \3:5 [ oy Whot3s3 


eingesetzt und ergibt nach kurzer Zwischenrechnung © 


U (Umx. Je, U \40 SPU. 50 SU Saha. ee 
cm =F (en C2 PE) Cee uw.) at) + eH f (qa) ae). (15) 


u 


Die in dieser Arbeit gemachten Angaben iiber die Lage der Umschlagstelle entsprechen dem 
gegenwartigen fragmentarischen Stand unserer Kenntnisse tber den Turbulenzbeginn. Ihren 
vollen Nutzen werden unsere neuen Formeln also erst dann bringen kénnen, wenn man hieriiber 
einmal besser unterrichtet sein wird. 


6. Zusammenfassung. Durch geeignete Abstimmung der in den Grenzschichtrechnungen 
~ auftretenden empirischen und theoretischen Konstanten aft sich die schlecht abzuschatzende 
Potentialgeschwindigkeit am Kérperheck aus der Widerstandsberechnung ganzlich entfernen 
und eine sehr tbersichtliche Widerstandsformel gewinnen. 


; (Eingegangen am 11, August 1948.) 
Anschrift des Verfassers: Professor Heinrich Bruno Helmbold, (13b) Miinchen 23, RémerstraBe 14 IT. 
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Uber den durch einen Verdichtungssto8 verursachten Zusatzauftrieb. 
Von H. B. Helmbold. 


1. Einleitung. Der Auftrieb eines unter kleinem Anstellwinkel « angeblasenen schlanken 
zylindrischen Kérpers in einer ebenen Unterschallstrémung folgt der Prandtlschen Regel, wo- 
nach bei festgehaltenem Anstellwinkel die Auftriebszahl ca = A/q.F' dem Prandtlschen Faktor 


R= 1/V 1—92, proportional ist: ca = Ba. [2 1-6 f(H)] - (1) 


Die Bedeutung der Formelzeichen ist folgende: A Auftrieb, q.5 =Q. U?/2 Staudruck mit Oc 
als ungestérter Dichte und U als Anblasgeschwindigkeit, F gréBte Projektionsflache des zylindri- 
schen Kérpers, 9. = U/a. Machzahl mit a, als ungestorter Schallgeschwindigkeit, 6 =d/L 
Dickenverhaltnis mit d als gréBter Dicke und L als gréBter Lange im Profil des zylindrischen 
Korpers, &% = U L/y Reynoldszahl mit y als kinematischer Zahigkeit. 

Bei hinreichend groBer Reynoldszahl verschwindet. die Abhangigkeit von derselben, also 
wird f()=konst. Zum Bau der experimentell! gesicherten Gleichung (1) ist von seiten der 
Theorie zu bemerken, da®B in der kompressiblen Potentialstrémung die Auftriebszah!l cq mit 
der Machzahl t,. zunehmend starker als proportional $$ wachst, doch besteht daneben in der 
wirklichen Strémung infolge der inneren Reibung ein den Auftrieb schwachender Grenzschicht- 
einfluB, welcher seinerseits der Prandtlschen Regel gehorcht und den UberschuB der theoretischen 
Auftriebszahl tiber 2 2a hinaus sogar tberwiegt; so kommt eine zufallig fast lineare Beziehung 
nach (1) zustande?. | Ces 

Die der Gleichung (1) zugrundeliegende Prandtlsche Regel setzt auBer fast paralleler Stré- 
mung (schlankes Profil bei kleinem Anstellwinkel) streng genommen noch voraus, daf die Stré- 
mung reine (unterkritische) Unterschallstrémung ist, daB also die Stromungsgeschwindigkeit 
nirgends die kritische Schallgeschwindigkeit erreicht oder tiberschreitet. Diejenige Machzahl M.., 
bei der gerade in (mindestens) einem Punkte der Strémung die kritische Schallgeschwindigkeit 


_erreicht ist, hei®t die kritische Machzah! tyr. Es ware demnach zu erwarten, daB die Abhangig- — | 


keit der Auftriebszahl von der Machzahl nach (1) oberhalb Nt ihre Giltigkeit verliert. Die Ex- 
“ perimente zeigen jedoch einwandfrei, da dies im allgemeinen nicht der Fall ist: bei gréeren 
Anstellwinkeln wurden an diinnen symmetrischen Profilen mit 40—42 %. Dickenriicklage be- 
trachtliche Uberschreitungen der kritischen Machzahl (um 56 % bei 6=0,0875 und «=4° nach 
H, B. Helmbold, um 47% bei 6=0,10 und %«=6° nach O. Knappe, um 40 % bei 6=0,12 und 
a=7° nach B. Géthert) ohne Abweichung der Auftriebszahl von der linearen GesetzmaBigkeit 
- nach (1) beobachtet. Nur bei kleinen Anstellwinkeln erfolgt eine Anderung im Auftriebsverhalten 
sofort oder kurz nach Uberschreiten der kritischen Machzahl. Diejenige Machzahl, von der ab 
die Auftriebszahl nicht mehr der Gleichung (1) gehorcht, wollen wir t4: nennen. Die Anderung 
im Auftriebsverhalten bei Uberschreiten von Maz besteht im plétzlichen Auftreten eines Zusatz- 
auftriebs zum Auftrieb nach (1), der etwa der empirischen GesetzmaBigkeit Aca 0 (Mo —Ntar)? 
gehorcht. Nicht unerwahnt darf bleiben, daB der Zusatzauftrieb seiner GréBe nach auBer vom 
Dickenverhaltnis 6 sehr stark durch die Nachbarschaft fester Wande oder freier Strahlgrenzen 
beeinfluBt wird; er ist unter tblichen Verhaltnissen im Kanal mit. festen Wanden mehrfach’ ~ 
(Gré8enordnung: zehnmal) gréfer als im Freistrahl. Das plétzliche Auftreten des Zusatzauftriebs. 
nach unter Umstanden betrachtlichem Uherschreiten der kritischen Machzahl wurde sofort 
mit dem gleichzeitigen Auftreten eines VerdichtungsstoBes in Verbindung gebracht, die Mach- 
zahl Mar mithin als ,,StoBmachzahl‘ (Ms:) gedeutet.. Der experimentelle Beweis konnte in | 
schlissiger Form jedoch erst 19453 durch die Auswertung der im FIA T-Report mitveréffentlichten _ 
Géthertschen Versuche erbracht werden. Hier soll nun der Beweis fir den aus den Versuchen 
zunachst nur erschlossenen, dann auch nachgewiesenen Zusammenhang zwischen Zusatzauftrieb 
und VerdichtungsstoB auf theoretischer Grundlage gefihrt werden, indem gezeigt wird, inwie- 
fern ein Verdichtungssto$ einen Zusatzauftrieb iiber den bei stoffreier Strémung bereits vor- 
handenen Auftrieb hinaus verursacht. etal 2 
1 H. B. Helmbold, Jahrb. 1942/43 d. Dt. Akad. d. Luftfahrtforschg., S. 367; O. Knappe, Ber. 156 d. Lilien- . 
thal-Ges. f. Luftfahrtforschg. (1943), S. 128; H. B. Helmbold, FIAT-Rep. 1092 (1947), - i 
* A. Busemann, Schriften d. Dt. Akad. d. Luftfahrtforschg., Heft 18 (1940), Diskussionsbeitrag S.22.; . 
H. B. Helmbold, Dt. Luftfahrtforschg., Unt. u. Mitt, 2502 (1944), : . , 
* H. B. Helmbold, Dt. Luftfahrtforschg., Unt. u. Mitt. 2505 (1945). 
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Zur Frage des. ,,StoBverzugs“, d. h. der Tatsache, daB ein StoB unter Umstanden erst nach 
sehr erheblicher Uberschreitung der kritischen Machzahl auftritt, sei hier nur auf die Stellung 
des Problems durch G. I. Taylor! und die Entdeckung mathematischer ,,Grenzlinien® fur die 
adiabatische Potentialstrémung durch W. Tollmien2 hingewiesen. Wenn gezeigt werden kann, 
dafs das Vorhandensein eines Verdichtungsstofes (auBer im Fall des symmetrisch angeblasenen 
symmetrischen Profils) einen Zusatzauftrieb bedingt, kann der Nachweis des Zusatzauftriebs 


den, unmittelbaren Nachweis des VerdichtungsstoBes und gegebenenfalls auch des StoBverzugs 
ersetzen. 


2. Die Methode der Expansionsquellen. Zur Darstellung der kompressiblen Potentialstrémung 
verwenden wir die von L. Poggi’ begriindete Methode der Expansionsquellen. Diese beruht darauf, 
da es méglich ist, die kompressible Gasstrémung kinematisch als Strémung einer inkompres- 
siblen Flissigkeit aufzufassen, welche in Beschleunigungsgebieten zusatzliche Flissigkeit her- 
vorbringt und in Verzégerungsgebieten Flissigkeit verschwinden laBt. Da mit der Beschleu- 
nigung des strémenden Gases Expansion, mit seiner Verzogerung Kompression verbunden ist, 
sprechen wir von Expansionsquellen und Kompressionssenken. Diese ortsfesten Quellen und Sen- 
ken, welche dazu dienen, die Dichteveranderungen der kompressiblen Gasstromung nachzuahmen, 
haben jedoch keine dynamische Bedeutung, sondern nur den Sinn eines kinematischen Ge- 
dankenmodells. Der VerdichtungsstoB wird sich uns als flachenhafte Kompressionssenke dar- 

' stellen; einen VorstoB in ahnlichem Sinne hatte bereits O. Knappe* unternommen. : 

Auf anschauliche Weise ergibt sich die vorstehend skizzierte Vorstellung folgendermaBen: 
Fur eine Stromréhre in der kompressiblen Gasstromung, deren als sehr klein-angenommener . ed 
-Querschnitt f langs der Stromlinie:s veranderlich ist, besagt die Kontinuitatsbedingung, daB | Oe 
durch jeden solchen Querschnitt in der Zeiteinheit die gleiche Masse hindurchstrémt; der Massen- Ne 
flu8 ist also ; I 

onic m=owf=konst., (2) 


wobei @ die Dichte und w die Geschwindigkeit bedeutet. Das durch die verschiedenen Quer- 
_schnitte in der Zeiteinheit hindurchstrémende Flissigkeitsvolumen ist jedoch verschieden, nam- 
lich an der Stelle s=sy der Stromlinie V}=w,)f, und an der stromabwarts darauffolgenden Stelle 


/s=s)+ds) : = 


- 


naa ACT weet Ha ae 

Coad ds ) Une } , : ~o ie e 

_ Also beiudet sich zwischen den bei sy und s)+ds gelegenen Querschnitten eine Volumen- er ea ; 

quelle von der Ergiebigkeit (Volumen je Zeiteinheit) : : . 

: -  d(wf) ,.__ d(mn/g) ido, wf do eh. 
daAVv= ga dee oe = o ds ds = 0 ag ds. (3) a 


-Indem wir nun aus der Bewegungsgleichung dp/o+wdw=0 und aus der Gleichung fir 
die Fortpflanzungsgeschwindigkeit schwacher Druckstérungen, die Schallgeschwindigkeit a, | 
dp/dg=a* den Druck p eliminieren, erhalten wir 
ee My: (4) 
< : e Le 
Dies in (3) eingesetzt ergibt zunachst : 
ae | ay a (\* SY peas 
ze pe Nee ds 


Nes / ——— 
Fay — 


‘und, wenn wir fortan die Quellergiebigkeit mit Q statt V, das Volumen der Stromréhre zwischen 
den bei s, und s)+ds gelegenen Querschnitten mit dr=f ds bezeichnen und beachten, dal im 
_ Ganzen der Stromung betrachtet w auch noch von den zu s orthogonalen Koordinaten abhangt, 


dx \a ds 


By 1, I. Taylor, 5. Convegno Volta 1935, S. 198, R. Ace. d’Italia, Roma 1936. 
2 W. Tollmien, Z. angew. Math. Mech. 17 (1937) S. 117; 21 (1941) S. 140. — G. Guderley, Z. angew. Math. 
Mech. 20 (1940) S. 185. os “ wet? | 
3 LZ. Poggi; Aerotecnica 12 (1932) 5.1579; 14 (1934) 5. 532. 


y : 
“Say O. Knappe, Ber. 156 d. Lilienthal-Ges. f. ‘Luftfabrtforschg. (1943), S. 174. 
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als Ausdruck fir die Dichte der Volumenquellen in der kompressiblen Strémung. Hieraus ist 
unmittelbar abzulesen, da® der Beschleunigung (dw/ds> 0) Expansionsquellen, der Verzégerung 
(dw/ds<0) Kompressionssenken entsprechen und da deren Dichte dem Quadrat der értlichen 
Machzahl 9% = w/a proportional ist. ; 

Ist die raumliche Verteilung der Quelldichte nach Gl. (5) bekannt, so kann die kompressible 
Strémung als inkompressible Potentialstromung berechnet werden. Da aber auf der rechten 
Seite. dieser Gleichung die gesuchten Geschwindigkeiten vorkommen, mu man iber eine Aus- 
gangslésung verfiigen, von der aus verbesserte Lésungen im Iterationsverfahren gewonnen wer- 
den kénnen. L. Poggi benutzte die quellenfreie inkompressible Strémung als Ausgangslésung, 
in neuerer Zeit verwendet man hierfir meist die Prandtlsche erste Naherung fir die fast parallele 
Unterschallstrémung in der von A. Busemann angegebenen Form!; wir werden uns diesem Ver- 
fahren anschlieBen. 


3. Die kompressible Str6mung mit Verdichtungssto8. Fir die Anwendung auf den Verdich- 
Paget tungssto& mége der einfacheren Darstellung halber angenommen werden, daf dieser die Strom- 
= . linien tiberall senkrecht durchsetzt (orthogonaler Sto), was jedenfalls mindestens angenahert 
eee zutrifft, solange die Strémung am Koérper anliegt; der StoB soll also schwach genug sein, um 
ae keine Ablésung der Grenzschicht (falls wir eine solche iberhaupt als vorhanden beriicksichtigen 
on wollen) hervorzurufen; andernfalls ergibt sich ein krummer oder schrager Sto%. Ubrigens ist 
die Annahme auf das Ergebnis ohne Einflu8. Der StoB fallt also innerhalb der betrachteten 
Stromréhre mit einem ihrer Querschnitte zusammen. Durch die Vorderflache f, des StoBes 
strémt in der Zeiteinheit das Volumen V, =w,f,, durch seine Rickenfliche f, das Volumen je 
Zeiteinheit V. =Wefo, die Flachen f, und f, fallen aber raumlich zusammen, so daB ff=fp=fs 
gilt. Da der Massenflu8 durch alle Querschnitte konstant gleich gwf ist, entspricht der sprung- 
haften Verdichtung von Q, auf 0, > @, bei konstantem Querschnitt fs eine sprunghafte Abnahme 
der Geschwindigkeit von w, auf w.<w, und eine ebensolche Abnahme des Volumenflusses © 
um den Betrag (w,—w,)fs:. Der StoB bedeutet also eine flachenhafte Senke von der (negativen) 
Ergiebigkeit ; 

f 


Qst = (w,—w,) for << 0. (6) 


_ Wir missen nun im Hinblick auf spater folgende Uberlegungen noch zeigen, daB der StoB 
eine flachenhafte Konzentration von Kompressionssenken darstellt, die bei fehlendem Sto8 
raumlich verteilt vorhanden zu denken waren. Mit oder ohne Sto8 gilt namlich der Satz: 

, Die Summe (Integral) ailer in der kompressiblen Strémung vorhandenen Expansionsquellen 
ist innerhalb jeder einzelnen Stromréhre und daher insgesamt gleich Null. mi) 


Dies folgt durch Integration der Gleichung dQ =dV = d(wf) bis zum StoB: 


$y 
dQ 5 
W, fs = J a4 woe fos 
So 


_— : d . 
Sere We fst ae 1 aa ds+ Wo foo roa lie 


_hieraus und aus (6) ergibt sich 


(7) 


, y . y sil 7 s2 > 
a —0s: = (wy — We) fst = { ae as = f oe ads:3 
salso oo): spe ae a. oe 2 7 
P : sy eo} ¢ 
dQ f dO aap Dees. . 
J GedstOut f Zpas=o, f ; (8) 
_ So x 


was zu beweisen war. Die Umwandlung der Strémung ohne Stof in eine solche mit Sto® kénnen 

wir gedanklich aber auch auf die Weise vornehmen, da8 wir mit der Einfiigung des StoBes in 

die stoBfreie Strémung gleichzeitig raumlich verteilte Expansionsquellen hinzufiigen; im Grunde 

ist das mit der vorigen Betrachtungsweise gleichbedeutend, denn das Fortnehmen einer Senke 
ist dem Hinzufigen einer entgegengesetzt gleichstarken Quelle gleichwertig. 


1 A, Bustmann, Jahtb, d. ‘Wiss. Ges. f. Luftfahrt 1928, S, 95: Artikel uae ae E 
-perimentalphysik, Bd. IV, 1, S. 410. Leipzig 1931. 9; Artikel ,,Gasdynamik“ im Handb. d. Ex- 
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Bei der Bestimmung der Integrationskonstanten w,, f.. in (7) haben wir stillschweigend eine 
Vereinfachung vorgenommen, deren Zulassigkeit durchaus noch erértert werden mu8. Ist nam- 
lich ein StoB vorhanden, so erfahrt beim Durchgang durch denselben der Ruhedruck (Pitot- 
druck) py eine Verringerung. Da die Strémung hinter dem StoB wieder als adiabatisch gelten darf 
(Warmeleitung wird vernachlassigt), ist die Dichte im Unendlichen hinter dem Sto kleiner als 


‘im Unendlichen davor 
x—1 
[ese en (Ba) 


und nach der Kontinuitatsbedingung ist dann (wf )+0>(wf)_.- Die Vereinfachung, welche 
(wf) +0 = (wf) o = Wo fo sétzt, ist nur dann erlaubt, wenn der Ruhedruckverlust im Stof ver- 
nachlassighar klein ist, wenn also der Sto8 so schwach ist, daB die Zustandsanderung des Gases 
im Sto8 noch als adiabatisch gelten darf. Wir miissen uns demnach, da wir an der Giltigkeit des 
soeben aufgestellten Satzes festhalten wollen, auf schwache Sté8e beschranken; wir dirfen es, 
weil starkere StéBe als die der Einschrankung unterworfenen die Grenzschicht der Strémung 
von der Kérperoberflache veranlassen wiirden, womit die Berechnungsméglichkeit ohnehin 
verloren geht. Wir gewinnen dabei weiter den Vorteil, da® wir das Vorhandensein von Wirbeln 
im Nachlauf des StoBes, d. h. auf den Stromlinien, die den Sto8 durchsetzen, auBer acht lassen 
und die Strémung weiterhin als adiabatische Potentialstromung behandeln dirfen; nur dirfen 
wir sie nicht als isentropisch auffassen, weil die Entropievermehrung im Sto8 zwar sehr klein, 
aber doch noch von Null verschieden ist ; wir sprechen also von einer nichtumkehrbaren Potential- 
strémung, da sie im StoB selber nicht umkehrbar ist. Ein .» VerdiinnungsstoB" ist ja wegen der 
_ mit ihm verbundenen Entropieverminderung physikalisch unmdglich. 


4. Der Zusatzauftrieb des Verdichtungssto8es. Vorbereitend bestimmen wir nun die Auf- 
triebsanderung, welche entsteht, wenn man in die bakomiprercile Strémung um den zylindrischen 
_ Kérper einen Quellpunkt einfiigt. ses : le 
Indem man sich die Expansionsquellen und Ronipaloianssraken als esha Punkte uber 
den Raum verteilt bzw. in der StoBflache konzentriert vorstellt, ergibt sich der Vorteil, daB sich 
ihr Feld der quellenfreien inkompressiblen Grundstrémung linear iberlagert und da® fir die 
uns hier interessierende ebene Strémung die hydrodynamischen Singularitaten bei einer kon- 
_ formen Abbildung erhalten bleiben. Die abgebildete Stromung stellt dann allerdings keine mog- 
_liche Gasstrémung mehr. dar, weil sie ja mit der Abbildung gegeniiber der Originalstrémung ix im 
-unendlich Kleinen eine Volumenverzerrung gema} 
a apiece af : dv’ dtp 
~erfahren hat und demnach der ‘Gieiehuiiy (5) nicht mehr gehorchen wird. In jedem Falle, bei 
der Originalstrémung wie bei der abgebildeten Strémung, hat man noch zu beachten, daB die 
‘Hinzufigung— von Quellen zur inkompressiblen Grundstrémung um die Kontur die Rand- 
bedingung verletzt, solange man nicht im Innern der Kontur passende weitere Singularitaten 
_ hinzufigt, welche die Randbedingung wiederherstellen. Um dies in einfacher Weise zu bewerk- 
_ stelligen, nehmen wir die konforme Abbildung der gegebenen Originalkoatur auf den Einheits- 
p rele in dessen Mittelpunkt der Koordinatenursprung liegen mége, zu Hilfe. 

In der Kreisebene (¢) befinde sich auBerhalb des Kreises im Punkte P,(&.7) ein Guelipunkt 
von der Ergiebigkeit Q; dieser verursacht gemeinsam mit den im Kreisinnern anzubringenden 
et aeclion im Punkte P(é,7) der ae oatnk eine Storgeschwindigkeit mit den Komponenten . 
(uw |[é, 0" y ss a, ; 

ba : 0 GViEF np IRF | 


CR a by, eS ah 
Perey A tee Beet m MI-P)— 048 i 
re - Oe bee rer 
ee ee eae WG ebm VI-@)—A+R9_ 


4 Poles ReaB by oe Die “Abbildung eee wir so, daB derjenige Punkt des eee 

der dem ZusammenfluBpunkt (Hinterkante) der gegebenen Kontur entspricht, die Borer tey 

i 1, J =0 hat; dann Eden ek aus (9) fir, diesen Punkt “if 
Buta gael ac urea panne ee Oita gin ail (10 

rt a a a ‘ a 0 cae v v(o= ee “3 URE (1&0? ’ * (1 ) 
ee * > c A : f ¥ hi vi ; 18 
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also gerade doppelt so viel, wie wenn der Quellpunkt ohne den umstroémten Einheitskreis allein 
vorhanden ware. Da dieser Punkt (C=1) jedoch ebenfalls ZusammenfluBpunkt sein muB, mu 
dort die Geschwindigkeit iberhaupt verschwinden, was ohne Stérung der Randbedingung nur 
dadurch méglich ist, daB ein im Kreismittelpunkt gelegener Wirbelpunkt hinzugefigt wird, 
welcher im Punkte €=1 die Geschwindigkeit A['/2 7=v' zusatzlich verursacht; er mu also 
nach (10) die Wirbelstarke (rechtsdrehend positiv) 
m1 
Ana ~ 2 OTE) 

haben. Nachdem die Koordinaten des Zusammenflu8punktes in der Kreisebene mit =1, 7=0 
festgesetzt wurden, mu die quellenfreie Grundstromung in dieser Ebene aus einer Parallel- 
stromung W’ mit den Komponenten U’=W’ cos a||£ und V’= W’ Sin a||7, einer entgegengesetzt 
zu W’ orientierten Dipolstromung mit dem Moment M=2 7U’ und einer Wirbelstr6mung mit | 
der Zirkulation [)=42V’ (Dipol und Wirbel liegen im Kreismittelpunkt) zusammengesetzt 
sein, damit auch in dieser Grundstrémung die Geschwindigkeit im Punkte =1, 7=0 ver- 
schwindet. Dementsprechend werden wir auch in der Originalstrémung um die Profilkontur 
das Koordinatensystem so wahlen, daf die Bezugslange der Kontur (Profiltiefe) L gleich dem 
Durchmesser D=2 des Einheitskreises ist, die Hinterkante die Koordinaten x=1, y=0 erhalt 
und U=W cos «||x sowie V= Wsin a|ly ist. Nach (11) verursacht der Quellpunkt in der Stré- 
mung (W’) um den Einheitskreis die zusatzliche Auftriebszahl 

, DAE 4n, Q : : 
| Aca = TT = > G80? LW’ (2) 
und der entsprechende Quellpunkt in der Originalstrémung (W) beispielsweise um eine Ellipse | 
vom Dickenverhaltnis 6 die zusatzliche Auftriebszahl 4| 


(11) 


§ 


gl UE An, Oxy 
« dees LW (SS) EW * th 
OE aseiarea dabei ist Aca um einen Faktor W’/W=4(1+-6) kleiner als 4c, ; denn um L (die groBe Achse-der 
ee Ellipse) gleich dem Durchmesser D=2 des Kreises zu machen, muf die Abbildungsfunktion 
Bi: _ 146 1-6 1 , 
hee Reg Sa aan ener: 
ace gewahlt werden, und dann wird im Unendlichen (= 0) 
[dz] _ |W] _ 1+8 
Rtg Ve, UW AS oh ae 
a mG Nachdem wir uns nun schon klar gemacht haben, daB der Verdichtungssto8 kinematisch als 
i Kompressionssenke aufzufassen ist, kénnen wir) das vorlaufige Ergebnis (13) als den Zusatz- >|} 


auftrieb deuten, den ein Verdichtungssto8 mit der negativen Ergiebigkeit Q unmittelbar, d. h. 
ohne Beriicksichtigung der dem Stof zugehérigen raumlich verteilten Expansionsquellen, ver- 
ursachen wurde, wenn wir die ber die StoBflache (in der ebenen Strémung: Stoflinie) verteilte 
gesamte Ergiebigkeit vereinfachend in einem Punkte, beispielsweise dem Schwerpunkt der 
Intensitatsverteilung mit den Koordinaten x,,y, konzentriert annehmen; €, und 7, sind dann 
die Koordinaten des abgehildeten Schwerpunktes in der Kreisebene. Natirlich ist das, auch 
von der Vereinfachung abgesehen, noch nicht der wahre Zusatzauftrieb. Diesen wirden wir 
erst erhalten, wenn wir durch Integration einer der Gleichung (13) entsprechenden Differential- |} 
formel den Einflu8 der raumlich verteilten stoBbedingten Expansionsquellen hinzunehmen 
wurden, welcher entgegengesetztes Vorzeichen. hat und somit den Zusatzauftrieb gegentiber 
unsrer vorlaufigen Abschatzung vermindern muff. Die genaue Ermittlung der zusatzlichen | 
Quellverteilung ist jedoch recht mithsam, so da® wir nach einem Verfahren suchen miissen 
welches den erwahnten Einflu8 wenigstens naherungsweise zu erfassen gestattet. . ; 
Hierfir wollen wir uns der von vornherein gemachten Voraussetzung erinnern, daB der um- 
stromte zylindrische Kérper ein schlankes Profil hat und unter kleinem Anstellwinkel angeblasen 
wird, so dafs die Strémung nur wenig von einer Parallelstromung abweicht. Dann wird das 
Geschwindigkeitspotential ® aus einem Grundpotential ® = Ux und einem kleinen Stor- 
potential y bestehen und die Stérgeschwindigkeiten u=dq/dx und v=09/0y werden ebenfalls 
klein gegen U=0@,/0x sein. | Die Kontinuitatsbedingung der kompressiblen Strémung nimmt | 
bei Verwendung kartesischer Koordinaten die bekannte Form. | . i 
d[e(U+u v) ; a) | 


0x ae 
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an, welche besagt da die Dichte der Massenquell i i i i 
quellen (das ist die Bedeutung der linken Seit 
gleich Null ist. Diese Gleichung formen wir in oa - oe 
‘du OV: fea 00 v a2) 

0x 0x e dere O Oxy 

um. Die linke Seite der neuen Gleichung ist dann aber nichts anderes als die Dichte der Volumen- 
quellen, die wir unter der Bezeichnung dQ/dt mit (5) bereits berechnet haben; nur geht deren 


rechte Seite fiir die jetzt betrachtete fast parallele Strémung in (U/a,,)?- du/d« tber, so daB 
wir (5) durch 


(14a) 


du O07 i Ne eat ig Ot 
ve tae go) Ge =e (15) 
oder nach Einfiihrung des Stérpotentials durch ic 
0? 0° 
(1—me,) = + =0 (15a) 


: 0 x? 0¥7 
ersetzen kénnen. Diese Jinearisierte Potentialgleichung der kompressiblen Strémung geht fir 
Mien <= Poin die Potentialgleichung einer inkompressiblen Strémung (i) tiber, wenn man nach 
L. Prandtlt die Koordinaten entsprechender Punkte gema® den Gleichungen xi = x und 


yi = V1—M, y verknipft. Nach A. Busemann! miissen wir die Stérpotentiale in entsprechenden 
_Punkten durch gi =(1—2,)y zueinander in Beziehung setzen, damit die Randstromlinien 
(Kontur) der kompressiblen Strémung in Randstromlinien der inkompressiblen Strémung iiber- 
gehen (,,Stromlinienanalogie‘‘), was sich leicht verifizieren laBt. Man erhalt demnach umgekehrt 
die Stérpotentiale pm in einer kompressiblen fast parallelen Unterschallstrémung, wenn man die 
Stérpotentiale gy; in einer inkompressiblen Vergleichsstrémung, deren Querkoordinaten yi 
gegenuber y um den Faktor 1/$ kontrahiert sind, mit dem Faktor $2 vergréBert. Fir die Stér- 
geschwindigkeiten in entsprechenden Punkten gilt dann u=$2u; und v= Bu. 
3 Wesentlich fir unseren Zweck ist, daB die Gleichungen (15) und (15a) in den Storgeschwindig- 
_ keiten und im Stérpotential linear sind, so da®B man ihre Lésungen tberlagern darf. Es sei nun 
in der kompressiblen Strémung langs des Stiickes h einer Aquipotentiallinie ein Verdichtungs- 
sto} mit dem (langs h verdnderlichen) Geschwindigkeitssprung w,—w, vorhanden. Ihm ent- 
spricht in der inkompressiblen Vergleichsstrémung nach der Stromlinienanalogie eine tiber- 
lagerte Senkenverteilung mit dem Geschwindigkeitssprung w;, —wj,=(w,—w,)/$? langs des — 
Stiickes h;=h/ der entsprechenden Aquipotentiallinie. Die (negative) Ergiebigkeit (Qi)s: dieser 
Senkenverteilung ist gegentiber Qs; um 1/%% verkleinert, weil die Ergiebigkeit preportional 
dem Produkt aus dem Geschwindigkeitssprung und der Stofflache (StoBlange h) ist. Nach (13) 
-verursacht dann_ die ‘Senkenverteilung -(Qi)s¢ in der inkompressiblen Strémung eine: Zusatz- 


auftrieb | ee: . ih Ose és) 

ei < : B (l—E)*+n" LW 
_wobei W=U ist; dieser Ausdruck ist positiv, weil Qs; negativ ist. Dem Vorhandensein der 
_ Senkenverteilung (Qi)s entspricht .ein gewisser Beitrag Au; zu den Stérgeschwindigkeiten ui 
und zu der in (15) rechts vorkommenden Ableitung dui/dx, aber kein Beitrag zur Dichte der 
Volumenquellen fae genet) oe Gt erakte ie ae 
om 2  aee GAT ae 
weil auf der rechten Seite von (15) bei inkompressibler Strémung M., =0 zu setzen ist. Anders 
‘ist dies nun in der nach der Stromlinienanalogie zugeordneten kompressiblen Originalstrémung: 
dort sind aufer der StoBsenke Qs; nach (15) zusatzliche stoBbedingte Quellen 
a _ (AQ) _ gpa ye 94 0 Pia? 

: Pe ay Mean se ax 

vorhanden, welche den Zusatzauftrieb abschwachen. Ihren Einflu® erfassen wir, indem wir 
den dem imkompressiblen Zusatzauftrieb A ca; der Gleichung (16) entsprechenden kompressiblen 
Zusatzauftrieb Aca nach der Stromlinienanalogie berechnen. Die Auftriebszahlen sind den 
Stérpotentialen proportional, beim Rickgang von der inkompressiblen Vergleichsstrémung zur 
kompressiblen Originalstroémung also mit $$? zu vergréBern. So erhalt man entsprechend (16) 
den Zusatzauftrieb des VerdichtungsstoBes ay ‘arte 


’ 


; Ata = 


aa ony wy—w, dn 17 
ip . (1)? +0? ee setade | ; a 


1 FuBnote 1 Seite 282. — 


; hy <i a 
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Der Faktor 1/9 stellt hier den abschwachenden Einflu8 der stoBbedingten Expansionsquellen 
dar. Auer den stoBbedingten Expansionsquellen sind natirlich von vornherein andere Ex- 
pansionsquellen nach (15) bereits vorhanden; ihr EinfluB ist aber in dem ,,stoBfreien“’ Auftrieb 
: nach (1) schon enthalten. Mit n ist der langs des StoBes (Aquipotentiallinie) gemessene Wand- 
abstand bezeichnet; der StoB beginnt: ja bekanntlich an der Kérperoberflache, wenn das Vor- 
handensein der Grenzschicht unbericksichtigt bleibt. Jedem Werte n ist ein Punkt (%4,7;) in 
der Originalebene und durch die konforme Abbildung ein Punkt (£,, 7,) in der Kreisebene zu-_ | 
| geordnet. Vine 


Profil: 8,75 Gii 4475-0 
Anstellwinkel o=3° 
Druckpunktage 


‘ 5 ‘ iJ 
a ele ty RRR NE Ta le ea Ni mee . 


. 


yy ea geeragae 


2 ath) Uber dem Prandilfaktor $B = 1/ YI Wa, « 


es UES 


4 

: 
ScHoriae ns. Da sich aber die i 
, i, epee act iat , + hs 
‘\ i F sitapcts : ae : ae | 
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_Integralaussagen der auf die Stromlinienanalogie begriindeten Prandtlschen Regel nicht nur 
fir den Auftrieb, sondern auch fir die ibrigen Stromungskrafte empirisch sogar im uberkriti- 
schen Bereich (Qt. > Mr), d.h. bei Vorhandensein eines drtlichen Uberschallgebietes, als an- 
wendbar erwiesen haben, ist auch hier zu erwarten, da®B der entstandene Fehler klein, wenn nicht 
uberhaupt zu vernachlassigen ist. Dafir spricht folgender Vergleich zwischen theoretischer Ab- 
schatzung und Versuchsergebnis (Abb. 1): 

An einem symmetrischen Profil von 144mm Tiefe mit dem Dickenverhaltnis 6=0,0875 
und der Dickenriicklage £;=0,418 (in Profiltiefen L gemessen) wurden beim Anstellwinkel ~=3° 
im 360mm hohen Kanal nachstehend aufgefiihrte und mit der Berechnung (th) verglichene 
Werte gemessen:! - 


ee es a 


Cao ; (4 Ca)th 
$ Meo - Gl. (1) | Aca | Gl. (17) | fa th | 
ee eer gee en ae ee De ee AE pence yp a Se 
1,400 , 0,700 . 0,368 0,368 0 0 0,368 Ca = Can + Aa 
4 1,433 0,716 0,392 0,376 0,016 0,006 0,382 
1,474 0,735 0,459 0,386 0,073 0,045 0,435 Ca th = Cap + (A Ca)th 
1,538 0,760 0,587 © 0,404 0,183 0,141 0,545 


Wegen der Nahe der Kanalwande liegen die gemessenen Werte ca gegeniiber den Auftriebs- 
zahlen in freier Strémung systematisch zu hoch; die aus den MeBwerten abgeleiteten Werte A ca 
sind im Hinblick auf analoge Freistrahlmessungen schatzungsweise mit einem Faktor der GroBen- 
 ordnung 0,55 zu korrigieren, den ‘abschwachenden Einflu8 der Grenzschicht auf (Aca) kann 
man auf eine GroSenordnung von 10 % veranschlagen; dann besteht angenaherte Ubereinstim- 
mung bei den beiden héchsten Machzahlen. 

Bei M.. =0,735 beginnt ein plétzlicher Wiederstandsanstieg (Abb. 1), der StoB befindet sich 
gerade an der dicksten Stelle des Profils (€¢=0,418) und-an dieser Stelle setzt Riickstrémen der | 

~Grenzschicht ein; oberhalb dieser Machzahl liegt die Stromung nicht mehr vollstandig an, was 

doch Voraussetzung fir die Berechnung von (Aca)m war; trotzdem bleibt die Rechnung bis 
M.=0.75 brauchbar. Bei M,.=0,760 wurde Riickstrémen der Grenzschicht an der Mefstelle 
£=0,70 beobachtet, wahrend der Sto erst bis £+0,62 vorgeriickt ist; bereits bei M,.=0,749 
hat die Grenzschichtablésung die Hinterkante erreicht (aus den statischen Driicken der Druck- 

_ seite zu erschlieBen), damit beginnt die vollstandige Ablésung der Strémung von der Saugseite, ... 
und von M,.=0,762 ab fallt der Auftrieb wieder ab. (Abb. 1). 

Der Berechnung von (A ca)m liegen zugrunde die Riicklage & des StoBes gegeniiber der Vorder- 
kante in Profiltiefen L gemessen, der Geschwindigkeitssprung (w,—w,)) des StoBes an der Kontur 
und die StoBlange nach der Formel — Gite a « be 

f i : 2 ; a 1Jo ker aa s 

: | pa eR, E 
elcher Ry den Krimmungsradius der -Kontur und a die kritische Schall- 


: (G. Richter), in ~ 


“ geschwindigkeit bezeichnet. Aus den Messungen ergab sich: 
Bi : ary eta ean IE 


{ 


eae a | (y= e)o he 
| Pee oe iG 
Pp tin,) ex 0,700. | 0,11. 0 es 
rts oars — 0,716 0,20 _ 0,130 0,121 
Sete mafOes ea tid 4 0,735 0,418 | 0,235 0,338 
os ad 0,760 | 0,62 | 0,338 0,558 


et or 5 
pa <4 


( ! us Zusammenfassung. Der Zusatzauftrieb des VerdichtungsstoBes 14Bt sich berechnen, indem 
man den Sto als Kompressionssenke auffa8t, mit welcher ein System raumlich verteilter Ex- 
-_pansionsquellen entgegengesetzt gleicher Gesamtergiebigkeit verbunden ist. Der auftriebvermin- 
: dernde EinfluB der Expansionsquellen kann in erster Naherung nach der Stromlinienanalogie 


ermittelt werden. . 
= core, ee oh aa ae ee oh by eee : 
aN Siehe Fufnote 1 von Seite 280. ; i { 
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pidwlans ir CL eu) | (Eingegangen am 11. August 1948.) ve 
a Anschrift des Verfassers: Professor Heinrich Bruno Helmbold, (13b) Miinchen 23, RémerstraBe 14 II. 
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Untersuchungen iiber die Wandschubspannung in turbulenten 
Reibungsschichten. 
Von H. Ludwieg und W. Tillmann. 
1. Einleitung. Die Wandschubspannungen in laminaren Grenzschichten lassen sich auf rein 


theoretischem Wege berechnen, da man hier den Zusammenhang von Geschwindigkeitsprofil 
und ubertragener Schubspannung kennt. Da dieser fiir die durch den turbulenten Austausch 


hervorgerufenen Schubspannungen noch fehlt, ist ein ahnliches Vorgehen bei den turbulenten | 


Reibungsschichten nicht méglich. Aus’diesem Grunde ist man bei der Ermittlung der Gesetze 
fir die turbulente Wandreibung auf experimentelle Untersuchungen angewiesen. Diese sind 
in zwei einfachen Fallen, an einer nicht angestellten, ebenen Platte — im folgenden kurz ,,Platten- 
stromung*’ genannt — ‘und in einem Rohr bzw. Kanal von flach rechteckigem, konstantem 
Querschnitt — im weiteren mit ,,Rohr- bzw. Kanalstrémung“ bezeichnet — zur Geniige durch- 
gefihrt worden. Aus den Versuchsergebnissen haben die einzelnen Bearbeiter einfache Nahe- 
rungsformeln fiir den Reibungswiderstand gefolgert, die einigermaBen tibereinstimmen. 

-Einige Untersuchungen sind auch an Reibungsschichten bei Druckanstieg oder -abfall an- 
gestellt worden, die aber teils keine, teils unbefriedigende Aussagen tiber die Wandreibung ent- 
halten. Diese Versuche wurden in einem Kanal von kreisférmigem! oder meistens rechteckigem 
Querschnitt ausgefihrt?, 3, 4, 5. Bei letzteren war eine Kanalwand als ebene MeBplatte her- | 
gestellt, an der die zu untersuchende Reibungsschicht vermessen wurde. Die gegeniiberliegende jf 
Wand konnte zur Einstellung des gewiinschten Druckverlaufes verstellt werden. Sie war von 
der Versuchsflache soweit entfernt, daB zwischen beiden Reibungsschichten noch ein Kern mit 
Potentialstromung erhalten blieb (freie Reibungsschicht). Aus den gemessenen Geschwindigkeits- 
profilen wurde mit Hilfe der von v. Kérmdén angegebenen Impulsgleichung die Wandschubspan- 
nung ermittelt. Der Vorteil der beschriebenen Versuchsanordnung liegt darin, da man mit ver- 
haltnismaBig kleinen Anlagen dicke Reibungsschichten erzeugen kann (groBe Reynoldssche 


Zahlen). Ein wesentlicher Nachteil ist die im Verhaltnis zur Reibungsschichtdicke geringe Breite jf 


der Versuchsstrémung. Hierdurch kénnen leicht Sekundarstromungen entstehen, welche die fir 
die Auswertung nach der v. Kérmdnschen Impulsgleichung vorausgesetzte Zweidimensionalitat 
der Strémung aufheben und vermutlich fir die unwahrscheinlichen Ergebnisse der oben er- 
wahnten Autoren verantwortlich sind. Einige von ihnen hatten namlich festgestellt, daB in stark 
verzogerten Strémungen der értliche Widerstandsheiwert c;(=Tw/Q, wo Tw die Wandschub- 
spannung und Q der Staudruck au®erhalb der Reibungsschicht ist) nach einer gewissen Lauf- 
strecke in Strémungsrichtung plétzlich auf ein Mehrfaches des Anfangswertes anstieg, und nicht, 
wie man eigentlich analog dem Verhalten der laminaren Grenzschicht erwarten sollte, abfiel. 
Da sich hierfir keine stichhaltige Begriindung finden lieB, ist von W. Tillmann® untersucht 
worden, ob dieser Effekt nicht durch Sekundarstrémungen vorgetauscht wurde. In einem Schnitt 


quer durch den Kanal wurde die Geschwindigkeit nach GréBe und Richtung bestimmt. Eine 


Abschatzung des Einflusses der hierbei festgestellten Sekundarstrémungen ergab einen um 40 % 
erniedrigten cz-Wert. Damit war das bisherige MeBverfahren in Frage gestellt. 

Neuerdings ist von H. Ludwieg’ ein einfaches MeBgerat entwickelt worden, mit dem man an 
beliebigen MeBpunkten mittels einer Warmeiibergangsmessung die Wandschubspannung be- 
stimmen-kann. Dieses direkte MeBverfahren ist von eventuell auftretenden Sekundarstrémungen 
unabhangig. Mit dem neuen Méfgerét wurden die Untersuchungen tber den Reibungswiderstand 
turbulenter Reibungsschichten bei Druckanstieg und -abfall wieder aufgenommen; wortber in 
der vorliegenden Arbeit berichtet werden soll. . 


1 A, Buri, Eine Berechnungsgrundlage fiir die turbulente Grenzschicht bei beschleunigter und verzégerter 
Grundstrémung, Diss. Ziirich 1931. 


* E. Gruschwitz, Ing.-Arch. 2 (1931) S. 321. 

° A. Kehl, Ing.-Arch. 13 (1943) S. 293. ; 
-. * W. Mangler, Das Verhalten der Wandschubspannung in turbulenten Reibungsschichten mit Druck- 
anstieg. ZWB, UM 3052-(1943). é : 

_ ® K. Wieghardt, Uber die Wandschubspannung in turbulenten Reib hichten bei veranderli - 

druck, ZW, UM 6603 (1913) Pp g tare ibungsschic Le aa ae AuBen: | 

: W. Tillmann, Uber die Wandschubspannung turbulenter Reibungsschichten bei Druckanstieg, Diplom- 
arbeit. Géttingen 1947. Gu a ; 

” H. Ludwieg, Ing.-Arch. 17 (1949) S. 207. . 
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i 2.. Theoretische Betrachtungen iiber den Reibungswiderstand. Die bisherigen Untersuchungen 
aber Platten- und Rohrstrémung fihrten zu folgenden Ergebnissen. 


a) Das Geschwindigkeitsprofil der Reibungsschicht ]a8t sich stets in der nach- 
stehenden Form darstellen!: 


Wi taee aa 
: a t8ls-) : (1) 
(w= Geschwindigkeit in der Reibungsschicht; U=Geschwindigkeit auBerhalb der Rei- 
4 . . i (ee) 
bungsschicht ; y = Wandabstand; 6,= f Z (1 a) dy = Impulsverlustdicke der Rei- 
I 0 EF je 1 


U 
bungsschicht ; re= U6,/y = Reynoldssche Zahl der Reibungsschicht, gebildet mit der Impuls- 
verlustdicke 6,; »=kinematische Zahigkeit). 
In (1) ist g eine feste Funktion, die aber natirlich fir die Platten-, Rohr- und Kanal- 
stromung je verschieden ist. Die Abhangigkeit u/U von re ist nur gering, d. h. die Ge- 
schwindigkeitsprofile ftir verschiedene re-Zahlen unterscheiden sich nur wenig. 


a Der értliche’ Reibungsbeiwert cy kann immer in der folgenden Form dargestellt 
werden: 


cy = F(re) (2) 
(cy = rol Oe r= Wandschubspannung; o= Dichte). 


' Fist wieder je eine feste Funktion fir Platten-, Rohr- und Kanalstromung. Bei der 
Plattenstromung lat sich F mittels des Impulssatzes berechnen, wenn g bekannt ist, da 
hier der gesamte Reibungswiderstand als Impulsverlust in der Reibungsschicht erscheint. 


_ te) Bir den wandnahen Teil der Geschwindigkeitsprofile ergab sich ferner die Be- 
ziehung eh 
; u u* y : 

: ae u* = f( v ) = ? (3) 
(ut =] tu/o= Schubspannungsgeschwindigkeit). 

Diese Beziehung gilt mit ein und derselben Funktion f sowohl fir den wandnahen 
Teil der Platten- wie fir den der Rohr- und Kanalstrémung. Fir nicht zu kleine 
u* y/y-Werte (voll turbulenter Bereich, u*y/y>50) kann (3) mit groBer Genauigkeit 
durch folgende Naherungsformel ersetzt werden: | 


u 
Tare a log 


ae te b, (3a) 


_ wobei a und 6 universelle Konstanten sind. Dieses logarithmische Gesetz laBt sich durch 
ein Potenzgesetz ees é 


ee ee 


annahern, wobei C-und n Konstanten sind, die noch ein wenig:von dem u* y/y-Bereich 
_abhangen, fiir den die Annaherung besonders gut sein soll. 


Wie in der erwahnten Arbeit von H. Ludwieg bereits ausgefiihrt wurde, ist zu erwarten, 
da® das universelle Gesetz (3) bzw. (3a) auBer fiir die Platten-, Rohr- und Kanalstroémung auch 
bei weit allgemeineren Reibungsschichtstrémungen in Wandnahe noch giiltig ist. Selbst bei 
Geschwindigkeitsprofilen, die weit von den Profilen der Plattenstrémung abweichen und z. B. 
bei Strémungen mit starkem Druckanstieg oder -abfall auftreten, kann man annehmen, daB in 
Wandnahe noch immer das Gesetz (3) bzw. (3a) gilt. Nie sh 
_ Eine gewisse experimentelle Bestatigung fir die allgemeine Giltigkeit des Gesetzes (3) baw. 
(3a) ergibt sich aus der folgenden, von K. Wieghardt? zuerst festgestellten Tatsache: Wenn man 
die Reibungsschichtprofile in der Weise log u/U uber log y/6, auftragt (Abb. 1), so erhalt man fir 
kleine y/d, lauter parallele Geraden. Also ist wu in allen Fallen der gleichen Potenz von y pro- 
portional. Aus der Neigung der Geraden bestimmt sich diese Potenz zu 0,13=1/7,7, was fiir 


1 In dieser Formel haben wir nicht, wie sonst iiblich, die mit der Lauflange, sondern die mit der Impuls- 
verlustdicke 0, gebildete Reynoldssche Zahl als kennzeichnende GréBe gewahlt, da dies fiir Reibungsschichten 
mit veranderlichem Auf endruck zweckmafiger ist. | 
2K. Wieghardt, Z. angew. Math. Mech. 25/27 (1947) S. 146. 
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den in Frage kommenden u* y/v-Bereich gut mit (3b) ibereinstimmt. Dies ist jedoch kein zwingen- 

der Beweis far die Giltigkeit der Gleichung (3b), da wir durch die Profilmessungen nur die Po- |} 

tenz von y, aber nicht die Konstante C prifen kénnen, weil u* unbekannt ist. a! 
‘Bei Giltigkeit von (3) fiir den wandnahen Teil der Reibungsschicht hangt u* und damit Tw |fj 

bzw. c; nur von dem Geschwindigkeitsprofil und den Stoffkonstanten des strémenden Mediums |} 

ab. Es laBt sich daher c; berechnen, wenn das Geschwindigkeitsprofil bekannt ist. Emme ent- 

sprechende Beziehung zwischen cy, der Reynoldsschen Zahl re und-einem noch naher zu defi- 

" nierenden Profilparameter soll im folgenden abgeleitet werden. 


AS 


konstanter Druck 


mabiger Druckanstieg 


starker Druckanstieg 
Druckabfall 


Biax = 


7 
- 


. f=-log y 
5 ee | 1 re=286-10? 
x ‘a 


M2 = 12833 Pp = 97 


2 = re=798+107 
Hye = 13375. YP = G685 


3 re=448-70% 
Hy = 1442; Y = Q623 
4 re=265 -70% 


Mg=15h; p= 058, 


-5 -re=325 +10" 


Me~ Wari = 047, 


— 06 


Abb. 1. Einparametrige Schar von Geschwindigkeitsprofilen der Reibungsschichten mit Druckgradienten 
R in dimensionsloser, doppeltlogarithmischer Auftragung. a e 


is 


___ Als Profilparameter fahren wir y=us,/U ein. Dabei ist us, folgendermaBen definiert: Wenn _ 
_ das Gesetz (3), (3a), (3b) bis zu y-Werten gilt, die gréBer als 6, sind, so ist 25, einfach der Wert ; 
von u an der Stelle y=6,. Gilt dagegen das Gesetz (3), (3a), (3b) nur bis zu 4 Werten, die kleiner | 
als 0, sind, so ist us, der Wert, den u annehmen wirde, wenn das Gesetz (3), (3a), (3b) bis zum 
 Punkt y=6, gelten wirde*. Man erhilt also aus der doppeltlogarithmischen Auftragung von 
u/U uber y/d, (siche Abb. 1) den Profilparameter y, indem man den geradlinigen Profilteil bis 
 y/6g=1 verlangert und den zugehérigen Wert von us,/U abliest. ai EN Ral 
- Zur Ableitung der oben erwahnten Beziehung zwischen cy, re und y gehen wir von (3) aus 
Zur Einfihrung des Profilparameters setzen wir hierin y=0, und u=up,. Es ergibt sich dain 


poserrersh 


so 


Ja 


Fans bh atime 


sy es 


ane 


Uz) 


wees, ° $i 92 : 
1 E, Gruschwitz, a. a. O., definierte als Profilparameter die GréBe n=1—(us,/U)*. Dort ist aber immer 


7. 


unter 45, der Wert von u an der Stelle y=6, zu verstehen, 
ay } . ‘ a, ie . i A Te i * i AD ne ‘ Oy 


/ i 


sonata eS NILE ane nTiRaE 


ae Gleichung besagt, daB ein direkter Zusammenhang zwischen Ru und eso besteht. Es 
u v 


59 
Shik (een): 


U5, v 
wobei die Funktion h durch die Funktion f bestimmt ist. Durch Einfihren von 


ae also gelten 


4 


u* = “ = z, | =y und — = Te 
erhalt man nach einfacher Umformung cea: 3 
OEE CF =2y*h?(re ae y? H(re y). (5) 


fice ist zur ‘Abkiirzung te Funktion 2 h?=H gesetzt. 

_ Far alle turbulenten Reibungsschichtprofile, deren wandnaher Teil durch das allgemeine 
Gea (3) dargestellt wird, ist also der Reibungsbeiwert in der Form ( (5) gegeben. Um in dieser 
Gleichung nun die Funktion h bzw. H zu bestimmen, kénnten wir (4) fir-den in Frage kommen- ie 
den Argummentbercich gemaB (3a) durch . ae 


ahs Ud, Us, 0 ne B a 
Pees cee a alog * St +b = log ( at 
rsetzen und diese Gleichung elas nach u*/us, =h auflésen. Da jedoch die Konstanten 


und bd in (3a) nicht besonders genau bekannt sind, scheint uns das folgende~ Verfahren’ zur G ss 


Nach (1) ist der Profilparameter y fur die Profile der Plattenstrémung, den wir mit y, be- 
zeichnen wollen, nur von re abhangig, also Yo= Yo (re). Setzen wir diesen Wert in (5) ein, so muB— 73 
ms diese Cleichung | den Widerstandsbeiwert c; fir die Plattenstrémung liefern. Unter Berick- of : 
si Pocus. von act bee eyes also fiir FH folgende Bee cehetet 


22H (ro) = Fle), 


2), ist die Konvergens ms K ame so gut, 
nun¢ | im : ten Fakto das | ae ae 


\ 
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Sie wurde unter der Voraussetzung abgeleitet, daB in Wandnahe das universelle Gesetz (3) gilt. | 
Die in ihr auftretenden Funktionen F und y, kénnen den experimentellen Untersuchungen der |} 
Plattenstrémung entnommen werden. ; 

Fir eine tiberschlagige Rechnung geniigt es, in (8) diejenigen Funktionen F und yp einzusetzen, 
die aus der Annahme des bekannten 1/7-Potenzgesetzes fiir das Geschwindigkeitsprofil folgen. | 
Wegen der Affinitat der Profile ist yg hierbei von re unabhangig und kann durch eine einfache 
Integration berechnet werden, die den Wert y)=0,717 liefert. Die entsprechende Funktion F 
mit der von E. Gruschwitz angegebenen Zahlenkonstanten lautet 


c; = F (re) = 0,0251 re-V/* . 
Setzt man diese Werte fir y, und F in (8) ein, so folgt 
Cf = 0,0449 y7/4 re“"4 . 


Da das 1/7-Potenzgesetz fiir die Geschwindigkeitsverteilung und das daraus folgende 1/4-Potenz- 
gesetz fir den Reibungsbeiwert nur in roher Naherung gilt, ist die eben abgeleitete c, -Formel |} 
verhaltnismaBig ungenau. Eine bessere Anpassung an die tatsachlich auftretenden Widerstands- | 
beiwerte erreicht man durch eine leichte Anderung der Zahlenkonstanten, wodurch sich folgende 


Formel ergibt: 
c = 0,0580 yi" re-%788. (9) 


Mit dieser Formel wird die Gleichung (8) recht gut angenadhert, wenn man in dieser fir die Funk- 
tion F das Schultz-Grunowsche Plattenreibungsgesetz (siehe folgenden Abschnitt) und fur die 
Funktion y, den in Abb. 2 dargestellten Verlauf einsetzt. In dem Bereich 1-108 S re [4-104 
bleiben die Abweichungen kleiner als 3 %. 


Abb, 2. Profilparameter yg der Plattenstrémung als Funktion der Reynoldsschen Zahl re, 


{ 


Aus der einfachen Naherungsformel (9) sieht man, daS bei konstantem re der Widerstands- 
beiwert cy proportional y+" ist. Da bei Reibungsschichten mit Druckanstieg langs der MeB- 
strecke y abfallt und re ansteigt, nimmt also c; stark ab, was ganz im Gegensatz zu den Er- 
gebnissen von W. Mangler und K. Wieghardt steht, die hier einen starken Anstieg des cy-Wertes 
festgestellt atten. Wir sahen uns daher veranlaSt, zur Bestatigung der in diesem Abschnitt — 
fir den Reibungsbeiwert cy abgeleiteten Beziehung (8) Versuche anzustellen, die im folgenden 
geschildert werden sollen. ; ‘ 

-3. Verhalten des Widerstandsbeiwertes in Reibungsschichten mit Druckgradienten. Die Rei- 
bungsschichten wurden in der gleichen MeBstrecke von rechteckigem Querschnitt untersucht, 
die bereits von F. Schultz-Grunow! beschrieben wurde und nach dem fir Reibungsschicht- 
messungen iblichen Prinzip, das schon in der.Einleitung erlautert wurde, gebaut war. Wahrend © 
bei F’. Schultz-Grunow die verstellbare Wand so eingestellt wurde, daB langs der gesamten Mef- 
strecke konstanter Druck herrschte, wurde sie hier so eingestellt, daB der gewunschte Druck- 
verlauf zustande kam. An 10—12 Stellen entlang der Mittellinie der ebenen glatten MeBplatte 
(1,4 mal 6 m) wurden die Geschwindigkeitsprofile der Reibungsschichten ermittelt. An den glei- 
chen Stellen wurde die Wandschubspannung mit Hilfe des bereits erwahnten, neuen, auf ther- 
mischer Grundlage beruhenden Mefgerates von H. Ludwieg gemessen. Insgesamt wurden vier 
verschiedene MeBreihen . 


1 F, Schultz-Grunow, Lufo 17 (1940) S. 239. 


\ 


1 
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a) bei konstantem Druck in Strémungsrichtung (Plattenstrémung), 
b) bei maBigem Druckanstieg, 
c) bei starkem Druckanstieg, 


d) bei Druckabfall 
durchgefiihrt. 

Bevor die Messungen ausgefiihrt wurden, ist das Schubspannungsmefgerat in der folgenden 
Weise geeicht worden. Der Kanal wurde auf Plattenstrémung eingestellt und das MeBgerat an 
zwei verschiedenen Stellen der MeBplatte eimgebaut. Zur Aufnahme der Eichkurve wurde jedes- 
mal der gesamte Geschwindigkeitsbereich der Versuchsanlage durchfahren. Hierbei wurde je- 
weils die zugehérige Eichschubspannung mittels des Reibungsgesetzes fiir die Plattenstrémung 
ermittelt. Aus den vorliegenden Naherungsformeln fir den Widerstandsbeiwert der Platten- 
stromung haben wir die von F. Schultz-Grunow aufgestellte} 


ue maicey 
\ “f = “Cog re) 1,838 


verwendet, da diese an Messungen in derselben Versuchseinrichtung gewonnen wurde. An- 
schlieBend haben wir dann zur Kontrolle bei konstant gehaltener Geschwindigkeit die Wand- 
schubspannung mit dem geeichten MefSgerat entlang der gesamten Versuchsstrecke gemessen. 
Aus den gleichzeitig gemessenen Geschwindigkeitsprofilen wurde die Funktion y)=y,(re), 
die wir bei der Prifung von (8) brauchen, ermittelt. Sie ist in Abb. 2 ber log re aufgetragen. 
In der gleichen Abbildung ist das aus den Schuliz-Grunowschen Messungen ermittelte yy ein- 
gezeichnet. Unsere MeBpunkte liegen bei kleinen re-Zahlen etwas oberhalb der Schultz-Grunow- 
schen Kurve. Den spateren Auswertungen haben wir die stark ausgezogene Kurve zugrunde- 
gelegt. ; : 
In Abb. 3a ist. fiir die vier verschiedenen MeBreihen der Widerstandsbeiwert c; tber der 
re-Zahl doppeltlogarithmisch aufgetragen. Zum Vergleich ist das Reibungsgesetz von Schultz- 
Grunow, mit eingezeichnet. Die MeBpunkte der zur Kontrolle durchgefihrten MeBreihe bei 
konstantem Druck (Plattenstrémung) fallen mit der Schultz-Grunowschen Kurve zusammen 
und bestatigen damit die Richtigkeit der Eichmessung. Die beiden MeBreihen bei maBigem bzw. 
starkem Druckanstieg ergeben am Ende der Mefstrecke cy-Werte, die erheblich unterhalb der 
Schultz-Grunowscheh Kurve fir die Plattenstrémung liegen. Dies bedeutet also, daB die An- 
-nahme eines nur von re abhangigen c;-Wertes, die 4. Buri, E. Gruschwitz, H. B. Squire und 

A. D. Young®, A. Kehl, A. E. v. Doenhoff und N. Tetervin® ihren Berechnungsverfahren fir 

turbulente Reibungsschichten mit Druckgradienten zugrundegelegt haben, nicht zutrifft. Erst 

recht wird das Ergebnis der Untersuchungen von W. Mangler und K. Wieghardt widerlegt, 
die bei Reibungsschichten in verzogerter Strémung nach einer gewissen Anlaufstrecke in Stré- 
“mungsrichtung ein erhebliches Ansteigen des c;-Wertes festgestellt hatten. In der Abbildung 

sind noch die MeBpunkte der Versuchsreihe bei Druckabfall eingetragen. Diese Punkte liegen 
-ein wenig oberhalb der Schultz-Grunowschen Kurve. 
Nach diesen Vorbereitungen gehen wir zur Priifung von (8) tiber. Zur Vereinfachung der ® 

-Schreibweise dieser Gleichung fihren wir noch folgende Abkirzung ein: 


e id ae 89 
ae odie) =m. 
Wenn wir (8) dann in die nachstehende Form bringen 
i : ia Ay é ss 
; ‘ ie YO a fie y 8 
ool) =F (re sotrey)” , > 


” besteht zwischen den -beiden Ausdricken cr (Yoly)? und re y/yo(re) der gleiche funktionale 
Zusammenhang wie zwischen den GréBen cz und re bei der Plattenstrémung. Wir haben daher 
‘im Abb.3b log les (2) ) iiber | log (re a) aufgetragen und das Plattengesetz von Schultz-Grunow 
=e: L oY ¢ EV QF: 
aes ra der Arbeit von F. Schultz-Crunow ist c’f als Funktion der mit der Lauflinge x gebildeten Reynoldsschen 
-Zahl angegeben. Wir haben diese Formel auf die mit der Impulsverlustdicke 6, gebildete Reynoldssche Zahl re 
‘umgerechnet. pe ce aes ! ; 
> 2 H. B. Squire and A. D. Young, The calculation: of the profile drag of aerofoils. Tech. Rep. of the Aeron. 
‘Res. Comm., R. and M. No. 1838 (1938).' _ ue ieee ay Sas 

8 A,B. 0. Doenhoff and N. Tetervin, Determination of general relations for the behaviour of turbulent boun- 
dary layers. NACA Rep. No. 772 (1943). — eS re keel 


ae ay ~ zm 
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mit eingezeichnet. Bei Richtigkeit unserer Gleichung (8) miissen jetzt die Punkte simtlicher 
MeBreihen auf die eingezeichnete Kurve fallen. Die Punkte der MeBreihe mit konstantem Druck | 
wurden einfach aus Abb. 3a tibertragen, da sich eine Umrechnung eribrigt. Die Punkte der |} 
restlichen MeBreihen liegen befriedigend auf der Kurve von Schultz-Grunow. Insbesondere fallen | 
die MeBpunkte der Versuchsreihen mit Druckanstieg, die zu erheblich kleineren cp-Werten als 
den entsprechenden Werten der Plattenstromung gefihrt hatten, gut mit der Schultz-CGrunow- | 
schen Kurve zusammen, womit unsere Gleichung (8) auf das beste bestatigt wird. Lediglich 
die umgerechneten Punkte der Reihe mit Druckabfall liegen ein wenig unterhalb der eingezeich- | 
neten Kurve. Dies braucht kein Versagen unserer Gleichung (8) zu bedeuten, da das Gesetz | 
von Schultz-Grunow bei kleinen re-Zahlen ziemlich unsicher ist, was aus dem Vergleich mit | 
Naherungsformeln anderer Verfasser}, 2, °, 4, > hervorgeht, die hier von der Schultz-Grunowschen | 
Kurve um die gleiche Gré®enordnung abweichen wie unsere Mefpunkte. 
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Abb, 3. Widerstandsbeiwerte ef in den Reibungsschichten mit ‘Druckgradienten 
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Aus unserer Gleichung (8) (oder einfacher nach. der Naherungsgleichung (9) einzusehen) 
folgt, da mit Annaherung an die Ablésestelle einer turbulenten Reibungsschicht (y—>0) der 
Widerstandsbeiwert cr gegen Null geht. Kurz vor der Ablésestelle miissen daher sehr kleine i 
cy-Werte auftreten, die wir mit der MeBreihe bei starkem Druckanstieg nachzuweisen versucht ee 
haben. Dabei erreichten wir einen cy -Wert von 0,0010 anstatt, eines CF von 0,0020 bei der Platten- ee 
stromung mit gleicher re-Zahl. Noch kleinere Widerstandsbeiwerte lieBen sich in der verwendeten . 
Versuchseinrichtung nicht einwandfrei messen, da die Versuchsstrémung in den Ecken des 
Kanalquerschnittes vorher abri&. 

Bei der Ableitung von (8) wurde vorausgesetzt, daf} das universelle Gesetz (3) far den wand- 
nahen Teil der Geschwindigkeitsprofile in der Reibungsschicht auch hei allgemeinen Reibungs- 
‘schichten giiltig ist. Es miissen also alle Profile in der Auftragung u/u* tiber log u* y/y in Wand- 
nahe zusammenfallen. Dieses Verhalten wird an einigen, beliebig aus den vier MeGreihen aus- 
gewahlten Profilen in Abb. 4 gut bestatigt. Die hierzu notwendigen Schubspannungsgeschwindig- 


© re= 798-10 


mabiger Druckanstieg 
+ re=265-70* 


starker Druckanstieg 
x  re=7148-70# 
o 6re=3,25-10* 


Drackabfall 
a re=386-103 


20 


: ; 20 35 30 35, 40 
ae Abb. 4, Einige Geach windiskeitsprotle von Reibungsschichten mit Druckgradienten mit universellem Gesetz (3a) in Wandnihe. 
aE aes ; ‘ : S ‘ : i ; 


m voll turbulenten Bereich’ des -Geschwindigkeitsprofiles, ‘sondern schon im Ubergangsgebiet 


witz gefundene 


in Wandnahe zusammen. Hieraus geht nun die Auftragung 
Fruschwitz simtliche Profile eine einparametrige Schar bilden sollen, 
Math Mech. 20 (1940) S0297.- Sa, . 
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durch affine Verzerrung der Ordinate mit u*/U und der Abszisse mit y/u* 6, hervor. Da die beiden 
GréBen u*/U- und »/u* 6, voneinander unabhangig sind, sollte man erwarten, das schon der 
wandnahe Teil des Geschwindigkeitsprofiles in der Auftragung u/U uber y/d. zweiparametrig 
wird. Nun kann man aber das allgemeine Geschwindigkeitsgesetz (3), das fur den in: Abb. 4 
zusammenfallenden wandnahen Profilteil gilt, durch ein einfaches Potenzgesetz (3b) auSerhalb 
der laminaren Unterschicht annahern. Ein Potenzgesetz geht bei beliebigen, voneinander un- 
abhangigen affinen Verzerrungen in Ordinaten- und Abszissenrichtung stets wieder in Potenz- 
gesetze mit gleichem Exponenten tber, die eine einparametrige Profilschar in Wandnahe ergeben.  }}j 

Der anschlieBende Teil des, Geschwindigkeitsprofiles, in dem sich die Geschwindigkeit nur 
noch langsam dndert, ist bei gegebenem wandnahen Teil weitgehend festgelegt, da der AnschluB_ 
stetig und mit stetigen Ableitungen erfolgen, das Profil sich asymptotisch dem Wert u/U=1_ 
annahern und das Integral 


(Definition der Impulsverlustdicke) sein muB. . 
Hiermit ist die angenadherte Einparametrigkeit fiir das ganze Geschwindigkeitsprofil ver- 
standlich gemacht. Diese Einparametrigkeit bezieht sich jedoch nur auf den turbulenten Teil 
des Profiles und schlie8t die laminare Unterschicht nicht mit ein. 
Da wegen der Einparametrigkeit sich die Geschwindigkeitsprofile durch einen beliebigen 
Parameter eindeutig kennzeichnen lassen, mu zwischen zwei Profilparametern ein eindeutiger 
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+ mabiger Druckanstieg 


. x storker Druckanstieg 
a A Druckabfoll 
2 * konstanter Druck 
‘ und Turbulenzgitter ins 
ee a konstanter Druck 
und Leisteé 


Se 


i 4 - 
Sok Sn, =2,333-10~ 99982 ees 
hg 
=] 
— 
Miz “ 
5 16 17 


Abb. 5. Zusammenhang zwischen den beiden Profilparametern y und eas 
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Zusammenhang bestehen. Dieser ist in Abb: 5 fiir y und die GréBe H,, dargestellt, wobei H. 

_ das Verhaltnis von der Verdrangungsdicke! 6, zur ImpulsverJustdicke 6, ist. In dieser Abbildeng 
¢ liegen die MeBpunkte, die den verschiedensten Profilen aus allen MeBreihen entsprechen, recht 
gut auf einer Kurve (mit Ausnahme. der Versuchsreihe ,.konstanter Druck und Leiste‘, siche 
folgenden Abschnitt), womit die Einparametrigkeit auf das beste bestatigt wird. Hy, ist deshalb 
ein. zweckmaBiger Profilparameter, weil diese GréBe in den Rechenverfahren fir turbulente 
Reibungsschichten haufig verwendet wird, und es liegt nahe, sie an Stelle von y in die Glei- 
chung (9) einzufihren. Der in Abb. 5 dargestellte Zusammenhang zwischen y und H,, laBt sich 
durch folgende empirische Formel-annahern: © aie cm 


yy = 2,333 - 10-9398 Hie, fe 


{ 
7 


‘ 


1 V d - Cay aes 22 u A : 
: er rangungsdicke Coeur ( 2b =) dy. : of er fee Ms 


Durch Einfihren dieser Beziehung in (9) erhalten wir die Naherungsformel 
Cp = 0;246.- 1079078 Biz pe 0208 | - (9a) 


die die Widerskandalichwsete cy in Abhangigkeit vom Profilparameter H,, und der Reynoldsschen 
Zahl re gut wiedergibt. 


5 Anwendungshereich der gefundenen GesetzmaBigkeiten. Nachdem die Ergebnisse, die in 
den vorhergehenden Abschnitten gefunden wurden, sich bei Reibungsschichten mit Druck- 
anstieg und -abfall durchaus bewahrt haben, wollen wir jetzt prifen, inwieweit die abgeleiteten 
Beziehungen bei irgendwelchen Stérungen der Reibungsschichtstromung noch giltig bleiben. 
Dazu haben wir noch zwei zusatzliche MeGreihen 4 


e) bei konstantem Druck und Turbulenzgitter, 
f) bei konstantem Druck und Leiste 


durchgefihrt. Fur die Messung der ersten Versuchsreihe wurde ein bereits von K. Wieghardt1 be Se 

_hbeschriebener Versuchsaufbau verwendet. Dabei war vor der Kanaldise ein aus Blechstreifen 

is angefertigtes, grobes Gitter angebracht, das die Turbulenz des Luftstromes betrachtlich erhéhte 
_(Erniedrigung der kritischen Kugelkennzahl U D/y von 3,75 - 105 auf 1,3- 10°). Der weitere Ver- ‘i 
-suchsablauf war der. Bote wie bei der aa mit konstantem Hs (a). In Abb. 6a ist wieder 
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misch aufgetragen. Die MeBpunkte dieser Versuchsreihe fallen dann mit den Mef{punkten der 
Plattenstrémung zusammen, wodurch auch die Anwendbarkeit unserer Gleichung (8) auf der- 
artig stark turbulente Strémungen erwiesen ist. 

Im Gegensatz zur soeben besprochenen MeBreihe e), wo die Reibungsschicht vom AuBenrand 
her gestért wurde, wurde sie bei der MeBreihe f) von der Wand her gestért. Dazu wurde auf 
der MeBplatte quer zur Strémungsrichtung eine durchgehende Vierkantleiste (13 mal 13 mm) 
in einer Entfernung x=3m von der Plattenvorderkante angebracht. Die Dicke der Reibungs- 
schicht betrug an dieser Stelle ca. 60 mm. In verschiedenen Abstanden hinter der Leiste wurden 
wieder Geschwindigkeitsprofile und Wandschubspannungen gemessen. In Abb. 6a sind die MeB- 
werte fiir c; eingezeichnet. In kiirzeren Abstanden hinter der Leiste liegen die MeBpunkte natur- 
gemaB tiefer als bei der Plattenstromung. Mit zunehmender Entfernung steigen dann die c;- 
Werte steil auf gréBere Werte als bei der Plattenstroémung an und fallen bis zum Ende der MeB- 
strecke noch nicht bis zu den Werten der Plattenstrémung ab. Die gegen die Plattenstrémung 
erhéhten c;-Werte sind durch die verstarkte Turbulenz bedingt, die durch die Ablésung an der 
Leiste verursacht wird. In Abb. 6b sind die entsprechend (8a) umgerechneten c;-Werte ein- 
getragen. Der erste Punkt dicht hinter der Leiste fallt nicht auf die Kurve des Plattengesetzes, 
wahrend schon nach gréBenordnungsmafig zehn Reibungsschichtdicken die MefSpunkte die 
Gleichung (8) wieder einwandfrei erfillen. 
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konstanter Druck 
und Turbulenzgitter 


a_ konstanter Oruck 


mit Leiste 


re=121+ 10% 
Hye = 124g ; P= 964 

2: re=800-103 
Hye =126; 3 p=O725 

3: re=798- 10? 
Hye =1357 5 P =9685 

4: re=102-70" 
Hye= 16045 p= 956s 


-05 


K 
Abb. 7, Einige Geschwindigkeitsprofile der MeGreihen e und f in dimensionsloser, doppeltlogarithmischer Auftragung. 


Der Grund far das Versagen der Gleichung (8) fiir die ersten MeBpunkte liegt darin, daB die 
stets vorausgesetzte Giltigkeit der Gleichung (3) in Wandnahe durch die Leiste gestért ist. 


_ Zur Bestatigung dieser Tatsache sind in Abb. 7 entsprechend der Abb. 1 einige Profile dar- — | 


gestellt. Profil Nr. 4 ist ein kurz hinter der Leiste gemessenes Profil. Ein Vergleich mit Abb. 1 
zeigt, daB hier in Wandnahe das Gesetz (3) nicht mehr gilt. AuBerdem sieht man, daB das Profil 
aus der einparametrigen Schar herausfallt. Profil_Nr. 1 ist in dem, Bereich gemessen, wo (8) 
Wieder zutrifft. Wir erkennen den geradlinigen Verlauf mit derselben Neigung wie in Abb. 1 - 
entsprechend der Giltigkeit von (3). Zum Vergleich ist das Plattenprofil aus Abb. 1 ubertragen. 
. Ferner ist in Abb. 7 ein Profil der MeBreihe mit Turbulenzgitter eingezeichnet, welches wiederum 
die Giltigkeit von (3) zeigt, wie es nach obigen Ausfithrungen nicht anders zu erwarten war. 

In Abb. 8 sind aus den letzten MeBreihen einige Profile, far die (8) wieder gilt, entsprechend 
der Abb. 4 aufgetragen und zum Vergleich ein Profil der Plattenstrémung mit eingezeichnet. Wir 
erkennen auch hier wieder, daB alle eingezeichneten Profile in Wandnahe Gleichung (3) erfiillen. 
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6. Zusammenfassung. VeranlaGbt durch die unbefriedigenden Kenntnisse des Verhaltens der 
Wandschubspannung in Reibungsschichten mit Druckgradienten wurde diese Aufgabe erneut 
behandelt. Es wurde gezeigt, daB bei allgemeinen turbulenten Reibungsschichten in Wandnahe, 
d.h. in der laminaren Unterschicht, im Ubergangsgebiet und im wandnahen Teil des voll tur- 
bulenten Bereiches, das gleiche universelle Geschwindigkeitsgesetz wie bei der Plattenstrémung 
gilt. Aus der allgemeinen Giltigkeit dieses Gesetzes wurde eine Formel fir den ortlichen Wider- 
standsbeiwert c-; abgeleitet, in der c; nur von der mit der Impulsverlustdicke gebildeten Reynolds-. 
schen Zahl re und einem Profilparameter y abhangt. Diese GesetzmaBigkeit wurde durch mehrere: 
Versuchsreihen, bei denen die Wandschubspannung mit Hilfe eines neuen MeBgerates unmittel-. 
bar gemessen wurde, gut bestatigt. Man kann also allein aus der Kenntnis des Geschwindigkeits- 
profiles den zugehérigen Reibungsbeiwert cy bestimmen. 


se 
u* 
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konstanter Druck 
© re=798-107 


konstanter Druck 
und Turbulenzgitter 

* re=800-10% 
konstanter Druck und Leiste 
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Abb. 8. Einige Geschwindigkeitsprofile der MeGreihen e und f mit universellem Gesetz (3a) in Wandnahe. 


Aus der Formel fir c; folgt in Ubereinstimmung mit den Versuchen, daB bei Reibungsschichten , 
mit Druckanstieg bzw. -abfal] kleinere bzw. gréBere c;-Werte als bei der Plattenstrémung gleicher 
_re-Zahl auftreten. Es lassen sich also bei gréBeren Reynoldsschen Zahlen kleine drtliche Wider- 
standsbeiwerte cy nicht nur durch Laminarerhaltung der Grenzschicht, sondern auch durch 
geeigneten Druckverlauf in turbulenten Reibun sschichten erreichen. Das von verschiedenen 
Bearbeitern festgestellte Ansteigen des Reibungsbeiwertes auf ein Mehrfaches des Beiwertes 
der Plattens.rémung in Reibungsschichten mit Drickanstieg ist durch diese Untersuchung 
widerlegt. 


Max Planck-Institut fiir Strémungsforschung in Géttingen. 


(Eingegangen am 11. August 1948.) 
Anschrift der Verfasser: Dr. Hubert Ludwieg, Dipl. Phys. Walter Tillmann, 
(20b) Géttingen, BéttingerstraBe 6/8. 
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Uber Werkstoffdimpfung bei Biegeschwingungen *. 
Von G. Denkhaus. 


1. Einleitung. Wenn man den Resonanzausschlag eines schwingend beanspruchten Bauteils 
vorausberechnen will, so scheitert diese Aufgabe haufig daran, daB man die Dampfung quanti- | 
tativ nicht hinreichend kennt. Nun entfallt ja bekanntlich ein beachtlicher Anteil der gesamten 
Dampfung auf die Werkstoffdampfung, und es ware somit schon ein Fortschritt, wenn es gelange, 
den Resonanzausschlag wenigstens unter Beriicksichtigung der Werkstoffdampfung zu errechnen. 

Unter der Werkstoffdampfung versteht-man ja die Erscheinung, da ein Teil der bei 
Schwingungsbeanspruchung in einen Bauteil hineingesteckten Formanderungsenergie durch 
nicht reversible Vorgange im Werkstoff in Warme umgesetzt wird. 


2. Kenngréfen. Will man die Werkstoffdampfung zahlenmaBig erfassen und damit der 
Rechnung zuganglich machen, so .muf man sich zunachst uber die festzusetzenden Kenn- | 
groBen klar werden. Naheliegend ist es, die absolute Dampfungsarbeit S in (cm kg)'zu nehmen, | 
die wahrend einer Schwingungsperiode in Warme umgesetzt wird. Ein Mafi fir diese Arbeit | 
ist bekanntlich der Flacheninhalt der Hysteresisschleife, die von der Kraft-Verformungsweg- _ 


bei Zug- Druck: P kg) & 
beiTorsion: —M4lkg:c 
be/bregung. P (kg) 


Schwingungsumkehriage 


a) 


ie 
is 


oN) 


LongendnderungAl {cm} 
Torstonswinkel @ {a} 
Liegepteil f (cm) 


Tm Dimptungsarbert S 
arminderungsarbeit A in der Schwingungsumkehriage 


bei Normalspannung a 
bel Schubspannung © 


Be: 


Verh i, 

ie hae rie tim L/ementdéimptongstahigheit — 

Dehnunge as 

Schiebung y =F : { 
es SSS) 502. Hementforminderungsarbeit 


a ve c= % 


As ag Beet ~aV Lea is aa | 
__ Abb. 1. Hysteresisschleife bei Schwingungsbeanspruchung _ a) fiir einen Bauteil b) fiir ein Volumenelement d V. 


ea 


_Kurve im Verlauf einer Schwingungsperiode beschrieben wird (Abb. 1a). Experimentell ha 
sich ergeben, daB diese Dampfungsarbeit von folgenden Faktoren abhangig ist: ily 2 
1. dem Werkstoff, . 4! Mave, Pes : ae 2 
2. der Beakepriohuedatt (Zug-Druck oder Torsion oder Biegung), 
3. der Amplitude der Beanspruchung, ca ad etc cae 
4. der Form und GréBe des Schwingers, | . Fhe | 
5. sonstigen Faktoren, wie Temperatur, Vorgeschichte des Werkstoffs usw., die fiir die vor-_ | 
liegende Betrachtung nicht wesentlich sind. Bear i | see ‘Aj || 


1 lng 


_* Die Arbeit entstand aus einer Anregung von Herrn Professor. Dr. K. Kloter und wurde im Mai 1948 "| 
auf der Tagung der Fachausschiisse fiir Schwingungs- und Schalltechnik des VDI. in Stuttgart vorgetragen. | 
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Wahlt man die Dampfungsarbeit S als KenngréBe, so stért vor allem die Tatsache, daB sie 
von Form und Gré8e des Bauteils abhangig ist. Man braucht also eine Unmenge langwieriger 
Versuche, um sie zu erfassen. Es gelingt aber auch, eine Kenngré8e zu definieren, die den 
stoérenden Einflu8 nicht aufweist, und zwar auf folgende Weise: Man erstreckt die Betrachtung 
nicht auf den ganzen Schwinger, sondern auf ein beliebig kleines Element desselben. In einem 
solchen Element wird die Dampfungsarbeit dS in (cm kg) vernichtet. Diese Elementdampfungs- 
arbeit ist im wesentlichen von den gleichen Faktoren abhangig wie die Gesamtarbeit S, nur ist 
zu Punkt 2 und 3 zu bemerken, da hierher die Amplitude der Beanspruchung des Elements zu 
setzen ist, gemessen durch Spannung oder Formanderung, und daB ferner unter Beanspruchungs- 
art nunmehr die beiden Typen Normal- und Schubspannung zu verstehen sind. Hinsichtlich 
Punkt 4 entfallt der EinfluB der Form, von der man bei einem beliebig kleinen Element ja 
nicht mehr sprechen kann. Es bleibt allerdings noch der EinfluB der GréBe, beschrieben durch 
das Volumen d V des Elements, bestehen. ; 

Den Einflu8 des Volumens kann man nun aber eliminieren, indem man mit dem Quotienten 
dS/dV arbeitet. Dieser Quotient # in (em kg/em®) ist somit eine nur von Werkstoff und Be- 
anspruchung abhangige GréBe und heif®t Elementdampfungsfahigkeit. Sie wird veranschaulicht 
durch den Inhalt der Hysteresisschleife, die von der Spannungs-Verformungs-Kurve fiir das 
betrachtete Element bei einer Schwingungsperiode durchfahren wird (Abb. 1b). e 

Man kann aus der Elementdampfungsfahigkeit @ eine dimensionslose Gré®e ableiten, indem 
man sie auf die aus der Festigkeitslehre bekannte spezifische Formanderungsarbeit a=dA/dV 
bezieht. Hierbei ist dA die Formanderungsarbeit, die in der Schwingungsumkehrlage in dem 
Element aufgespeichert ist (vgl. Abb. 1a). Die neu definierte GroBe heiBt Elementdampfung: 

. o dS 
§ eat | Wi aa (1) 
_ Analog dazu kann man auch zu der gesamten Dampfungsarbeit S eine dimensionslose GréBe 
definieren, indem man sie auf die im ganzen Bauteil in der Schwingungsumkehrlage auf- 
_ gespeicherte Formanderungsenergie A bezieht. Diese KenngréSe heiBt Nenndampfung und 
_ gibt an, welcher Bruchteil der Formanderungsenergie im Bauteil durch Dampfung vernichtet 


wird: : | te 
ag Ad - wh { 


ee en es oe (2) 


Man muB also grundsatzlich zwei Gréfen unterscheiden: die WerkstiickkenngréBe p, und 
die WerkstoffkenngréBe y. DaB im allgemeinen ein Unterschied zwischen den beiden GréBen 
 besteht, liegt daran, da die Spannungsverteilung im Bauteil meist nicht gleichmaBig ist und 
- somit auch nicht alle Elemente gleichmaBig zur Gesamtdampfungsarbeit beitragen. 

__ Ermittelt man die Werkstoffdampfung experimentell, so erhalt man stets eine Nenndampfung, 
denn auch ein Probestab ist ja ein Bauteil mit bestimmter Form und Gréfe. Die Ermittlung 
-geschieht entweder durch Dauerversuche mit Dampfungsmessung oder durch Ausschwing- 
versuche. Die letztgenannte Methode ist recht bequem, denn bei schwacher Dampfung besteht 
‘ein einfacher Zusammenhang zwischen dem bekannten logarithmischen Dekrement der Aus- 


der Formanderungsenergien bezeichnen kann. Es gilt?: Su aad | 
e Shs BT oN SOE i ceca ian ag. 2.8 4: “ e (3) 
3. Einfluf der Beanspruchung. Wie schon betont, ist die Werkstoffdampfung keine Konstante, 


sogenannte Dampfungskurven darstellt, pflegt man die Beanspruchungsamplitude durch eine 
rformungsgréBe auszudricken, also bei Normalspannung durch die Dehnungsamplitude ¢ 
nd bei: Schubspannung durch die Schiebungsamplitude y. . 


Es ist natirlich erwinscht, diese Kurven in Formeln zu fassen. Die experimentellen Befunde 
rh ben ergeben, da® man sie durch Potenzfunktionen von folgender Form darstellen kann: 


By RA Si BBP BET pe RG ae e re PP aa 
Pe ea aa aden EO DEM, Pacha RO 


‘sind o und B ZW. a* und p* reine Werkstoffkonstanten. Die Werkstoffdampfung ist 


Go. Weydekanpf, Die Dauerpriifung der Werkstoffe. Berlin 1929. 


has = +h. > Pak 


schlage p=n Q, — In Q, 41 und der N enndampfung Wn» die man als das logarithmische Dekrement 


sondern im allgemeinen von der Beanspruchung abhangig. Wenn man diese Abhangigkeit durch - 


(4a) 


einer Potenz des Ausschlags proportional, man darf ihre KenngréBen also nicht ohne 
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weiteres in die bekannten Schwingungsdifferentialgleichungen fir geschwindigkeitsproportionale || 
Dampfung einfihren. 
Bei Bauteilen, in denen die Beanspruchung sich von Element.zu Element andert, hat es sich | 
als zweckmaBig erwiesen, als ,,NenngroBe* die Beanspruchung der am starksten verformten 
Elemente anzugeben; das sind bei Torsionsbeanspruchung alle Randelemente, und bei Biegung 
die Randelemente an der Stelle der gro8ten Kriimmung, also des maximalen Biegemoments. 
Fur die Elementdampfung, die sich ja nach (4a) mit der Beanspruchung andert, empfiehlt sich 
ein entsprechendes Vorgehen. Auch sie hat ja ihr Maximum in denselben Randelementen wie 
die Beanspruchung, und man kennt daher als Sonderfall der Elementdampfung die Rand- | 
dampfung!: 
PR= OR (4b) 


4, Werkstiickdaten und Werkstoffdaten. Die eingangs skizzierte Aufgabe, den Resonanz- | 
ausschlag unter Bericksichtigung der Werkstoffdampfung zu ermitteln, lieBe sich nunmehr mit | 
einem Mindestaufwand von Versuchen in folgender Weise lésen: Man ermittelt. zundchst ex- | 
perimentell die Kurve y,;=Yp (éx) fiir einen (z. B. zylindrischen) Probestab und rechnet diese | 
_Werkstiickkurve in eine Kurve y=wr (er) =a ef um, die nur vom Werkstoff abhangt. 

Es miSte nun méglich sein, aus dieser Werkstoffkurve fiir jeden Schwinger von beliebig 
gegebener Form und GréfSe wieder die entsprechende Kurve wrs=Yn (Er) zu errechnen und 
damit zahlenmafig die Dampfung anzugeben, die zur Beanspruchung des Bauteils gehért. 

Daraus ergibt sich die Notwendigkeit, Formeln aufzustellen, die es gestatten, Nenndamp- | 
fungskurven in Elementdampfungskurven umzurechnen und umgekehrt. Diese Formeln werden — 
letzten Endes bestimmt durch die Art der Verteilung der Spannung und damit auch der Damp- 
fung auf die einzelnen Elemente des Schwin- 
gers (Abb. 2). 

Bei Zug- Druck-Schwingern ist die 
Beanspruchung und damit auch die Ele- 
mentdampfung gleichmaBig verteilt, und 
somit gilt hier: 


Torsion P=Zz=W=?- (5) 


PO» 
DY 
or Biegung 


Bei Torsions-Schwingern ist die 
Beanspruchung zwar tiber die Lange hin 
gleichmaBig, aber den Querschnitt aber 
linear-polarsymmetrisch verteilt.- Somit 
verteilt sich auch die Elementdimpfung 
polarsymmetrisch, allerdings nicht linear, 
sondern wegen (4a) mit der Potenz f. 


Abb. 2.‘ Spannungsverteilung bei verschiedenen Beanspruchungsarten. Elementdampfung und Nenndampfung 
I werden hier also nicht tibereinstimmen. 
Die Ableitungen der Umrechnungsformeln fir Torsions-Schwinger findet man im Zusammen- 
hang bei K. Klotter?. Danach gilt fir Bauteile mit Kreisringquerschnitt 
ck 4 1—ottB 
Vn WR 4 +B 5 ea of 


mit yYr=a*yR- (6) . 


Der Ausdruck g=r,/r, ist das Verhaltnis von Innenhalbmesseér r, zu AuBenhalbmesser r, des — 
Ringes. Man sieht tbrigens, da nicht die absoluten Abmessungen in die Formeln eingehen, 
sondern lediglich der Verhaltniswert 0. = rN 

Aus der angegebenen Formel ergibt sich als Sonderfall fir den Vollkreisquerschnitt mit 9¢=0— 


A enti 
Po TPR eet « 3 ee (Y) 


Fir den Vollkreisquerschnitt entfallt also’ jeglicher Einflu8 der Abmessungen, und auch die 
Nenndampfung wird zu einer reinen WerkstoffgréBe, bleibt allerdings verschieden von der 
Elementdampfung. Ze ; 


1 Q, Féppl, E. Becker u. G. v. Heydekampf, Die Dauerpriifung der Werkstoffe. Berlin 1929... 
2 K. Klotter, Einfiihrung in die technische Schwingungslehre, Bd. 1, S. 184-186. Berlin 1938. 


é 
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5. Umrechnungsformeln bei Biegeschwingern. Umrechnungsformeln in der eben besprochenen 
Art waren fiir Biegeschwinger bisher noch nicht bekannt. Bei Biegebeanspruchung andert sich 
die Spannung und damit auch die Elementdampfung ja nicht nur itiber den Querschnitt, sondern 
auch tiber die Lange hin, und auBerdem ist die Anderung iiber die Lange noch vom Belastungs- 
fall, d.h. der Art der Auflagerung, Einspannung und Belastung, abhangig. 

Unter Beachtung der soeben erlauterten Einflisse konnen nun folgende Uberlegungen zur 


- Ableitung der gesuchten Umrechnungsformeln fihren: Definitionsgema8 gilt nach (2) fiir die 


Nenndaémpfung der Quotient 


S 
deg Be 


Der Zahler, die Dampfungsarbeit S des Schwingers, ergibt sich als Summe der Dampfungs- 

arbeiten dS aller Elemente. Aus der Definition der Elementdampfung in (1) folgt 
dS=ypdA=ywadV 

und damit fir den ganzen Schwinger unter Beachtung von y=a ef und a=} E &? 


S=f yadVauy [ ear. 
; P 
(V) (V) 


| Mit dV=dF dx und ¢e= ee wird Henne 


Bade 
) 
Sax (ez) J merase vit! dF. as (8) 


Der Nenner, die Formanderungsenergie A, ermittelt sich nach der Lehre von der Balken- 


- biegung zu 


l 
Az zz | M%x) dx. (9) 


Damit erhalt man fiir den ganzen Quotienten den Ausdruck: 


1 


F , " je J M?+8(x) dx 
Bl 
3 Wn = & (+5) + y2tB dF —_—______ ; : (10) 
f: (F) S M(x) dx 
ae 0 a 


Der EinfluB des Bauteils und des Belastungsfalls steckt nun noch in den Integralen iber 
M(x). Man kann ihn eliminieren, indem man die Integrale durch eine VerformungsgréBe, zweck- 


_ maBig die groBte im Schwinger auftretende Dehnung er, ausdriickt. Sie tritt bekanntlich am 


_ Rande an der ‘Stelle der gré8ten Krimmung auf. - 


: einfach durchfihren. Es empfiehlt sich daher, nunmehr die Belastungsfalle einzeln zu behandeln. - 


aay 


In der bisherigen, allgemein fir Biegung geltenden Form laBt sich diese Operation aber nicht 


6. Der beiderseits frei aufliegende Biegeschwinger. Als Musterbeispiel sei die Ableitung an 
dem frei aufliegenden Schwinger (Abb. 3) durchgefihrt. 


Man braucht zunachst die Funktion M=M(x), die Momentlinie. Dabei ist der Ursprung 


der Koordinate x zweckmaBig so zu wahlen, da M(x) auch in der Potenz (2+ 8), wo B eine 
-gebrochene Zahl sein kann, noch einfach integrierbar ist. In dem vorliegenden Beispiel kommt 


noch hinzu, da® die Momentenlinie infolge der Einzellast als Unstetigkeitsstelle der Belastung 
durch verschiedene Ausdriicke wiedergegeben wird und daher in zwei Intervallen integriert 


werden muB. 


Wegen M(x,) = Sue x, fir OS x%,Sb und M(x) = ey fir 0 S x, Sc wird das Inte- 


L 
as: = 1 z 4 ay 3 4 | ‘ 

gral { M?(x)dx zu : ee a ae ad ‘ 

Seer ge ; ‘ ‘ 7 

6b c Pb ; be Ly 1 Xp 

| ; c c 

f M? (x1) dx, + if M?(x2) dx_= ( 1 ) eC ii Abb. 3. Beiderseits frei 

; 0 Te Oe te fis aufliegender Biegeschwinger mit Einzelmasse. 
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und entsprechend gilt fir das Integral in der Potenz (2+) 
b 


. Pbc\2+8 b+e 
ff m2 (x)dx,+ f M2*(x)day= (=$4) Soe. 
0 0 ra 


Damit ergibt sich durch Einsetzen in die allgemeine Gleichung (10) fir das vorliegende Beispiel 


1 \o(Pbe\e 3 1 f ioe | 
oe Ziel der Ableitung ist es nun, einen Zusammenhang zwischen dieser Nenndampfung », und der 
(05 Elementdampfung wy zu finden. Da sich w von Element zu Element andert, empfiehlt es sich, 

a mit einer ausgezeichnet2n Elementdampfung, namlich der Randdampfung ypr=a eb, zu arbeiten. 
Fir den vorliegenden Belastungsfall ermittelt sich bekanntlich die gréBte Dehnung zu 
Mynax Pbe ; 
FRSC OR Dao er ded bey eee es 
Hierbei ist Mmax das groBte Biegemoment und y, der Abstand der auf ersten Randfaser des 
Querschnitts von der neutralen Faser. Setzt man den Wert Pbc/l aus (12) in (11) ein, so er- 
halt man 


Saks 
= B 2+6 dF. 


Mit «%=~pr ist die Randdampfung in dieser Formel enthalten. Der Einflu8 von Form und 
GréBe des Balkenquerschnitts steckt nun nur noch in dem Ausdruck | 


ee oe (13) | 


ye J. 
(F) 
Dieser Wert ist eine (wegen 8) werkstoffbedingte QuerschnittsgréBe. 
Die Umrechnungsformel zwischen Element- und ee lautet also fir den beider- , 
seits frei a ee toe: Biegeschwinger 


Yn = va aap n- : (14) 


In ahnlicher Weise sind die Formeln fir die wichtigsten Belastungsfalle abate: und in Tafel 1 
gusammengestellt. 


7. QuerschnittseinfluB. Die QuerschnittskenngréBe 74 kann nach (13) fir jeden: beliebigen || 
Querschnitt ermittelt werden. Fiir die beiden wichtigsten Falle, den Rechteck- und den Kreis- 
-ringquerschnitt sei die ‘Ableitung kurz skizziert. 
Fir den Rechteckquerschnitt (Abb. 4a) ist der Randfaserahstand rea und das axiale 
Flachentragheitsmoment J=dh3/12. Fur das Integral gilt 


dhs+6 | 

3 Be WOS TS eb, Bandanas | 
‘ es F=2 25+B(3+B) © | 
Damit erhalt man durch Einsetzen in (13) i 

den Wert | } || 

3 ses | 

Marae © ge (1S) 


Beim Kreisringquerschnitt (Abb. 4b) _ || 
empfiehlt es sich, von den kartesischen | 
‘Koordinaten y, z auf Polarkoordinaten r, Y 
uberzugehen. Der Randfaserabstand ist 

 ¥r=Tq, das axiale F Pace aiagacte 
ist 


“Abb. ia: a) RasdiebLguasschiies mit kartesischen Koordinaten y, z. TET" ; 
» (ee mit Polarkoordinaten T, Q. A ; J= a3 a at —o*) mit e= 4 =T/ 2 ie 


Ai ‘ 


sta fiir das Integral ergibt sich damit 


~ 


4+68 (1—o4+8 
| freer g TattP(1—et*8) a Lf Gin gh dy 
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a 


“a - - % 
_ Der Ausdruck 2 oe gp) +® dp=K ist ein Zahlenwert, der nur von f abhangt und dessen 


Ermittlung eit grundsatzlichen Schwierigkeiten bereitet. Die endgiltige Formel far den Kreis- 
ringquerschnitt lautet dann 


4K 1—ot+8 
~ g(4+p) 1—o | (16) 


Die Formel fiir den Vollkreisquerschnitt erhalt man als Sonderfall fur o=0. 
Siehe Tafel 1. 


8. Anwendung der Tafel. In der Tafel 1 sind die Umrechnungsformeln zwischen Nenn- 
dampfung und Elementdampfung fir die wichtigsten Belastungsfalle und Querschnittsformen 
zusammengestellt. Man sieht, daB bei Rechteck- und Vollkreisquerschnitt die Abmessungen 
keine Rolle spielen und daB8 beim Kreisringquerschnitt nur das Halbmesserverhiltnis @ in die 
Formel eingeht. Ein Einflu8 der Balkenlange I ist nicht vorhanden. 

Nichtsdestoweniger schlagt sich der Einflu8 von Querschnittsform und Belastungsfall in 

den Umrechnungsformeln nieder. ' 

An einem Beispiel mége nun gezeigt werden, was die Tafel leistet und in welcher Weise sie 
angewendet werden kann. 

Es ist zunachst notwendig, fir einen gegebenen Werkstoff die Dampfungskurve y=« é 
aufzunehmen. Dies mége z. B. durch Dauerversuche mit zylindrischen Probestaben an einer 
Prifmaschine, bei der das Biegemoment tber die MeBlange konstant bleibt, geschehen. Man 
nimmt so experimentell eine Kurve pvotnreis=Y,(ER) auf. Durch Umkehrung der entsprechenden 
Formel in Tafel 1 kann man daraus die Elementdampfung des Werkstoffs errechnen: 


— ; re 7 ¢ 
= Mee ne YR= tah (€x) = rae — ae PYollkreis * 


' 


e. Es mége nun ein beiderseits eingespannter Schwinger mit Bingelniasce j in dee Mitte gegeben 
sein; der Balken habe Rechteckquerschnitt und bestehe aus dem eben besprochenen Werkstoff. 
: Dann 148t sich in einfacher Weise aus der Tafel die Formel ablesen, mit deren Hilfe man die 
_ Nenndampfung dieses gegebenen Bauteils aus der Elementdampfung y=« e% errechnen kann: 


i 
‘ 


: 9 
h, » ¥ ie : VRechteck > YR BiB a rane ‘ER (3 +f)? ; 


: Will man nun den Resonanzausschlag dev Schwingers sontittell: so itberfuhrt : man die Werk- 
stoffdampfung, die ja einer Potenz des Ausschlags proportional ist, in eine fiktive geschwindig- 
keitsproportionale Dampfung. Bei dieser Reduktion muf die Darapfangsarbett je Schwingungs- 
periode erhalten bleiben. 

Far die geschwindigkeitsproportionale Dampfung tah sich nach bekannten EF ormeln diese 
_ Dampfungsarbeit zu Bis 


E eae. ist be der Dampfungsfaktor, Q die| Erregerfrequenz und f die Schwingungsamplitude 
aa ‘(Biegepfeil), gemessen an der Stelle der schwingenden. Masse. 


Fa ur ae Peto kdamrtune ergibt sich nach (2) die Dampfungsarbeit zu 
Bere eS UA eis se Oecd th 418) 


t, ermittelt sich nach (9) und ergibt auch je nach Se eees verschiedene Ausdriicke. Auch 


iese- Formeln sind in Tafel af zusammengestellt. 
 Sowohl Vn als auch A sind dort in Abhangigkeit von’ der Groftdehnung § ER angegehen. Der 
usammenhang zwischen GréBtdehnung ép und Schwingungsamplitude f hangt ebenfalls vom 


AN 


ee a S,=b, Qa f?. a foie fe) 3 (17) 


ie. pet ndcravsecncrcet ei in der Schwingungsumkehrlage im Schwinger aufgespeichert “ 


lastungsfall ab und ist ebenfalls i in Tafel 1 zu finden. Es ist nun also méglich, die Dampfungs- es 
r eit s als Funktion der eae roan lnieheny und gwar gilt fir das vorliegende Beispiel _ feos 
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‘und damit nach (18): 
E dhl E dhl 12h \2+8 
ase Sees ES 2 oul 2 
Pemee ae R= a" (ga) ( P ) Jans (19) 


a ae Gleichsetzung der Formeln (17) und (19) ergibt sich dann fir den Dampfungsfaktor 
er Wert 


ak dhl 12h \2+6 
= B 

1222 (3+pP ( P ) ose 
Hiermit kann man nun in die bekannten Differentialgleichungen fiir geschwindigkeitspropor- , 
tional gedampfte Schwingung eingehen und den Resonanzausschlag errechnen. 


(20) 


9. Kritik. Die abgeleiteten Beziehungen gelten unter folgenden Voraussetzungen : 


1. Es gilt die Annahme von Navier und Bernoulli, daB die Querschnitte der Bauteile vor 
und nach der Beanspruchung eben bleiben. Diese Annahme trifft bekanntlich bei Torsion von 
Staben mit nicht kreisférmigem Querschnitt nicht mehr zu, weshalb hierfiir auch keine Formeln 
abgeleitet wurden. 


2. Nicht bericksichtigt wurde eine zusatzliche Dampfung, die durch gegenseitige thermisch- 
elastische Beeinflussung der Elemente bei inhomogener Spannungsverteilung (Warmefluk- 
tuationsdampfung) auftritt. Diese Dampfung ist bei normalen Frequenzen vernachlassigbar und 
erhalt eine gewisse Bedeutung erst bei Frequenzen in der GréSenordnung von 1000 Hz., Es 
sei hier auf die einschlagigen Arbeiten von K. Bennewitz und H. Rétger! sowie von C. Zener? 
verwiesen. 


4 


3. Hinsichtlich der Formeln fiir Biegeschwingungen wurde die Dampfungsarhbeit, die von den 
wechselnden Querkraften hervorgerufen wird, vernachlassigt. Sie ist sehr gering. 


Es kann nun die Frage aufgeworfen werden, ob nicht eine Umrechnung der Dampfungs- : 
kurve p=a* y6* fir Schubspannungen in die Kurve y=« e fir Normalspannungen eines Werk- at 
stoffs méglich ist. In dieser Richtung hat J. Schmidt? Versuche durchgefihrt und fiir verschiedene $ 
Werkstoffe einen Umrechnungsfaktor ermittelt. Dieser Faktor ist jedoch fiir jeden Werkstoff 
ein anderer, und es fehlt eine Begriindung dieses Faktors auf allgemeinen Werkstoffkonstanten. 

Man miBte also bei einem neuen Werkstoff doch einmal sowohl Schub- als auch Normalspannungs- zh 
versuche durchfiihren, um den Faktor zu finden; und.da kommt es auf das gleiche hinaus, wenn 
man auf die beiden Dampfungskurven direkt zurickgreift. 


10. Zusammenfassung. Es wurde hervorgehoben, daB man zwei Kennwerte zur Darstellung 
_ der Werkstoffdampfung unterscheiden mu. Der eine, die Nenndampfung. ist eine Funktion von 
b: Werkstoff und Bauteil; der andere, die Elementdampfung, ist eine reine WerkstoffgréBe. Aller- 
dings hangen beide von der Beanspruchung ab. Man kann diese Abhangigkeit fir einen gegebenen 
Werkstoff im allgemeinen durch eine Potenzfunktion darstellen. ee 
: Wahrend bei Zug-Druck Nenndampfung und Elementdampfung tbereinstimmen, findet man the 
bei Torsion, daB die Nenndampfung im allgemeinen von Querschnittsform und -gréBe abhangt. 
Bei Biegung ist die Nenndampfung auBer yom Querschnitt auch noch vom Belastungsfall (Auf- 
lagerung, Einspannung, Belastung) abhangig. 
Es wurden Umrechnungsformeln abgeleitet, die es gestatten, aus der Elementdampfung 
_ des Werkstoffs fir die einfachen Biegeschwinger die Nenndampfung des Bauteils zu ermitteln. 
-SchlieBlich wurde gezeigt, ‘wie man daraus die Dampfungsarbeit errechnen und nach Reduktion 
auf eine fiktive geschwindigkeitsproportionale Dampfung den Resonanzausschlag unter Berick- 
‘sichtigung der Werkstoffdampfung ermitteln kann. — 


1 K, Bennewitz u. H. Rétger, Phys. Z. 37 (1936) S. 578; 39 (1938) S. 835; Z. techn. Phys. 19 (1938) S. 521. 
by: 2 C. Zener, Phys. Rev. 52 (1937) S, 230; 53 (1938) S. 90; Proc, Phys. Soc. London 52 (1940) S. 152.) 
8 J, Schmidt, Die Dampfungsfahigkeit von Eisen- und Nichteisenmetallen. Mitt. Wohler Inst. 9 (1931). 


ow 


/ 


bene (Eingegangen am 19. August 1948.) ; 
Anschrift des Verfassers: Dr.-Ing. Giinter Denkhaus, (17a) Karlsruhe, KanonierstraBe iy 
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Die Bahn des Momentanpols und die Kardanlage. 
Von W. Meyer zur Capellen. ‘ 


1. Vorbemerkung. Der Begriff der Kardanlage dirfte zuerst von Rauh' benutzt worden 
sein und wurde von Alt? naher untersucht. Da hierbei verschiedene, in friheren Darstellungen®, 4 
gebrachte Gesichtspunkte tbersehen wurden und einige wesentliche Erganzungen zweckmabig 
sind, sei das Problem nochmals beleuchtet. 


2. Die der Kardanbewegung momentan ahnliche Bewegung. Bei der Kardanbewegun eau 
rollt bekanntlich ein Kreis (Gangpolbahn p) im Inneren eines doppelt so groBen festen Kreises 
(Rastpolbahn x) ab, wobei getrieblich auch zur Erzeugung der gleichschenklige Schubkurbel- 
trieb oder die Verschiebung einer Strecke mit ihren Endpunkten auf je einer Geraden gewahlt __ 
werden kann. Hierbei ist p identisch mit dem Wendekreis, dessen augenblicklich mit dem 
bewegten System verbundene Punkte Wendepunkte ihrer Bahn, hier aber doppelt zahlende 
Strecken, und zwar Durchmesser von gz beschreiben. 

Hier beschreibt also auch der Momentanpol P, d.h. der Berthrungspunkt der Polkurven, 
als Systempunkt eine Bahnstelle mit unendlich groBem Krimmungskreis, wahrend im all- 
gemeinen der Momentanpol eine Bahnstelle mit dem Kriimmungshalbmesser Null, eine gen . 
woéhnliche Spitze, beschreibt. Von einer momentan Ahnlich gearteten Bewegung wird | 


man verlangen miissen, daB der Momentanpol nicht gerade eine Bahnstelle mit unendlich groBem, | 


immerhin mit endlichem Krimmungshalbmesser beschreibt, d.h., daB die von P be- 
schriebene Bahnkurve eine Schnabelspitze besitzt. : 
Wie R. Miiller gezeigt hat, ergibt sich fir den Kriimmungshalbmesser der bei beliebiger 
ebener Bewegung vom Momentanpol beschriebenen Bahnstelle der Wert 


(2b) | 


ae db4+2dtid$42d*t 
Bir dig eae eaten ce — Q) 
worin ds das Bogenelement der festen Polbahn x, d8, d?+d2%,... die unendlich kleinen Dreh- 
ungen um die Momentanpole P, Q,... und dt der Kontingenzwinkel der Polbahn x bedeuten. 
Im allgemeinen wird hiernach R=0; wenn jedoch dt=d@ oder : 

Ras tae ee . ; (2a) 
wird (Striche bedeuten Ableitungen nach der Bogenlange s), so wird R nicht gleich Null, sondern _ 
hat den im allgemeinen endlichen Wert ; ; - 

: pele 
Mo Vics en eee ue 


Wie R. Miiller>*° bereits gezeigt hat, berihrt der Krimmungskreis hierbei finfpunktig, d.h., 
' P ist ein Burmesterscher Punkt*. Der Krimmungsmittelpunkt liegt auf der Polbahntangente, 
_und eine solche Getriebestellung liegt beim Gelenkviereck, wie ohne die vorstehenden Formeln 
leicht zu erkennen ist, in den Umkehrlagen, Abb. la, b, vor, da dann der Punkt B=P den 


auBersten Punkt seiner Kreisbahn erreicht hat. AuBerdem ist eine solche Getriebestellung dann |} 


vorhanden, wie bereits R. Miiller gezeigt hat und wie unten elementargeometrisch nachgewiesen 
_ wird, wenn der FuBpunkt des von P auf den Steg A, B, gefallten Lotes mit dem Schnittpunkt H_ |} 
von Steg und Koppel zusammenfalit (vgl..Abb. 2). Der Kriimmungshalbmesser hat dann den 


Wert jas aetae Pica = | pa | 
~ “cos(B—a) cosy —sin(atp)siny crn ee 


Nun sind zwei bemerkenswerte Sonderfille zu nennen®: 


; 1K. Rauh, H. Marks, W. Biindgens, K. Otto, Praktische Getriebetechnik Bd. 2, ‘Kardanbewegung und 
_ Koppelbewegung. VDI-Verlag 1948. #4 Rates ie veep eR Wk 
2 H. Alt, Ing.-Arch. 14 (1944) 319. eae rhea Pe Pa Garpewietire| |, 
3 a) Reinhold Miiller, Einfiihrung in die theoretische Kinematik. Berlin 1932. — b) Ders., Z. Math. Phys. 48 a) 
(1902) 224. — c) Ders., Z. Math. Phys. 42 (1897) 269. Php aeis ah s 
‘ W. Meyer zur Capellen, Mtsh. Math. Phys. 41 (1934) 285. As errs si ae At he foo 
Q Von Durchschlags- und Verzweigungslagen, ebenso von Lagen mit unendlich fernem Pol soll abgesehen © 
werden, : Met ha fy Wo 


' 


{ 


= 
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a) Der Krimmungshalbmesser wird unendlich gro; dann beschreibt P. eine gewohnliche 
Spitze, aber mit finfpunktig berihrender Tangente, was der oben skizzierten Kardanbewegung 
naher kommen dirfte. Nach (2) tritt dies ein, wenn neben 1’ = 9" noch 29” — t’’ ist (vgl. unten), 
und beim ebenen Gelenkviereck erfordert dies das Verschwinden des Nenners in (3) (vgl. unten). 


Abb. 2. Allgemeine Stellung 
Abb. 1, Umkehrstellungen beim Gelenkviereck als Kardanlagen. des Gelenkvierecks in der Kardanlage. 


b) Der Kriimmungskreis hat nicht nur fiinf, sondern sechs unendlich benachbarte Punkte 
mit der Kurve gemeinsam (die Héchstzahl bei der Koppelkurve.als trizirkularer Kurve sechster 
Ordnung). Dann muff R einen stationéren Wert haben. Ein gewéhnlicher Kriimmungskreis 
vom Kriimmungshalbmesser ry berihrt dreipunktig, eine vierpunktige Berihrung erfordert 
dr =0, eine fiinfpunktige d?r,=0 und eine sechspunktige zudem d37,=0. Da hier in unserem 
Zusammenhang bereits eine finfpunktige Berithrung vorliegt, muB zudem d?R=0 sein, und 
dies fihrt nach (1) mit t’=# auf die Bedingung 


t’=0, (4) 
was 
3K 
pet (5) 
liefert!- und nach (3) beim Gelenkviereck 
. tg ytg(a+6)=—1 (6) 
fordert; der Kriimmungshalbmesser R hat nach (3) den Wert 
Sj eae (7) 


. 2 cosa cosf ’ 
d.h., der Kriimmungskreis ist identisch mit dem Umkreis des Dreiecks Aj By Po, worauf unten 
noch naher eingegangen wird. 


3. Die Polkurven in der Kardanlage. Bis jetzt haben wir von der , Ahnlichkeit’* mit einer 


 Kardanbewegung im betrachteten Augenblick, d..h. von einer ,/Kardanlage“ nur gesagt, daB sie 


durch eine endliche Kriimmung oder durch die Kriimmung Null der vom Momentanpol be- 


i schriebenen Bahnstelle charakterisiert sei. Eine weitere Abnlichkeit ergibt sich durch die Be- 


trachtung der Polkurven und vor allem der Kreispunkt- und Mittelpunktkurven. ; 
Die Krimmungshalbmesser 0,= PM, der festen Polbahn x und rp = PM, der beweglichen 


 Polkurve p sind bekanntlich? durch 


a= 1 Ty 
Ox aes ia 3 Tp = Wat (8) 
gegeben und durch die Euler-Savarysche Gleichung 
DOL 
Tp On Puri 


verknipft, wobei die 1/ gleich dem Durchmesser des Wendekreises ist. Die Punkte M,, Mp 
liegen mit W auf der gleichen Seite der Polbahntangente, wenn t' positiv ist (vgl. Abb. 3a far 


die allgemeine Lage). 


Z a Glchg. (4) ist in der Glchg. (16a) der in FuBnote 4 von S. 308 zitierten Arbeit falsch angegeben, ohne 
-» daB die weiteren Ergebnisse dort davon beriihrt werden. 


2 Vgi. FuBnote 3 von S.1 Arbeit a). : 


Rog Serre ots ee Sta 
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Die Euler-Savarysche Gleichung = - 7 = -- vermittelt die Beziehung zwischen den Abstan- 
den PA = rund PA, = 9 des Bahnpunktes A und seines Krimmungsmittelpunktes 4, vom Momen- 
tanpol P; sie bildet die A-Ebene auf die A,-Ebene ab, und wenn die den Wendepol W, d. h. die 
positive Polbahnnormale enthaltende Halbebene, als obere Halbebene bezeichnet wird, so gilt, 
wenn die in die untere Halbebene weisenden Radienvektoren r, @ das negative Vorzeichen er- 
halten, die folgende Zuordnung, die zwar bekannt! ist, aber in Hinblick auf das folgende noch-- 
mals hervorgehoben werden muB: Der Halbwendekreis bildet sich in den Wendekreis, das zwi- 
schen Halbwendekreis und Wendekreis gelegene Gebiet in die obere Halbebene, die obere Halb- 
ebene auGerhalb des Wendekreises in die untere Halbebene ohne das Gebiet des Rickkehrkreises 
und die untere Halbebene in das vom Rickkehrkreis umschlossene Gebiet ab. 


Abb. 3. Die Mittelpunkte der Polkurven und die Kriimmungskreise der Kreispunkt- und Mittelpunktkurven 
a) in allgemeiner Getriebelage, b) in der Kardanlage. 


Nehmen wir nun zu den Gleichungen (8) noch die Bedingung (2a) hinzu, so heifbt dies aber, 
daB 0,=dw» und rp=}d, ist, d.h., daB der Kriimmungskreis der beweglichen Polbahn halb so 
groB ist wie der der festen Polbahn, daf dieser jenen umschlieSt und daf ferner der Wendekreis 
Krimmungskreis der beweglichen Polkurve p ist, wie bereits bei R. Miiller angegeben (vgl. 
Abb. 3b fir diesen besonderen Fall). Die erste Tatsache wurde von Rauh als Definition der 
Kardanlage, die zweite zusatzlich von Alt gegeben. 


4. Die Entartung der Kreispunktkurve und die Mittelpunktkurve. Der Ort derjenigen System- 
punkte, die im Augenblick eine Bahnstelle mit vierpunktig beriihrendem Kriimmungskreis be- 
schreiben, ist die Kreispunktkurve k, und der Ort der zugehérigen Krimmungsmittelpunkte, 
d. h. Kreispunktkurve fiir die Umkehrung der Bewegung, ist die Mittelpunktkurve x, beides | 
Fokalkurven dritter Ordnung mit je einem Doppelpunkt in P und der Polbahntangente und jf 
Polbahnnormalen als Tangenten in P. Seien PL=l1, PK=m die Halbmesser der Kriimmungs- ~ 
kreise von k in P (vgl. Abb. 3a) und PL*=I*, m* die entsprechenden Halbmesser fir x, so gilt 
bekanntlich Bea es 

; 3 30 3 : 
Be 20-E 7" 2 me == mn* oe? lis Gl Eg * ; (9) 4 


_ wobei l_und I* wieder durch die Euler-Savarysche Gleichung verkniipft sind. 

Fur den besonderen hier betrachteten Fall t’=& wird aber l= 1/9’ =d, und [*= «0, wahrend 
m bleibt. Dies heiBt zunachst, da® der Wendekreis auch zugleich noch Krimmungskreis der 
Kreispunktkurve in P ist. Wenn aber 1* nach o strebt, so bedeutet dies, daB die Mittelpunkt- 
kurve in den Kreis % und die Polbahntangente t zerfallt®. Von dieser kommt allerdings nur 
ein Punkt in Frage, namlich der Kriimmungsmittelpunkt P, der von P beschriebenen Bahn- 
stelle, wahrend bei der kinematischen Umkehrung alle Punkte der Polbahntangente Bahnstellen 


mit vierpunktig, der Punkt P, aber eine Bahnstelle mit finfpunktig berithrendem Krimmungs- 
kreis beschreiben. 5 v 


> W. Meyer zur Capellen, Z. angew. Math. Mech. 17 (1937) 288; Getriebtechn. 9 (1941) 489. 
2 P ist hiernach auch der Ballsche Punkt der Systemlage.. 
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Diese Entartung der Mittelpunktkurve in den Kreis x und die Polbahntangente ist ge- 
wissermafen der Zentralpunkt des Problems, und von diesem aus kénnen die Zusammenhange 
einfach und tbersichtlich dargestellt werden. 


Geben wir uns (Abb. 4) zwei Punkte Ay, By auf dem angenommenen Kreis x, so liegen die 
zugehérigen Bahnpunkte A, B auf den Strahlen PA), PB. Nun wird doch bei den wblichen 
Konstruktionen! der Kollineationspunkt H so bestimmt, da® der Strahl PB, um den gleichen 
Winkel gedreht wird, mit dem die Polbahntangente in die Richtung PA, gebracht wird; d.h., 
es ist Winkel tPH=«+f/. Andererseits bildet, da M, A, und Mi: B, mit t die Winkel 2« und 
26 bilden, auch die Mittelsenkrechte von A, By (Abb. 5) mit t den Winkel a+, d.h., 
PH steht auch senkrecht A, By. Alle méglichen Punkte 4, B liegen somit auf den durch H 
gezogenen Geraden. Wahlen wir einen Punkt A, so ist B damit bestimmt und damit auch der 


. Di tartete Mittelpunktkurve x und zwei 5 Abb. 5. Ermittlung des Gelenkvierecks, bei dem der Momentanpol P 
Be ea RE tae Av As By, B. eine Bahnstelle mit fiinfpunktig beriihrender Tangente beschreibt. 


- Wendekreisdurchmesser d,, = 1/; geben wir uns diesen, so sind damit A und B bestimmt, wobei 

AB durch H geht, da ja A, A, und B, B, durch die Euler-Savarysche Gleichung miteinander 

2 verknipft sind. a : ce | 
Die Bedingung, daB PH A, ein rechter Winkel ist, ist fir jedes Getriebe in der Kardanlage og 
kennzeichnend und nur eine Folge der Entartung der Mittelpunktkurve (s. jedoch unten). Durch ae 
die Wahl des Wendekreises w, d. h. durch Wahl von ®, und durch Wahl von #”, d.h., durch bya 


P i 

Wahl von m ist auch x bestimmt; erst die Annahme von t”, d.h. der Krimmungsanderung a 
_ der festen Polkurve, bestimmt R bzw. Po, oder die Annahme von P, bestimmt gemaB (2b) den a 
Wert von t”. So kann z. B. A gerade den Scheitel einer Ellipse oder eine zyklische Kurve be- a 
schreiben oder auch eine Koppelkurve irgenduines. Gelenkviereckes ; mur miBten A bzw. A, wp 
auf der zugehérigen Kreis- bzw. Mittelpunktkurve liegen. Ebenso kénnen die Bahnen von A oe 


und B als feste Kurven ausgefiihrt sein (Nocken, Schlitz), aber in der betrachteten Systemlage Re, 
“miissen diese einen’ vierpunktig bertihrenden Kriimmungskreis vom Mittelpunkt Ay bzw. By 
ae nun A,, B, feste Punkte sein, so stellt 4, By BA ein Gelenkviereck in der Kardanlage 
dar, und der dem Dreieck A,B, P umbeschriebene Kreis ist neben der Polbahntangente ¢ die 
_Mittelpunktkurve. Da dann A, B Bahnstellen’ mit 
mehr als finfpunktig berihrendem Krimmungskreis 
* beschreiben, lat sich mit Hilfe der héheren Wendepole, 
wie R. Miiller gezeigt hat, der Wert von 1” angeben 
“und damit die Formel (3) -herleiten. | 


Ps. Sonderlagen. Bevor wir zu Sonderformen der 
-Getriebe wtbergehen, wollen wir, ausgehend vom 
Kreis Hs nochmals die obenerwahnten Sonderlagen 
_betrachten. are 

a) Wir wablen Ay auf x, aber By auf t (Abb. 6). 
: Damit geben wir uns aber (da t ja auch zur Mittelpunkt-. 
-kurve gehért) den Krammungsmittelpunkt Aer By 


“1 Vgl. FuBnote 3a von S. 308. 


Abb. 6. Abb. 4, aber B, auf t. 


i 
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der von P beschriebenen Bahnstelle, d. h. den Wert nach (2b), und es gilt mit dem Wert m (9) | 
1 1 dt” 


m R 34 


WwW noch gewahlt wird oder durch A bestimmt ist, so kann 1” angegeben werden. Sehen | 
te phe ae “ae Punkte an, so erhalten wir ein Gelenkviereck in der Umkehrlage gemab | 
Abb. 2. Dain Abb. 5 die Punkte A und A, auf verschiedenen Seiten der Polbahntangente liegen, | 
mu n nach unten weisen, d.h., K liegt auf dem negativen Teil der Polbahntangente, m ist |} 
negativ, was auch, da 0 > 0, bedeutet, daB 3” <0 ist. Da 3” die Anderung des reziproken | 
Wendekreisdurchmesser angibt, gilt: K liegt auf der positiven (negativen) t-Achse, Je nachdem | 
ob &’>0 (8”< 0) oder ob der Wendekreisdurchmesser abnimmt Rus | | 

b) Wir wahlen A, auf dem Kreis x und lassen By mit K zusammenfallen. Dies bedeutet, da | 
R=m, also nach (2b) und (9), daB t”’=0” wird. Sind Ay, By feste Punkte, so haben wir ein |] 
Gelenkviereck in der soeben betrachteten Umkehrlage, nur da Steg und Koppel aufeinander | 
senkrecht stehen. Reet 


\; babn von P 


Bahn von E 
\ Bahn von 6 ; 


Abb. 7. Ermittlung eines Gelenkvierecks, bei dem der Momentanpol P eine Bahnstelle 
mit sechspunktig beriihrendem Kriimmungskreis beschreibt. 


c) Wir fragen nach dem Gelenkviereck, das bei gegebenem Kreis und gegebenem Steg A, By 
fir R den Wert © liefert, d. h. nach der geometrischen Bedingung fiir das Verschwinden des 
Nenners in (3). Hiernach mu8 mit y=y, als dem gesuchten Grenzwinkel cos (B —a) cos yo= | 
_ sin (« +) sin yy sein: Man fallt (Abb. 5) von Mx, dem Mittelpunkt des Kreises x, also des dem | 
Dreieck A,B, P umbeschriebenen Kreises, vom Halbmesser m/2 die Lote MP auf HP und 
M.T’ auf HA,. Das Lot von H auf P’T” liefert die gewinschte, Richtung. Es ist namlich 


Mi P’=HT’=} msin (B+a) und My T’= P’H=} moos (B—«), also HU=HP’ cos ‘Yo: aber 
auch gleich HT’ sin Yo d.h., die vorstehende Bedingung ist erfillt. oe : 

_ Bei anderer Lage des Strahles y fallen die Projektionen P’’ von P’ und T” von T’ auf den | 
Strahl yp nicht zusammen, die Strecke P’” T” ware ein MaB far den Nenner von R gema (3) | 


und wiirde auch auf eine Konstruktion von R hinweisen, welche anzugeben zu weit faihrt. — 

d) Wir fragen schlieBlich nach den Deutungen der Gleichungen (5) und (6) fir einen sechs- 
punktig bertihrenden Krimmungskreis in der Kardanlage. Die erste laBt im Vergleich mit dem 
Wert m aus (9) erkennen, daB R=} ™m wird, d.h., daB hier die Mittelpunktskurve x oder der 
Kreis tiber A) B,P zugleich der Krimmungskreis ist. Sei nun T der Schnittpunkt des Steges | 


= 
TO Le 1 > 
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Ay By mit t (Abb. 7), so mu8 der die Punkte A, B enthaltende Strahl Tangente in H an den 
Halbkreis iber PT sein!.. Denn man liest ab: y= 90° +(a +8), d.h. tg y=—etg («+f) oder 
tg a tg («+f)=—1, wie nach (6) verlangt wird. ‘ 


\ 


6. Sonderformen des Gelenkvierecks. Wenn auch der Entwurf von Gelenkvierecken in der 
Kardanlage (ohne daB dies Wort allerdings benutzt wurde) und die Sonderfalle dieser Getriebe 
an anderer Stelle? verfolgt wurden, miissen wir auf diese nochmals eingehen: 


a) Geben wir uns wie oben P, K, A), By, so gibt es einen Strahl durch H, der zu PA, parallel 
ist, der also A als unendlich fernen Punkt liefert. Dies heiBt, daB A, auf dem Rickkehrkreis 
liegt. Sehen wir A, und B, als fest an, so entsteht eine exzentrische Schleifkurbel (der Sonderfall 
einer zentrischen ist méglich?; Abb. 8), bei der auch eine sechspunktige Berithrung? (Abb. 7) 
moglich ist. 7 

b) Bis jetzt haben wir den Schubkurbeltrieb zu Mz 
noch nicht erhalten. Bei diesem beschreibt ein Rr 
Punkt, z. B. B, auf w gelegen, eine Gerade durch 
den Wendepol W, da ihr Krimmungsmittel- 
punkt unendlich fern ist. Weil in der Kardanlage 
erstens die Mittelpunktkurve in t und den Kreis x 
vom Halbmesser m/2 entartet, zweitens aber der 
auf % gelegenen Punkt B, unendlich fern liegen 
soll, muB x in eine Gerade, und zwar die Pol- 
bahnnormale entarten, d.h., die gesamte Mittel- 
punktkurve entartet in t,n und die unendlich 
ferne Gerade. Da aber der Kreis % auch Kriim- 
mungskreis fiir die Kreispunktkurve war, entar- 


tet auch diese,. und zwar in n und den Kreis vom ES eons : 
Abb. 8. Die exzentrische Schleifkurbel in der Kardanlage. 


Durchmesser J. Da hier in der Kardanlage fer- 
ner L mit dem Wendekreis zusammenfallt, ist ; 
dieser neben der Polbahnnormalen die Mittelpunktkurve; also missen alle Punkte von wim Augen- 
blick Bahnstellen mit vierpunktig beriihrender Tangente, d. h. Undulationspunkte beschreiben, 
sein. Rechnerisch fordert m=0, daB 0’ =0 wird. ae. : 

In Abb. 9 ist gegeben der Wendekreis w und ist ferner B auf ihm angenommen. Da auch, 
hier immer noch PH senkrecht A, By sein muf, lauft PH parallel WB und HA, demzufolge 
parallel PB. Bei gewahltem A, ist damit H bestimmt, und A ergibt sich als Schnittpunkt von 


ABP. umbeschrieben ist, zugleich das dritte Fokalzentrum C, und den weiteren 
-P durchlaufenen Doppelpunkt enthalt; ferner berithren sich die drei Polkurvenpaare 
lrei Kreispunktkurven?. Abb. 10 zeigt die betreffende Konfiguration mit den Gelenk- 


i i deal 


_P und deren lachpunkt in K liegt. | 
eh (ae nait a Pies soe Caaiaa area 


SOG Ri age 


den hierdurch méglichen Wertepaaren Ap, By ist der fiir die Folge der Punkte 13 By 
punkte der von P auf die Kreistangenten gefallten Lote bestimmt, d.h. durch — 
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Abb. 10. Die Robertssche Konfiguration fiir die Kardanlage. 


- 


die eingetragene Winkelgleichheit bestehen mu8. Beim’ Entwurf 
B, und C, angenommen worden. — : BO ae 


. 


, 


- vierecken 4) By) BA, A,C, C, A), B,C, C, By, wobei A,A =PA, und ‘By B,= PB Sein und” |. 
der Abb, 10 waren P, K, A, 
WN eee ye iis 


l 


|i 
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Die Wendekreisdurchmesser brauchen nicht, wie auch aus Abb. 10 leicht zu erkennen ist, 


ibereinzustimmen; doch miissen, da die Werte m und R tbereinstimmen, gemaB (2b) und (9) 
die Proportionen %: 03:05 = O12 OF Bf = tl orf sty! gelten. 


8. Stationire Werte. Zuletzt besteht noch die Aufgabe, festzustellen, ob irgendeine der 
betrachteten kennzeichnenden GréBen in der Kardanlage einen stationaren Wert hat: 


a) Bei séechspunktig berithrendem Kraimmungskreis in P wird 1’”=0 (Ziff. 3 u.5 d). 
Da v’ die Krimmung der festen Polbahn ist, hat diese einen Extremalwert. 


Abb. 11. Stellungen des Gelenkvierecks fiir stationiren Wendekreisdurchmesser. 


bp) Ist R=m (Ziff. 5b), so ist 7’”=0’. Da aber die Kriimmung der Polbahnen im nachsten 
Augenblick durch t’+7’ds bzw. & +0’ ds+t’+1’'ds gegeben ist, t’ aber gleich ® war, so bleibt 
las Verhaltnis der Kriimmungen das gleiche, ist also das System auch im nachsten 
\ugenblick in einer Kardanlage. 

ce) Rauh1 sagt, daB der Wendekreisdurchmesser in der Kardanlage beim Gelenkviereck 
in Minimum hat, wahrend Alt? feststellt, da dies im allgemeinen nicht dér Fall ist. Die Be- 
ingung fir ein Extremum ist ”’=0, und diese ist (Ziff.6b) nur beim Schubkurbel- 
rieb (und dem daraus abgeleiteten Schleifschieber), aber nicht beim allgemeinen Gelenk- 
iereck (einschlieBlich der Kurbelschleife) méglich. 

- Zur Klarung mu8 jedoch kurz, unabhangig von der Kardanlage, auf die Bedingung 0” =0 
ir ein Extremum in d, eingegangen werden. Nach (9) bedeutet dies, daB lediglich m—> oo geht, 
Iso die Kurve k in die Polbahntangenté t und den Kreis L und die Kurve % in t und den Kreis L* 
ntartet (Abb. 11). Geben wir uns t, W, L und damit L*, so zeigt sich, daB bei festen Punkten 
1 Vel. FuBnote 1 von S. 308. eee nae cs . 

® Vgl. FuBnote 2 von S.308. : : . ‘ 


a 


sity 
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ye B, Gelenkvierecke entstehen, bei denen Steg und Keppel parallel sind, also H mended | 
fern eae oder bei denen, wenn B, auf der Polbahnnormalen angenommen wird (Abb. 11b), | 
die Kurbel A,A auf dieser liegt, oike bei denen eine Strecklage vorliegt, wenn alle vier Punkte | 


auf n fallen. "| 


Durch / und I* sind nach obigem die Mittelpunkte M,, Mp gemaB 


i 
| 


Pes ORE erat FAMERS NTE 
ee eee 


bestimmt oder folgen umgekehrt J und I* aus 


1 ah, 1 1 Laifse 1 
on aes a kad pai a 
Sind z. B. 0, und rp konstant, d.h., liegt eine zyklische Bewegung vor, so sind hiernach I und /* | 
neben d, leicht zu erhalten!. Bei den in Abb. 11 fir 0’=0 herzuleitenden Gelenkvierecken | 
kann sich aber keines in einer Kardanlage befinden; denn dies erfordert ein Zusammenfallen ‘ 
von W und L, so daB eben, wie oben ausgefihrt, nur Schubkurbeltriebe entstehen kénnen, | 
es sei denn, wenn man nur die Bedingung 3” =0, also ein Extremum des Wendekreisdurch- | 
messers als Grundlage einer Kardanlage definiert. 
Die Kardanlage, zu der wir, ausgehend von der Bahn des Momentanpols, gelangten, ist nach | 
vorstehendem durch vier unendlich benachbarte Lagen gekennzeichnet. Zur Bestimmung des | 
a aac ger tone der von P beschriebenen Bahnstelle ist aber noch der Wert von 
” erforderlich, d. h. aber eine finfte unendlich benachbarte Lage. Diese ist z. B. dadurch fest- | 
Belen da wir einem Bahnpunkt einen finfpunktig beriihrenden Krimmungskreis zuschreibens | | 
wie dies beim Gelenkviereck fir die Bahnpunkte A und B der Fall ist. | 


ai Vel. z. B. W. Meyer zur Capellen, Cheeni 4 (1936) Ad. 


(Eingegangen am 10, September 1948.) 
_ Anschrift des Verfassers: Dr.-Ing. Walther Meyer zur Capellen, (22c) Aachen, Maria-Theresia-Allee 22. 
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‘Der homogene, durchlaufende Trager auf unverschieblichen Stiitzen. 


Von J. Fadle. 
1. Einleitung!. Eine Grundaufgabe der Baustatik bildet die Behandlung des durchlaufenden 
elastischen Tragers, der auf n+1 Stitzen unverschieblich, aber frei drehbar gelagert ist. Die 


Stetigkeitsbedingungen fir die Biegewinkel an jeder einzelnen Stitze liefern fiir einen Trager 
mit feldweise konstanter Biegesteifigkeit EJ, die Clapeyronsche Gleichung? in der Form 


M,.,£+2M, GF + Bre fo oe ee Pe) Eat Pa) 


v1 hs a | Hat i 2) Bs Ih 
aan Pyii Px+1(le41—Px+1) (2 Ly 4.1—Px+1) 
Es bedeuten hierin Jew ae 


M,_,, M,, M,,, Stitzmomente tiber den aufeinander folgenden Stiitzen x —1, x, x+-1 (positiv, 
wenn am linken abgeschnittenen Balkenteil entgegen dem Uhrzeigersiny 


drehend); Be |: : 
De Pes Einzellasten im Felde x bzw.- x+1; 
J Jatt Tragheitsmomente im Feld x bzw. x+1; 
Ves dey _ . Sttitzweiten der Felder x bzw. x+1; 
ip Dick _ Abstande der Einzellasten P,, von der Stiitze x —1 bzw. P,,, von der Stiitze x. 


Die Clapeyronsche Gleichung, in die auBer dem Stiitzmoment M, noch die benachbarten Stiitz- 
momente M,_, und M,.,; eingehen, stellt einen Zusammenhang her zwischen je drei aufeinander 
folgenden Stiitzmomenten und den Belastungen der zur Stiitze x benachbarten Felder und wird 
daher auch als Dreimomentengleichung bezeichnet. Fir einen an den Endstitzen frei aufgelager- 
ten Trager sind die Randstiitzmomente M, und M, gleich Null und es folgen n —1 Dreimomenten- 
gleichungen, die zur Berechnung der n —1 statisch unbestimmten Stiitzmomente M,, M2,...Mn-1 
ausreichen. _ ; 3 
Fir den Fall eingespannter Endstiitzen treten zwei weitere Gleichungen hinzu 


2M. + M, = P,l,-4* (1 T) (2 i) 


Ban | 7 Le Pn Pn Pn : 
Metta rie B) 2), 
welche die Bedingung fir die horizontalen Biegewinkel an den Endstiitzen aussprechen. Es | 
folgen somit n+1 Gleichungen zur Berechnung der n-+1 Stiitzmomente My, M,,... Mz. 


2. Die allgemeine Lésung der Differenzengleichung des homogenen, durchlaufenden Triigers. 
Ein durchlaufender Trager mit gleichen Stiitzweiten | und gleicher Belastung aller Felder wird 
als homogener durchlaufender Trager bezeichnet, fir den\ die Clapeyronsche Drei- 
momentengleichung in die einfache Form tbergeht es 


—— 3 P taea eh far Einzellasten. 


a ie 
a 4 


po | (1) 
|~ 71 Sé(l—-2 q(& dé far Streckenlasten. 
; 0 ¥ ¢ 1 : r) 


en te (uf 
ome, ‘ j 


= ci 4 es) 


>t deem oe : See eal 9 Se oe Z ‘ 3 : ; : 
das Glied auf der rechten Seite den Einflu8 der Belastung durch Einzelkrafte oder Strecken- 
asten darstellt, wird es sinngema® als Belastungsglied B bezeichnet. Gleichung (1) ist eine in- 
homogene Differenzengleichung mit konstanten Koeffizienten, deren Belastungsglied fir den 
vorliegenden Fall gleicher Feldbelastungen eine Konstante ist. Ihre allgemeine Lésung wurde 
schon von Clebsch® durch den Ansatz neo pat 


re, Hs : x : te Z me ee sf Ms ay Sis IES oe ae Cr* 7 Asin : Fake ¢ (2) 


4 


oe diesem ‘Aufsatz werdanke deh ‘Herrn Professor Péschl, der in seiner Vorlesung tuber 
es auf die in der bisher bekannten Literatur nicht erschépfende_ Behandlung des 


er, $.45, Stuttgart 1946. ns 
ester Kérper, S. 394. Leipzig 1862. Bleich-Melan, Die gewohn- 
en der Baustatik, S, 167. Berlin 1927. — aan 


80 Berlin 1948 — Biezeno-Grammel, Technische Dynamik, S. 231, Berlin 1939 Ss 


ae 
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angegeben. Geht man mit dem Ansatz (2) in (1) ein, so erhalt man nach einfacher Umformung |} 
die charakteristische Gleichung 


rP+4r+1=0 (3) |] 


iim? 8, Wa-2—VS, irl. io (4) | 
Die allgemeine Lésung von (1) lautet nun mit den beiden Auflésungskonstanten C; und C, und 
wenn noch r,=r gesetzt wird : 


mit den beiden Wurzeln 


Mak Bie Cia mer) (5) | 


wobei Ci und Cj aus den beiden Randauflagerbedingungen zu bestimmen sind. 


3. Rekursionsformeln und Additionstheoreme konjugierter Zahlen zur Darstellung der Po- | 
tenzen von r™ und r—™. Setzt man allgemein fir ganze Zahlen m | 


rm (24 Y3)" = am + V3 bm, r-™ = (—2+ V3)" =(—2- J 3)" = on — YB bm,” (6) 


so wird 
: | ge SS A I de ce ee (6a) 
Nach Multiplikation der charakteristischen Gleichung (3) mit aes kommt weiter in Zusammen- 
hang mit (6) 

rm—2 4. Apel pm — (a, 9 +4 amit os vs V3 (Om—2 + 4 Om—1 dba) = 0; 
woraus sich die Rekursionsformel zur Berechnung von am und bm ergibt 


adm = —(4 Om—1 + Gm—2) > bm = —(4 Bm—1 + bm—2) » (6b) 


| 


Nun ist 
ro =] Sat 3b), dubs a = 1, b= 05 
i= —24- (3 = a+ 3b): d. hea, = —2,6, = 1, 


= (B41) =+7, b= -(+ 440) ——4, usw. 
In der folgenden Tabelle (6c) sind die Werte fir am und bm bis m=12 zusammengestellt. 


und es kommt 


m am | bin 
0 Bat 0 : 
1 312 sa | 
2 4.4 ae 
3 SG) 415 
4 Oy rs: — 56 
5 <"362-> + 209 . -. (6c) |] 
6 +1351 — 780 ca 
1 — 5 042 i ee OTT 
8 + 18817 ' — 10 864 
9 ~ 10,226 . + 40545 
10 + 262 087 . ¢ — 151 316 
2 ook — 978 122 ' + 564719 
12 7. + 3659401. — 2107560 
Aus (6c) kénnen wir sofort folgendes Vorzeichengesetz feststellen: nie 
m= (1) fan], b= (1 fbf, (6a) 


Aus der Definitionsgleichung (3a) ergeben sich fiir zwei vemaniedene Potenzexponenten m 
und n folgende Zusammenhange 


a ft i) (x ae r= (rmtr pom zB) + las Te 
eee 2 Gm Gn = Onin + Gm—n : 6 a bn = = ae ' 

(rm ab Tn”) (r” a ro") = (wire —T gees Ee (rm—? — rere), 

% am bn =bmin =) bane 2 On bm = = Onan + baxt . Z f 


\ 
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Durch geeignete Zusammenfassung folgen schlieBlich die Additionstheoreme der konjugierten 
‘Zahlen in der Form 


Om+n = Om An + 3 bm ba, Bmtn = bm an + bn Om. (7) 
Im besonderen wird fiir n=1 . 
Z Omti+2dm=+3bm, bmi t2bn = tam. (7a) 


Weiter leiten wir noch folgende Zusammenhange fir drei ganze Zahlen m, n und k her: 


GQ, Om — 3 bn bn, = Ok An_ms 


bes: bm a bn bm—k = by On, ’ 
| : 


Seay Qm — An bea, = ak bse a 


Mit den entwickelten Formeln lassen sich nunmehr Stiitzmomente und Auflagerkrafte fiir 


die drei méglichen Randauflagerungen bei verschiedenen Belastungsarten besonders anschaulich 
ableiten. 


4. Der homogene durchlaufende Triger. a) Freie Auflagerung der Endstitzen: 


M,=0, Mn=0. Es erweist sich zweckmafig, Gleichung (5) durch geeignete Umformung auf die 
konjugierten Zahlen a, und b, zurickzufiihren. Zu diesem Zweck setzen wir : 


B C’,+C’ C’,—C’ B 
Meat SE pr) ¢ Se 9) = 2 4 Cae ¢ Cbs. (8) 
_C, und C, folgen aus den Randbedingungen fir x=0 und x=n 
x=0: M,=0=2+G,, q= 24; 


B B B 1—a,° 
| bpd a Mr = 0 = = — an + Cy bn, ee ae 
Unter Beachtung von (7) kommt schlieBlich g 
G Pen B. by Ontee Bb OW bab, ook 
M, = (1—ae— 5, bs) = => "i Ee ORB ke () 


Durch Vertauschung von x mit n —x geht (9) in sich selbst uber, woraus die Symmetriebedingung 


M,_.—M: folgt. 


b) Starre Einspannung beider Endstitzen. Bezeichnen wir mit B, und B, das 
Belastungsglied fur ein linkes bzw. fir ein rechtes Auflager, wobei B,+ B,=B ist, so lauten 


bekanntlich die Bedingungen fir eine horizontale Tangente der Biegelinie an den EKinspann- ~ 


stellen 
. | a 
x= 0: 2M, + My = By= + G, (2 4)-+ 0) + C, (2 by +b) = + 3C, by + Cy ay, 


1 as 
: C,= ~y (B,—B,); : 
d 1 
x=n: M,14+2 Mn = B,= ++, (an—1 + 2 dn) + > (By —Bz) (Bn—1+ 2 bn) = 
L a; B : 


x n—l 
(= 2-366,-F(B,-B) a, C=-S—(B,-B). 


| 
P 


Fahrt man C, und C, in (8) ein, so kommt 


B, BB: fF ; byt @—On~x by Gxt Gn—s , 
Mz = A a (dn —3 bn bs) — ae] aa TUE b= Paar § eat (10) 
Wertauschen wir in (10) x mit n—x, so gehen die Beiwerte von B, in diejenigen von B, iiber, 
dwh. es liegt Polarsymmetrie vor, wie es auch sein muB, wenn wir die Zahlung nicht von links, 
‘sondern von rechts beginnen. = . : 
Weiter erkennen wir die Dualitat der konjugierten Zahlen a, und b,; denn die bz und be, 
sind ‘wegen by=0 charakteristisch fir die freie Randstitzenlagerung (M,=M,=0), die a, und 
Rie: hingegen sind wegen 2 a,-+-a,=0 charakteristisch fir die starre Einspannung (2 M)+ M,=B,, 


— 


M,-1+2 Mn=B,). Diese Dualitat wird besonders schén bei der Behandlung des Tragers mit 
gleichen Feldlangen und einem einzigen, beliebig belasteten F eld in Erscheinung treten. 

c) Ghiniselte Randstitzenlagerung: erste Stiitze starr eingespannt, letzte 

“Stitze frei au fgelag ert. Da hier wegen der Einspannung der ersten Stitze dieselbe Rand- 

_bedingung wie in Fall 6) gilt, abernehmen wir sofort das dortige Ergebnis mit C,=} (By, — B,). 


! 
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Die freie Auflagerung der letzten Stiitze hingegen erfordert die ee 
ee arene My = 0 = +, n+ > (B,—Bs) bn» C= (By By 1+ By) = 


Mit diesen Werten fir C, und C, wird das Sttitzmoment : 
bn dy—by Gn An —A—3b 


Gy —Ax+3 bp—x 
2an Tape iO aes 


6a, 


Bey 


Me (B+ B,) (B= 8) 5 


sae In den folgenden Tabellen 1 bis 3 sind die nach a) bis c) ermittelten Stiitzmomente fiir 
homogene Trager bis zu 12 Feldern in Form von Zahlenpyramiden zusammengestellt. , 


Mio bbb 


1, Tabelle der Stiitzmomente Bees 6b n* fiir freie Auflagerung beider Endstiitzen. — 
n t 
: ! 
n=1| 0 0 
! n=2|) 0 . 
n=3 | 0 
3 
4 
| == 15 | 0 TE = 
11 8 eee | 
eae ‘ 62 52 
15 11 i 
n mand 0 ma | 7 0 3 
41 30 33 . 
eke 194 1v4 v4 Be 
() 56 41 45 56 
; =n AL kha 
. | o | ves | tes | B65 265 ee 
eee ee 163,-|, 342 ||) 198 120 Le eae 
7244. | 124 T2Z4 724 724 : 
Seat 209 | 153 | 168 | 164 153 209 te Aaa 
989 989e | 93y. 989 Yod You 
ea 571 | 418 | 459 448 451 451 | 448 459 | 418 | 571 |- ha 
- 2702 | 2702 | 2702 | 2702 | 2702 2702 | 2702 | 2702 | 2702 | 2702 ene 
; ; : , 3 | 
d) Ermittlung der Auflagerkrafte. Denken 1 wir uns in bekannter Weise nach Abb. 1 zwel | | 
_benachbarte Felder x und x+1 durch Schnitte aber den Stitzen in Einzelfelder aufgelést, wobei | | 
die Stitzmomente als auBere Momente einzufihren sind, so erhalten wir fir die Anteile AG 
von pueld x und AZ von F eld eee zur Auflagerkraft A, ke & 
2 
| <Mea 1 a 
yy Zi —1 i eS 
G4 eG Cs yf enpae. | 
‘i ay 
era eae (ay 
sae M,—Me.1. > 1 > Ne he eee | 7 
J , AVON ee ax+1 j ee) sa & % 
Aj = SS aes alé) ae ae 
: my , eN@) he Gel pees J es ‘ ae | 
pie: eat | | Sahoo a _M, ns 
Le hs Then! Bestirstanns der Abtiapeseriits HETREY: (Waeesayy a Ag oe eee ogra x+1 / EA 
RA eo des homogenen Trigers. Az = Spies ee ES ue fava. ce - 


~ 
oobi 


" SchlieBlich kommt mit (1). lay einfacher Umformung las ORIN ak Shlae e: ; 


LoL te =318-4 82) (8) d= 6 r Head flee 


} i : , a Sf - 


te 


¥ 
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2--_90T8 2 9018 z- 9018 2--_20T8 2 90T8 29018 Zz, 9018 z—. 90T8 2 _90T8 z,_90T8 / 
ae ate + air ae AG i a ote = 
a T69€ eo PSL = 61ST s 90ET # S9&T ei 6ST gq 968T be €8tr 8 8L6T Z. 686 gl =u 
T_90T8 T_90T8 T,_90T8 T9018 T~_90T8 T_90T8 Tt 9018 12018 T,_ 90T8 T_90T8 
oa sé 
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Ge a a ee 


Pir die erste und letzte Stiitzkraft A, und A, erhalten wir entsprechend den zwei Lagerungs- 
méglichkeiten als freie Auflagerung oder starre Einspannung zwei verschiedene Ausdricke. 


Zunachst folgt mit A aus (12) fir x=0 
1 


M,— M. 1 
A= 5 taf @ £) q(é) dé, 
und wir erhalten weiter fir 1 
a) freie Lagerung: M,=0, d.h. Ay =— - ey C=) q(é) dé, 
rt 


I 


b) starre Einspannung: 2M,+M,=8,= ai E (I— €) (21—&) q(&) dé, 
0 


u 1 
doh. Ag= 37% — 4247 f (8 of) d8=3 4% f U—2 a8) dé. 
0 0 


Weiter kommt mit A‘, aus (12) fir x=n 


} M,,—Mn-1 1 
; CE ae aa J &ag) a8. 
und wir erhalten weiter fiir 1 
Bre é \ My_1 1 
a) freie Lagerung: M,=0, d.h. 4n = — ; ~ + [ &4() dé, 
; 0 Z 


U \ 


b) starre Einspannung: M,-1+ 2 Mn = B, = = fea-8 (1+-&) q(é) dé, 
mes 


J I 
dh. An = 3" <p pf baled 3-7 toe | fa) aE. 
: 3 


Wir erhalten somit folgendes Schema, in welches die M, dem betreffenden Lagerungsfall ent- 
sprechend einzufihren sind. 


1 


c oak 1 u 
freie Auf- M. 1 : iM cal . : Maa. ok 
oe. —7 tT] owe co md @=ShE+ 3a) a(6) dé mist +7 fea@as 
4 i ; epee 1 
starre Ein- M, Tt Max lest M, 1 
vpanaung | 2p te f C—O 9Ods| 62 + ze | [P+0-9 aed | aes taf se? 
i ee 0 


e) Beispiel. Als Anwendungsbeispiel wahlen wir einen homogenen Trager auf 7 Stiitzen 


(n=6) mit einer trapezférmigen Feldbelastung gema8 q(£)=qi,+ i &. Eine einfache Rech- 


aS 


: nung liefert folgende Belastungsintegrale: 


Af c-aqat=Featm. — +f fa d= GF (at2m!, 
= Be ek atee 0 
S Af )8q(é) dé = a (4ata)l, 1 fe (8) dE = so (at4q)l, 
ie B , 4 eyed oo 2D = qr)"; Bs q 20. > 


Z 1 : Vesa : 
1 f sie 4 3eyqeyde=z f (92+ 8] a8) d= Gata) 
OF, ; Oe 


See : 
By =z f El-8) 21-8 8) dé = Batt HP, 
; 40° : 


: ; See vu ; 
B= fed-n rode as=Zatsn?, 
5 a 0 ars? post 


: ; Eee [ee 
B=B+ Ba Zeta). 


; 
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Aus den Tabellen 1 bis 3 entnehmen wir 
oe = A . a nun fir n=6 mit den berechneten Belastungs- 
= S 2 E | gliedern B, und B, die Stitzmomente und k6on- 
° rd BY rn SPS, eo. 1 418 nen damit auch die Auflagerkrafte ermitteln. 
S ee | ss In nebenstehender Tabelle sind die Ergebnisse fiir 
E " SF S die drei Randauflagerungen zusammengestellt. 
5. Triger gleicher Feldlangen mit einer ein- 
Q af ue x = zigen Belastung im i-ten Feld. Nach dem Super- 
a Bee a S s a positionsgesetz kann eine beliebige Belastung 
my = 8 S 2 ° &\s Blo des Tragers gleicher Feldlangen in n Unter- 
= Sah es ela {st z a 5 belastungsfalle aufgespalten werden, wobei in 
al2 3 |x 2 & y a jedem Unterbelastungsfall nur ein einziges Feld 
RO ve & a mit seiner Feldbelastung belastet ist. Die allge- 
Bs meine Belastung aller n Felder kann dann durch 

: ~ Uberlagerung aller n Unterbelastungsfalle zusam- 
a eared Wes oe < = mengesetzt werden. 
= A = | eebe ei 2|8le Bezeichnen wir mit i das einzige belastete 
S&S & ra a E 3 +\8 es 5 Feld auf den Stiitzen i —1 und 7, so kénnen wir fiir 

@|8 2/8 a 2 = a die Stiitzmomente der unbelasteten Felder 
a ~ BR 2 folgenden Ansatz nach (8) mit vier Auflésungs- 
ae konstanten C,, C,, Cj, C; machen (siehe Abb. 2): 
a x a OS * Sral My =) a, 4G 6.5 | (14) 
oi a a ey eS 8 | iSxZn: MA = G, 424 C, t. 
PE GaN eliget cae, Sp eat 82 18 ; 
Ss aS) aS a i 4Je/ 4/8 Die vier Konstanten aus (14) reichen aus, um 
2/2 a\8 3 ~*~ Jo Ss die zwei Auflagerbedingungen an den beiden 
nN 62 | 0 aad ; g o is ‘ ‘ pean = 

s 2 Endstitzen und die zwei Ubergangsbedingungen 

an an den Stiitzen i—1 und i zu_befriedigen. 

i Wir haben auch hier wieder gemaB Abschnitt 4a) 
ee De fas = a $ bis c) drei charakteristische Lagerungsfalle zu 
s $ S aS 8), |S. unterscheiden, die sich mit Hilfe der dualen Zah- 
S i pa 5 af Z e\s | ais len a, und bz leicht und tbersichtlich erledigen 

fr eae 2 B= So Se | ss lassen. 
“8/58 |g] |8| [8] |g Zunichst befriedi ir die far alle drei 
a <= 8| | 8 s gen wir die fir alle drei 
aha aos Lagerungsfalle gleichen Ubergangsbedingungen 

i ae durch Aufstellung der Dreimomentengleichung 
eS Ae he ~ a ie fir die Stiitzen i—1 und i und erhalten fir 
as aS S Ss hee ; 
= = Ble | Sle § 3 5 g x=t—1: Were Me™ (15) 
ra a S| Sr ts f &|5 C—t- M;_\ +4 M;+ Mi, = By. ; : 

SU A Se 8 % Nun miissen aber wegen der fir unbelas- |} 

i _  tete Felder geltenden Dreimomentengleichung _ |] 
es dat Rar ~) to dl Ma-2+4 M,_1+ M,=0 mit den nicht mehr vor- ; 
set & = SZ a kommenden Stiitzmomenten (Mj) und (Mj_,). ; 
pees Lae Si = 3\518le folgende Gleichungen erfillt sein: Seeane 
o | 4/81 5/8 | SIE | Sis | SIS ve ae ~ 
j s S|? 5 ale % § a eso Mt 2 M;_) + (Mi) cen gee oe a ; 
S| 13) [8] 18) [2 (Mia) + 4 Mi ML = 0. | ae 

Sane us ee : Féhrt man (15a) in (15) ein, so kommt a | 
Pt el ae 

A | Sag ih Oe eer My, — (M4) = Beater hee ‘ 

a = 7) — (Ch—Cy) ai + (Cz —C,) bs = B,, | 


v 


gc | (CG -G)a;4 =(G—G)by= AG (15b) ; 


I } Rie : 
. % 
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af et aaa ee 3 cane ar = 


Aus (15b) erhalten wir zunachst wegen 6; a;_; — bj; a; = 6; = 1 
G= C,— (B,b;_1 + B,bi), 
Cy = C.+ (Bya;1+ B,ai) 
und nach einfachen Zwischenrechnungen unter Beachtung von (7) 
M,, = C,a.+ C,bz; 
Wi Crack Cy bx + b,41-; By + b,_; By, 


wobei die ibrig bleibenden Konstanten C, und C, aus den Randauflagerbedingungen zu be- 
stimmen sind. 


(16) 


g(6/) a6. 


Mr=( Az + Co bx 
+hj-2 


Me = CO, + Gby 


+/y-7 


Abb. 2. Trager mit gleichen Feldlangen und einziger Feldbelastung in Feld 7. 


a) Freie Auflagerung beider Endstitzen. Aus (16) folgt mit 
a= 0: \M,=0=Ca,+ 636, dh. Cy =0, 
x =n: My’ =0 = Cybn+ Bny1—-iBy + bn: Bo, 
1 
pete Bn. (bn41-1 By + bn: Ba). 
Mit Beniitzung von (7b) und des Vorzeichengesetzes (6 d) wird schlieBlich 
by oe 1 
M, = — — (bn41-i.By + b,_i By) = (—1)*-1-* |b. br (|bn+1—i| By—|bn—i| By) , 


bn 
; 1 
M, =— bn Or aaig berlin b,+1~1) By + (bn—ibz— bn by_:) By] = (17) 
bn—x . <i 1 ‘ ; 
=~ (bi-1 By + 6: By) = (—1)**** [b,_.| bE (|b:-1| B, — |bi| By) . 


b) Starre Einspannung beider Endstitzen. Unter Beniitzung von (7a) folgt aus (16) 
mit den Randbedingungen fir 


“x=0: 2M,+M,=0, C,(2a)-+4a,)+C,(2b)+b,)=0, db. C,=0, 
x=n: ‘M’,,-,+2 M’,=0, C,(an—1+2 an)-+(b,-:+2 6,44—3) By+(ba-1-i +2 b,_:) Bs=0 


oder : 
1 
aaay 3 bn C, —@n41-i B,; ani B= 0; d. h. C, = Soa (G02; B,+4,_: B,) . 
Mit Beniitzung von (7b) und des Vorzeichengesetzes (6d) wird schlieflich 
< 1 
Mas — ae (Gn41—1 By + an—; By) = (—1)'-1~* |a,,| 3 [bal (Jan 41d B, — a, B,) ? 


M,=— sie [(Gn41—-i 42 —3 bn b,41-i) By + (@n—i 4x —3 bn by i) B, |= (18) 


n—xX x —t 1 . 
a a 2 (a;_1 By + a: By) = (—1)1 8 ana] (|a:-1| By — |ai| By) - 

ce) Starre Einspannung der ersten und freie Auflagerung der letzten Stitze. 
Da fir die erste Stiitze dieselbe Randbedingung wie im vorigen Fall b) zu erfillen ist, uber- 
nehmen wir wieder das dortige Ergebnis mit C,=0 und erhalten weiter aus der Randbedingung 
fur die letzte Stiitze 


1 
x=: M’,= 0= Can + Brea: B, + (ane B, : C, = ra (Bn41-i B, ++ bn B,). 
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4. Tabelle der Stiitzmomente fiir freie Auflagerung beider Endstiitzen, geordnet in der Reihenfolge 
der belasteten Felder. 
EATER Gree Ce ee | 


: | M, = (—1)-}-* |b.| 7 (|Pn41—<| By — |bn—s] Ba) 


a(S |b |bi1| B, — [bi| By) 


Ta | 


L 


1 
Lei =—q ([6s-i|B1 — [boil Be) Roi= (bi 1| B1— |bi| Ba) 


: a x 
M,, || La; =< (lbg_i|By—|83_:|Bo) | Raz =< (18:1 Bi — |2:| B: 
| my | ro 2 coacstar-tla-cbBo) | Roam ag (Pia i fl 2) | 
7 | 
R31= —q5 Be . 
1 L % 
Ls,2= 4, (4 B,—B) R3,o= 75 (Pi -4 Bs) 
. 1 . — 


1 1 
Lai = 56 (l0s—i| B1— bail Ba) | Rai = 5g ([be-1|B1— |i] Ba) 


Royse aS , 
4,1 “$56 2 
1 : 1 
L4,9= 5, (15 B,—4 By) Ryo = 54 (Bi—4 Be) 
5 1 1 ; | 
E45= 5g (4 Bi—Ba) Ry,3 = 5g (4 Bi —15 Bo) 
1 
laa 56 


1 1 
L5,i = 399 ([be—i|B1—lbs_ilB2) | Rs,i = Gog (l2e-1] Bi |i] Ba) 
Pauaueeg wan 
eL ~ 209 By 
=2 0 

t +1L5,9 15,2= 909 555 (58 B,—15.B,) Rs,2= 909 555 (B1- 4 By) 
1=3 0 —1L 

eats 15,3 = 909 sa5 (18 B,—4 By) R53 = 909 aaa ( 4B,— 15 B,) 
i=4 0 1L > 

+1L5,4 15,4 = 309 =; (4B, —B,) R5,4 = 909 = (15 B,— 56 Bs) 

=p |} 0 | ths, | +4055 | —15L5,5 | +50L5,5 rake = B, 


1 an a 4 
L6,i= 790 (|bz_ i|By— =|bg_ i\B,)| ee 780 san bE |B ie 13 


: = 209.1] +56R 1 ae oh 
I . —— 
pa SY eae Ce "730-2 
\ i 0 1 
Leo= ie 7g (209B1-56 By) | Res a! (B,-4By) | 
1=3 0 . 
76.3 = 780 ap (888-183) | Rea 780 zap (4B, —15 B,) | 
—___—______i__ 
i=4 0 z arg ; 
Pea = a oleae 42s) Rea = ae (15 B,—56 By) | 
el} poles Pe a eos | 
655 Sane B, —Bp) Res = 7R0 (56 B, —209 B,) | 
i=6 0 gral : : 
L6.6= 799 Ft 
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‘ @a 602 —'a 98) — =x! (%g 602 —'g 082) we sa by 0 
(®@ 9¢ —'g gt) —— =r) (®¢ os, —'a 1162) eS ere OPT TG 
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5. Tabelle der Stiitzmomente fiir starre Einspannung beider Endstiitzen, 
der belasteten Felder. 
re ze a a oy PR 
M,, — eco heater la.| 3 By leat By a la,—i| B,) 


My = (-1)*#!-![a ay rh ae | Ba) | 


geordnet in der Reihenfolge 


1 1 
zi M, Lei = Z5 (a5_i|B1— |ao_i|Ba) Rai = 79 (lial B1— |a;| Be) 
1 1 
+1Ro,1 | Lo = 75 (7 By— 2 Be) Ro,1= 9 (Bi— 2 Be) 
1 1 i 
—1Ry,2 || Le,a= 75 (2 Bi— Bo) Ro,9= yp (2 By— 7 Ba) 


at : DL) 
L3,i= G5 (@a_i|B1—|@5_i|Be) | Ra,i= gg (lae_1! B= ail Ba) 


1 
Ls= ye (26 By — 7 By) 


1 Ra) air 
14i= zs ({4s_s|81—|24—1|Bz) Ryi= 168 (jaz_4|B1— la; |B.) 


—2B,) 


y; al 
| +1Ra | Pa igs (97 B1— 26 Bo) Ry = asi 


iss 


+ 2Ly 2 +2Ry,9 ; —1Ry,o || La,2=Fe0 + (26 By— 7B) Rye = 


168 a @ 


By— 7 By) 


Pas = (4 Ay 4 By) R43 = 


Alar oe ee 7) 71B,—26B 
—2L 433) +1L4,3 Sane ia a a 2) 


‘| Lya= as —— (2 By lle Rya= 168 — (26 ey 97 By) 


j 


* (362 ‘ee 97 ay 


Bob i i —1R5,1 ‘eal fone ee 


Wibcae ean (27 B,— 26 B) — (2 Bie gon 


ca me 


~~ (26 B,- ~ 26 Ba) 


Lag ee zi By) 


a 


=} - 


~ ea 2G By 22) 


a7 OT Bs 362 Ba) se 


~ 627 


A eM eet ra ia 
CaN MSIE f cn 1 


zug (1951 By —862.5,) 


B,- 
| Fo1= re a0 


2B) 


Roa = “2340 


% 


1 
=— (2 
Rea aa tisk 


~ cea 
UO) : 


raat By 


“Be 


a5 82 By = 1351 Bs) 


~ 9340‘ 


Ar 


Ingenieur- 


hen Stiitzen.. 


t 


i 
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6. Tabelle der Stiitzmomente fiir starre Einspannung der ersten Stiitze und freie Auflagerung 
\ der letzten Stiitze. 


N 


a Mes (-11 Ia a (Boe B, ~ [bni| Bs) 


| | eg be ae ee = (lal B,-— |a;| B,) 


L di 
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+ i i 
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ky ; 1 z 
: 0 Lo,9= 7 B 


[oe I 
36 (15 By — 4 By) Rs1 = 56 (Bi — 2 Be) 
a 
— Bs) R39 = "56 (2 Bi— 7 Bs) 
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7 tea Metre hor eb Sas eoer nO Laa= 35 B, Be) ae roe vt 
oer : ae 
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SchlieBlich kommt 
ay, 
M, =— “as. (6,472 =% B, + bps i B,):= —- (—1)i- I-a/g,|—* lan) r (busi il B,— |b n-i4| B,) ’ 


, 1 
M, are ae as | (Bit1i@x—anb bxi1—:) B, ai (b n—i Ax —An by i) B,| — (19) 


Baas : : 
a (a;_1 B, + a; B,) = (—1)*4t|bn—-| Jan) ( 


In den folgenden Tabellen 4 bis 6 sind die nach 5a bis c ermittelten Stiitzmomente fiir homo- 
gene Trager bis zu 10 Feldern in iibersichtlichen Zahlentafeln zusammengestellt. Sie enthalten, 
fortlautend geordnet nach Einfeldbelastungen vom ersten bis letzten Feld, die Zahlenwerte fiir 
die Stitzmomente in Abhangigkeit von neuen, durch zweckmabige Zusammenfassung gebildeten 
Belastungsgliedern Ln, und Rn. 


a;_| B, Fane |ai] B,) . 


d) Bestimmung der Auflagerkrafte. Beriicksichtigen wir, daB nur das i-te Feld be- 
lastet ist, so kénnen wir mit dieser Einschrankung die Ergebnisse aus (12) ce und 
erhalten folgendes Schema: 


Ao Ax A;_y A; A’, A, 
freie Auf- M, : My 
lagerung | 1 M M; 1 5 Me i ee ee | 20 
67 | 6+ f d-9raeae | 65+ 1 f pqeas Ce 
-starre Ein- 3 My Rare G , 0 M, 
spannung , ts tye 1a ; at ips 


- Hierin aad wieder die Stiitzmomente fir den entsprechenden Lagerungsfall einzusetzen, wobei 


die Werte aus den Tabellen 4, 5 oder 6 fiir die betreffende Felderzahl aus der Spalte fir das 
belastete Feld entnommen werden kénnen. < 


6. 7 wenieccane Der von Clebsch angegebene Ansatz zur Lésung der Differenzen- 
gleichung fir die Sttitzmomente des homogenen durchlaufenden Tragers wird durch Einfihrung 


konjugierter Zahlen a, und b, in eine fiir die Rechnung bequemere Form M. =7 B+ Cy ax+ Cy bx 


gebracht. Nach dieser Darstellung werden die drei Lagerungsfalle des homogenen durchlaufenden 
Tragers behandelt und die Stiitzmomente in Abhangigkeit der Belastungsglieder in Form von 


_ Zahlenpyramiden itbersichtlich zusammengestellt. 


SchlieBlich werden nach derselben Methode die drei Lagerungsfalle eines Tragers mit sliehen 


; Stiitzweiten und einer einzigen Feldbelastung behandelt. Mit dieser Belastungsart als Unter- 


belastungsfall kann ein Trager mit gleichen Stiitzweiten und beliebigen Feldbelastungen 
nach dem Superpositionsgesetz durch Uberlagerung aller Unterbelastungsfalle geschlossen be- 


-handelt werden. 


‘tr * 
a: 


: (Eingegangen am 16. September 1948.) _ 
Saas “des Verfassers: Dr.-Ing. Johann Fadle, (17 ee tee One: Geors Biel -Strabe 17. 
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Der Einflu8 der Dampfung bei nichtlinearen Schwingungen. 
I. Gedimpfte Pendelschwingungen. 
Von F. K. Rubbert. 


endlicher Amplitude ist bisher nicht befriedigend beantwortet worden. Und zwar handelt es 
sich dabei um die der Geschwindigkeit und dem Quadrat der Geschwindigkeit proportionale 
Dampfung und die Festreibung. Die analytischen Schwierigkeiten dieser nichtlinearen Differen- 
tialgleichungen sind derart, daB selbst rohe Naherungen noch nicht angegeben werden konnten. 


Dampfung numerische Methoden angewandt. Bemerkenswert ist die Arbeit von Zech, der zu 
Aussagen uber das Abklingen nichtlinearer Schwingungen mit Hilfe der ,, UmriBlinien“ gelangte. 

Wahrend bei der: Differentialgleichung fir quadratische Dampfung durch Anwendung der 
Identitat g = 4 d(@2)/dy die Trennung der Veranderlichen ausfithrbar und erst die daraus folgende 
Differentialgleichung erster Ordnung nicht in geschlossener Form lésbar ist, besitzt die Pendel- 
gleichung fir lineare Dampfung iiberhaupt keinen Angriffspunkt. Diese Schwierigkeiten lassen 
sich umgehen, indem man die direkte Integration anwendet. Das Ergebnis erscheint hierbei in 
Form einer Reihe, die nach Potenzen eines Parameters fortschreitet und deren Koeffizienten 


ist von Poincaré* untersucht worden. , ea 
Die Differentialgleichung des Pendels mit endlicher Amplitude und linearer Dampfung ist 
bereits von Moigno®, Broch® und Horn’ auf diese Weise behandelt worden. Aber alle diese Lé- 
sungen kénnen nicht richtig sein, da sie fir verschwindende Dampfung ungiltig werden. Damit 
ist aber auch schon der Punkt beriihrt, dessen Beachtung die wahre Formel liefert. 
In einer zweiten Mitteilung soll die Wirkung der Dampfung bei unsymmetrischen Schwin- 
_ gungen behandelt werden. 


2. Das Pendel mit quadratischer Dimpfung. Auf ein ebenes Pendel mit dem Tragheits- 
moment @ und dem Schweremoment S wirke eine dem Quadrat der Geschwindigkeit pro- 
portionale Dampfung; seine Bewegungsgleichung lautet, wenn den gegen die Gleichgewichts- 
lage gemessenen Auslenkungswinkel bedeutet, i 


P —OG-+ (sgn G) Ry? + Ssing=0. ty ah) 
Hieraus erhalt man mit den Anfangsbedingungen _ 
p(0)=a, — o(0)=0. : (2) 


_ fir die erste Halbschwingung die Differentialgleichung ; 
Pe ane fe SR TT SID eo MES iss, oe (3) 
hee f, mit ten: ; “ 
de 8) 1 
; ‘ re ' 2 Sf. 23 
re pay ss fe: Bis Sas : i 
Wir fihren zunachst mit einem noch offen gelassenen Parameter x eine neue unabhangige 
_ Verdnderliche ein: | Ga . Cen ge Hea ke fats 


One i ry! Via ; r ; beh : ‘ (5) 
und setzen — : ey £ nog: Se 
A : qq = aa : (6) 


1 F, A, Willers, Z. Instrumentenkunde 53 (1933) S. 504. 
2 A. Spaeth, Dissertation, Gottingen 1940. i 
8 Th, Zech, Ing.-Arch. 13 (1942) S.21. : 
S. 157 #f. : bas 
5 M. Moigno, Calcul intégral (1844), -S; 584485), tS | 
° 0. J. Broch, Lehrbuch der Mechanik (1854), S. 207 ff. . ‘i 
” J. Horn, Z. Math. Phys. 47 (1902), $.400, §§5,6. % 


ra 


nacheinander aus linearen Differentialgleichungen berechnet werden. Leichte Umformung ergibt — 
schlieBlich die Fourierentwicklung der gesuchten Lésung. Die Konvergenz dieses Verfahrens_ 


Mae “oR aca Ong 
Audie erg gh=6- : ae eae 


« H. Poincaré, Méthodes nouvelles, Bd. I (1892) S. 58 ff.; auch E, Picard, Traité d’analyse, Bd. III (1896) _ 
Wl j , ie S ad t ai ee ne Po apes 


1. Einleitung. Die wichtige Frage nach dem EinfluB der Dampfung bei Pendelschwingungen | 


Willers! hat deshalb bei der quadratischen Dampfung graphische und Spaeth? bei der linearen 


. ne ; \ 
a fanart Aho as wrcdistarn ROT Ne ON 
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Die Differentialgleichung (3) geht dann tber in 
P+ (+x) p= agp + (1+) (3 98- a p+ +), (7) 


wobei die Striche nun Ableitungen nach t bedeuten. 
Mit den Ansatzen 


p= «4D, (t) + o> DB; (t) + a5 (zt) +--+, (8) 
%= a eo+ atey+ -:- (9) 
ergeben sich fiir die ®;(t) die linearen Differentialgleichungen 
O/ +0, =0, - (10) 
@/ +O, =— 0,0, + ¢02+1 63 
3 3 eer g Py 6k - (1d) 
” Hip eu ee | 1 
Di +, = — ¢Ds— eqD, + 2q HO, + = Oi O, + 5c, OB} — J GF (12) 
_usw. mit den Anfangsbedingungen 
@,(0)=1, @,(0)=0, (13) 
@,(0)=0, @,(0)=0, (14) 
@,(0)=0, @;(0)=0 (15) 
usw. 
Hieraus erhalten wir zunachst 
@, =cos T ae (16) 
und damit fir ®, die Differentialgleichung 
@, +O, = sat (s ~ ¢,) cos t~ 5 q cos 24 + 5. cos 37. ~ (17) 
Damit auf der rechten Seite das Glied mit cos t verschwindet, muB sein 
aa Se) 
Als allgemeine Lésung der COAG ara (17) ergibt sich dann 
1 1 1 
®,= 54+ (a5- 3 yao ets qcos 27 — 755 cos 37, (19) 


und hiermit bekommen wir fir ®, die Ay ey 


O, + O, = (oI 9") ee ; date q *— ¢,| cos T + (Gata) cos 27 + 


(20) 
1 1 I 3 
+ (563 I-57) cos Brt oh qcos 4r— a560 ©°S Dit 
‘Die Bedingung fir eine periodische Lésung lautet jetzt 
List : 3 
a Lae Ne Gh 
und damit die allgemeine Lésung 
eeut Bou 17 27 ei the 
O° (int 3) + Caras matt ms')om* oat 5 q°) cos 21 a) 


1 
aor sé Ts: a 1 *) cos 3T — a q cos 47 + iit cos ST. 


‘Durch Einsetzen der Wasdracks (16), (19) und (22) in (8) erhalten wir daher in dritter N aherung 
far die Fourierentwicklung der gesuchten Lésung 


yg =A,+ Ay 00 ¢ bly one 247 4 Ay cos 37+ Ajcos'dr + Ay os ae (23) 


bis den eS 


1 Role ee ae ey 
27 
tone Fah sgt ht 1024 ath + siaay oe (25) 


Ag = ye 2h sna She 7 ad he , (26) 


eee Ss : ab tae Lea Sie. || 
: 4 z ae || 
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a er On SS SE SS ET ETE SESS (5 > TEES SSDS ih 
A Fo + oath? +o ath ag 3 ez || 
3 192 ~ 6 155 3072 || 
me (28) | 
se 4 
Ay= — 3975 0 > 
= I 5 ; (29) } 
As = xyag0 *" ° | 


Der erste Umkehrpunkt wird erreicht fir t=2. Fiir die Dauer der ersten Halbschwingung 
folgt daher aus (5) nach Einsetzen von (18) und (21) 


5 i 2 i Ties (30) 
T= = 14 7a? (13M) — yah tpg 

oder 
cd 1 2 32 1 3 yal 30a). 
T= 5 [1+ 7g 07 (1+) — ao h + ao7g %" |: . Stee | 


Die Amplitude im ersten Umkehrpunkt ist 
. Pr = Ay— Ay + Ap— Ag + Ag— As- 
: Die Ausrechnung ergibt : 
Baer it ps|=e—S ath +o eM) eg Ce) (31) 


Fur h= 0 reduziert sich Gleichung (3) auf die der ungedampften Pend diesen oe endlicher | 
Amplitude | 
g+n'sing=0. (32), 
Ihre strenge Lésung bei den Anfangsbedingungen (2) ist bekannt und l4Bt sich durch Anwendung |f 
elliptischer Funktionen in geschlossener Form angeben: 


: cn nt | 
gy = 2 arcsin G ae =) : (33) 
“mit 
k= sin. r (34) 


Die ersten Glieder ihrer Fourierentwicklung erhalten wir, wenn wir in (23) bis (29) h=0 setzen; |} 
es folgt 


Geer 15 Sed 2 oe 
y= (eta 84 ag a®) cos + (a55 a * 3073 % a8) cos 31 + xara 0 cos 5T (35) | 


es | — — ——. (36) 
eS, Jiteetuet 

‘Auch die Dauer einer Halbschwingung ]a8t sich in einfacher Form darstellen, und zwar 
durch Anwendung des vollstandigen elliptischen Normalintegrals erster Gattung: 


pe es i). ) | (37) 


pte 
ae he ea (Go 


Die Entwicklung der hypergeometrischen Reihe ergibt in dritter Naherung 


/ 


ee aoe oa a oie? ek at) 


(38) 
was genau mit Formel (30a) fir h=0 uibereinstimmt. 7 
Ein weiterer Sonderfall, der in unserer Lésung enthalten ist, ist der der Kleinen Schwingungen 7 


mit pee deaciether Dampfung. Ihre Diftopmnaletchaek lautet § 
(a aa 5 h (sgn @) G2 + n2g SOE RAC EE See 


Vernachlissigen wir in den durch (24) bis (29) aniceienees ibefiecascn alle Glieder, die von ; 
der endlichen Amplitude herrihren, so folgt bei den panfabesbedinBahuen: (2) fir die ersten Glieder 
ihrer Lésung als Fourierreihe _ ; 


=( a ya es ath) +(1-+e bee iB oon et 4 
© 6 8 ne (40) ; 
+ (jgab— 7p athe) cos 27+ a a2h? cos cos 3T pes 


18 


ee DAVIE“ < : nt } 


: ee t= (41) 
foie - \2 ee as =o h? | 

_ Far die Spee der ersten Halbschwingung ne wir = | i, 
oe eae RS a* h?) fs (42) eae 
_ und fir die apres im ersten pags nel 
‘4 - Pk : 
3 Sane eae eal = & tes ~Sah+ > 5 ath?) | | (48) 7 anes 
__ Die Formeln (40) bis (43) sind bereits auf eine etwas andere Art von Burkhardt! abgeleitet worden. *)\) = so y - 
3 EES pepe eed mit Reibung. Seine Bewegungsgleichung lautet | eye ees: ey 
7 a es 9 +n? TS Ae (sgn =, 0, Ete 3 (1) é 
: und. ie Anfangsbedingungen seien se: fis whee : Jew 
4 e pO)=2, p0)=0. : Ql ee 
3 Wir fahren wieder eine neue unabhingige Veranderliche ein” durch. A Ze gee 

Py pete = ~ nt ‘ arr 
ao ; Regs Set 3);! 
ae 4 : c Ky Vile Ee 23 ; oa ; : ; { ( ) ; a 
2 Die ‘Differentialgleichung fur dic ersté’ Halbschwingung ist dann oT a= cores : 
c ” 1 
(gp F140) paacllt%)+(1+%) (se —aort-) At) ie 
a ; Fé. ; 8 opus ; : p oN j . a. As 
Mit a itis p= 46,(t ( 1) + «2 @. (ut) + 25D tye eee Bes (6) 
DN OO daa gat yd NE be oe 
4 bekommen v wir far die @;(z) die linearen Dilfercnsalgtictinngen’ tS a cae a “ 
. he SC See ee Me pay. tec 
we (WO =a Dit 4 1 gs, me SWEET re se re 
eee eu 4 ionita Fa ere 2 A 
OO =o %— Ears e Sa: O45 cy cD} + 52, Pon (10) 
{ 2 Borges : 


s —— 


eo Y ec US ay 
BPO) BL SS inh ee el CR) ty 
Xs ee me Oy areas. | 5. te ak some (13) 


re a: Sel) po y saath alts ayy £ 
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e 


Durch Einsetzen von (14) und (17) in den Ansatz (6) erhalten wir fir die zweite Naherung 
der gesuchten Lésung als Fourierreihe 


: | | ee 2 s)? 4 Ba s)? + og ( SNe (as) |) 


| 3 s 
3 T— t ry 


Der erste Umkehrpunkt wird erreicht fiir t=. Die Dauer der ersten ifalbaohwiewene ist | 
daher nach (3) ; 


; 1 i 
se T= 2/14 ; (% —s)?+ 5 Se (19) 


oder 


aye : gies | i 19a 
T, [1+ 36 (e s) Bey a : | ( ) 
Fur die Amplitude im ersten Umkehrpunkt ergibt sich 


lm] = 0-28 — Fs (a= — 5s (20) 
Fur die praktische Anwendung wird die zweite Naherung stets ausreichen. Es macht aber 
durchaus keine Schwierigkeit, auch die weiteren Glieder zu berechnen; nur die Rechenarbeit — 
nimmt mit jedem weiteren Glied erheblich zu. 
Auch die Differentialgleichung eines Pendels mit Reibung und quadvatiocher Dimpfung 
laBt sich auf diese Weise mit jeder gewiinschten Genauigkeit integrieren. 


4. Das Pendel mit linearer Dampfung. Die Bewegungsgleichung lautet jetzt — 
o+2kp+nising=0. (i) 


: : Die Anfangsbedingungen seien wieder : i 
eG : g(0)=a, p(0)=0. . ree eae (2) 
3 LU ee ee a (3), (4) 


2 
i 


und bei Beschrinkung auf die zweite Naherung geht (1) tiber in 


. d? d 1 ae 
Sha ih 7 qe +26 Ft ae ae eens ihe SSN | 


» Wir en an einer Lésung dieser Differentialgleichung, die fir kleine « in das bekannte a 
Integral der linearen Gleichung tthergeht und sich far 6=0 auf die Formel (35) von Ziff. 2 far 
die ungedampfte Schwingung endlicher Amplitude reduziert. Das kénnen wir aber nur er-— 
5 Bee - reichen, indem wir eine neue unabhangige Veranderliche @ einfihren durch 4 
Wreuaey on Saat ae T=O+feqX(l—e-29%), . eave ye | 
4 wo die Roun c As der Integration in der gewiunschten Weise bestimmt wird. 
a Weees., c a ane id site 


: ] dp. ~ d2@ d2¢ } rte for 
Oy MeO Ae ep. een bp rae aes 
Gir Be ee dh (sae ir YUN rhein we orn ea ee 
| RRC SE EN a OR Neer ieee 
t in zweiter Naherung ey . a se POC 
Wo aw he 2 —260 Necriae ee 
a ee (l 2a2%cde ).: : ies Eafe 
ee dp d Q ‘ it a Dae 
a foe oe Anith Soe FE (datedt 200), 
£ 1 ‘ t 4 
3 
=<9 % ’ a 
oe ae on 
ip se yea i ; 
ui ¥ tc docs % 
is ¥ . ? 
; : io ci dN a 


eee a a, a en ae 


- 
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Zu ihrer Integration machen wir den Ansatz 


p= 4D, (8) + a Bs(9) (8) 
und erhalten fiir ®, und ®, die linearen Differentialgleichungen 
@/+260,4+0,=0, (9) 
O, +269, +B, =1O8 + 4050-200 0% (10) 
mit dee Anfangsbedingungen | 
P, (0) = ®,(0) = 0, (11) 
@,(0)=0, O,(0)' > 0". | (12) 
Hieraus folgt fir ®, 
, = 0-9” (cos yO + <sin v9) (13) 
“mit 
y=) 1-8. (14) 


Einsetzen von (13) in (10) ergibt fir ®, die neue Differentialgleichung 


e399 (Di + 260@;4+@,) = (2-400) cos ee mers 2 ge) tee 


15) 

1—40? —4 se S ( 
yo ue sin3y0. 

Damit unsere Lésung auch fiir 6=0 giiltig bleibt, muB der Py enent von cos y? verschwinden, 


d. h., eS mu sein 
1 
To GR ean (16) 


Gleichung (15) reduziert sich nun auf 


Ye 5MD oD, =e 800(, sin yP + pcos 3 + ae sin 3). (17) 


Ihr allgemeines Integral lautet bei Beachtung der Aetopaingne (12) 
11+5*—126' 6 4T—74.62+2468 


—e—60)| i e368 : 
ie Masel tenet yo] + ioe (cos 70-+| (is) 
ORS, —360] 1-108 4126 6 9—320+42464 | | 
+<;sin 78) e532 prosody y 9672 (4—362) sin fol : 
-Mit (13) und (18) bekommen wir als Lésung unserer Differentialgleichung 
= (Ay, e999 + Aj, e—? 9”) cos yO + Azz e329? cos 3 ae (1s) 
+ (By, e— 99+ Bys e~99®) sin y+ Bg, e~3°% sin 3 yO, 
wo . 
114 6?—12 64 6, ATO TAB L240 
ie any 48 y?(4—3 6) a”, Bia 9673(4—-3 8)’ 
1 6 
Ax; = 1672 a 3 oo ? By — 167 a > 
, 11024126! 4 i 96-3224 08g 
ce he Te 308) Osh ee eg ys (4 388) 
ist und oO durhades Gleichung ‘ 
t= O+ 5yr = (1—e-? 49) : (6a) 


bestimmt wird. 
Der erste Umkehrpunkt werde erreicht fir 0=0,, und zwar ist dann 


Oy) = 2 H(O1) +2 BO.) = 0. 


Setzen wir nun 


4h: 4 - Sat en = 2 wh at = - a 5 Sud | 
nee x oy ? cables pak. " . Tt 
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so folgt, da ¢ als sehr klein angesehen werden kann, 


B(O,) =e! (FZ), a2 H,(O,) = 2H, (F), 


also 
/ 4 
ae Siac ee 
wr[ 
Pika) 
Die Ausrechnung der hier auftretenden Ableitungen von ®, und ®, fihrt schlieBlich zu 
hg 28.6248 20 | 20 
O1 = | + 0" a5 aa — 3%) (1 i): 20) 
und damit bestimmt sich nach (6a) die eee: der ersten Halbschwingung durch die Formel |} 
_ ; 
4 — 
re T, = [Oy + 45 Rie (Ie?) |. (21) 
Die Amplitude im ersten Umkehrpunkt ist mit ausreichender Genauigkeit 
i] 
= (0,) = be 3p, (2). 22 
. ~(0;) aD, (=) + « ®; (=) (22) 
‘ Entsprechend besteht die Beziehung , 
es a nm . | 
: Pe+i = VePy =) + VPs en ; (23) | 
wenn g, die apnea im es ten Umkehrpunkt bezeichnet. Hieraus folgt das allgemeine Ampli- 
tudengesetz | | 
af Peale. » 11+6*—12 54 ene . 
Pk ae BS aa 48 2 (4—3 62) (1 e 7) |. (24) q 


Nachdem g, bekannt ist, ergibt sich aus (20) und (21), wenn auf den réchten Seiten « mit ¢, 
vertauscht wird, 02-0, bzw. T,. Auf diese Weise kénnen nacheinander die q, und Ty, be- 
rechnet werden. 


Zum SchluB sei gezeigt, daB diese Ergebnisse auch die bekannten Grenzfalle enthalten. Fir 
6+ 0 geht (20) itber in O, =, und damit ergibt (21) 


, T,=T=2(1+ 7502), 


also die Schwingungsdauer fiir ein ungedimpftes Pendel endlicher Amplitude in zweiter Na- 
herung. Auch der fiir  gefundene Ausdruck (19) erhalt dann die gewiinschte Gestalt — 


\ 


Q= («+ jogo? )e0s 9 — 355 0° cos 30 


mit 
2 nt 
\ oy 0 = eager 4 
; ‘ a2 
LF 5h 


ae a <1 kann age zweite Term in (20) und (21) verneébliceigy werden; die Formel fir die 2 
Schwingungsdauen.nimmt daher wegen (4) die einfache Gestalt an 


re sos set ; 
Entsprechend reduzieren sich fiir kleine « das Amplitudengesetz (24) ) und die Bewegungs- 
sploiching (19) auf die bekannten Formeln aus der Theorie der kleinen gedimpften Schwingungen. 


y 
4 
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Die elastisch-plastische Biegung des krummen Stabes 
unter Beriicksichtigung der Materialverfestigung. 


Von W. Swida. 


1. Einleitung. In einer unserer friiheren Arbeiten! haben wir die elastisch-plastische Biegung 
eines krummen Stabes unter Zugrundelegung eines idealisierten Spannungs- Dehnungs-Dia- 
grammes mit einem langen horizontalen FlieBteil betrachtet (Abb. 1). 

___ Bekanntlich hat aber die Spannungs-Dehnungs-Linie fiir viele Werkstoffe keinen horizontalen 

FlieBteil. Das bezieht sich insbesondere auf einige Stahlsorten, Bronze, Aluminiumlegierungen. 
‘Die drei in der Abb. 2 angegebenen Kurven stellen die normalen Spannungs-Dehnungs-Linien 
fir Manganstahl (a), Nickelstahl (c) und Bronze (b) dar. 


Abb. 1. , Abb. 2. 


Fur die Untersuchung der Spannungsverteilung bei der Biegung der Stabe aus solchen Werk- 
stoffen kann man zwecks Vereinfachung der Berechnung ein idealisiertes Spannungs-Dehnungs- 
-Diagramm (Abb. 3) annehmen. aoe : 
Es wird angenommen, daB Om=Op und mn 
parallel pq ist. Dann kénnen die Zusammen- 
hange zwischen Spannungen und Dehnungen bei 
Zug und Druck als identisch betrachtet werden. 
Man ersieht aus Abb. 3, daB tg «=E und tg a, 
=F, ist. Dabei ist EF, der Elastizitatsmodul far 
_diejenigen Spannungen, welche die Spannung os 
iibersteigen. Diese entspricht den Knickpunkten 
‘der Spannungs-Dehnungs-Linie?. Auf diese Weise 
wird als Ersatzstoff ein idealisierter Stoff _ 
“mit einer Kennlinie nach Abb. 3 angenommen. 
Da die Entlastungslinie rs (Abb. 2und 3) im- 
‘mer als parallel der urspriinglichen Belastungslinie 
“om angenommen werden kann, ist-der Abschnitt 
-os Jie bleibende (plastische) Formanderung, und _ 
der ganze Vorgang kann als elastisch-plastisch — 
-bezeichnet werden. are 
Die Bericksichtigung der Spannungen, die o; tberschreiten, fihrt, wie gezeigt wird, zur 
richtigen Beurteilung der Tragfahigkeit des Stabes und kann Materialeinsparungen bei der Be- 


messung bringen.. 


- Die Gleichungen der Geraden pm, mn und qp sind 


Abb. 3. 


Z 


Ss o= Ee, \ o=60+ Eye, o=—Oot+ Eye. ‘ (1) 
‘TW. Swida, Ing.-Arch. 16 (1948), S. 364. 
2 Die Untersuchung eines geraden Stabes 
linie wurde von A. Néddai, ‘Der bildsame Zustand der Werkstoffe, 


e : . rf a 


S. 121, Berlin 1927, durchgefiihrt. 


unter Zugrundelegung einer Ahnlichen Spannungs-Dehnungs- 
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Wir nehmen ferner an, daB die Piepuny even e in der Krimmungsebene wirken, daB eine | 

der Symmetrieachse eines beliebigen Querschnittes ebenfalls in dieser Ebene liegt und betrachten | 

den Fall der reinen Biegung. Dabei wird vorausgesetzt, da man die seitliche Pressung zwischen | 
den Langsfasern auBer acht lassen kann und daB die Querschnitte nach der Biegung eben bleiben. |}} 

Auf Grund dieser letzten Voraussetzungen ergibt sich fir die relative Dehnung ¢ in einer be- || 

} liebigen Faser (Abb. 4) derselbe Ausdruck wie bei der elastischen Biegung: 
a ody y : Hl} 
a C= dg r+y ’ (2) : 
worin r der Kriimmungsradius der neutralen Flache ist. 
Die Spannungen bei der konkaven Oberflache sind in der Regel gréfer als bei der konvexen || 
Oberflache. Deshalb werden zwei Spannungszustinde unterschieden: 
1. Die plastische Zone (Bereich der Spannungen, die o; tiberschreiten) ents nur auf einer | 
Seite der Nullinie, und zwar an der konkaven Oberflache. . | 
2. Die plastischen Zonen treten beiderseits der Nullinie auf. 


2. Erster Spannungszustand. Auf Grund von (1) und (2) ergeben sich folgende Werte far | 


die Spannungen o bzw. oz in der elastischen bzw. in der plastischen Zone: 
é 


Eddp y 
a 9 ‘ 3 
dp or+y (3) 
E,édg —y¥ \ 
COur= a 4 
oe dp r+y (4) 


NY” 


Abb. 4. etic Abb. 5. 


Ist x= ie (Abb. 5) so wird jede von diesen Spannungender absoluten GréBen nach gleich | 
os sein. Also ist 
édp et os(r—y) _ (6.=6a) 7=91) 


EN ag, ih i dae” Ey E : ; Cy mas (5) 4 
Aus (5) folgt ij es AIO ee 
Senge AG Asse oie O;—Oa Rees pe Pees 7 fete : eee | 
ES net ; Shasteen ACE oe 5 ; eo: r | Tag (6) | 
Man kann auch. roa letzte Ergebnis unmittelbar aus Abb. 3 bthalten: ON Bes 4 


f - Bildet man die Summe der Projektionen samtlicher auf das Element des gekrimmten Stabes a) 
ig _ wirkenden Krafte auf die pee der x-Achse (Abb. 5), so verhalt man 


t 
ee 


ae ae ray eee Hoe Rg rey t)ar=. Be les oss: 
Daraus folet sae high Bee ; Cela ; 
aos 8 ody By On Baye a ; 
| dg pyar oe pyar (7) 
\ ‘ 7 ry [Ss 3 
Ga Gy 


Dabei i ist F, der. Flicheninhalt des Querschnittes as clastischen Bereiches; F ist derselbe fiir ‘ 
den plastischen Bereich. : : eet 


~ if a 


{ 


- XVII. Band 1949 


Swida: Die elastisch-plastische Biegung des krummen Stabes. 345 


‘ 


Das resultierende Kraftepaar der tiber den Querschnitt verteilten Spannungen mu mit 
dem auferen Kraftepaar M im Gleichgewicht stehen. Daher ist 


Edd 2dF E, 6d 2 Eéd 
M =f Ne eee) ab Dey tema aed us Y | 
eae: ea r ( ues roy) oF ine 


° 


dp 
(F,,) (F.) 
E, ody rodg | yd F ydF\-: 
s(Oigt foot (aes J). 


Die ersten zwei der erhalténen Integrale stellen die statischen Momente der Querschnitts- 
flachen des elastischen und plastischen Bereiches in bezug auf die Nullinie dar. Man bezeichnet 
sie mit S, und S,. Mit Hilfe von (7) findet man aus der letzten Gleichung, daB 


ddp  M+oq(Sy+Fu1) | : 
dg ES.+E, S, +48) 


ist. en 
Man kann nun die Gleichungen (7) und (8) auf Grund von (15) folgendermafen umformen: 


Og F,, os Os (r—¥3) (9) 
eee fa ey ve 
eee Ey, 
M-+oa(Su+F ur) 65 (r—41) (10) 
ES, +E S., Bye 


» Aus den Gleichungen (9) und (10) kann man die Werte von r und y, ermitteln und auf diese 
Weise die Lage der neutralen Achse und der Grenzen zwischen der elastischen und der plastischen 
“Zone bestimmen. ~ a 
Setzt man ddqy/dy aus (8) in (3) und (4) ein, so ergeben sich folgende endgiltige Formeln 
fir die Spannungen im ersten Spannungszustand: 2 a 
| M+oa(SutFur) 


| (eee Gok BRE Somat : oe 
Ss ee eis ; & [M-+o.(Su+Fy r)} 3 ; : 

ae : Ou = Oa + S.42S, rhy ’ (12) 
i . ‘s E, : ; 

 worin c= Se ist. 


Setzt man in (10) S,=0, F,=0 und y, =h, ein, so ergibt sich der Wert des Biegungsmomentes, 

_-der dem Entstehen der Spannungen o; in den auBersten unteren Fasern entspricht: 
RR oe ercomzags See, 
ie : | hs 
wobei S das statische Moment des ganzen Querschnittes in bezug auf die Nullinie ist. In diesem 
Fall ist r nach den bekannten Formeln fir die elastische Biegung zu bestimmen. 
7 ‘Setzt man in (11) y =h, und o =o, ein, und ermittelt man M, so erhalt man folgenden Wert 
- des Biegungsmomentes, der das Entstehen der Spannungen o, in den auBersten oberen Fasern 

bedingt: Sets 


(13) 


| Mm ou (Se+ £ 8.) E04 (Sut Fur). (14) 


¢ Gleichungen (14), (9) und (10) sind nach M, r und y, aufzuldsen. 


‘, 


Laan Eire Ee) 405) 
3 dp r+y ; 
“An der unteren und oberen Grenze awischen der elastischen und den plastischen Zonen 


as y be ve ie vb H \ | Oy = Oa 


- miissen die Spannungen der absoluten GréBe nach gleich sein. Folglich ist 
ee a i sf aa ah es carat A i ' 
Bin SME No oa it eee | (16) 
pba en as Dew R eat 6 ek ee 
Wiehe oa ee ee Beis Nis 2s CNG’) 
Wentie ead Oe 
re Ai s 22 
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Bildet man die Summe der Projektionen samtlicher auf das Element des gekrimmten Stabes | / 
(Abb. 6) wirkender Krafte auf die Richtung der x-Achse, so ergibt sich | 


Eddp ydF 1 Evodp oy Vp eee Pe Ee ae 
petite Pfu “Jl 


dp rty | dp r+y/ dg ry 
(F) () | a 
wobei Fy die Querschnittsflache der oberen plastischen Zone ist. Daraus folgt 
(<8 dp Oa (Fu—F) (17) 
dg E f: yale / ydF 4 yd F 7 
rey! aa hit r+y ) 
(F,) (F.) (F) 


Die zweite Gch oh Beit ROvIMRL NS erhalt die Form 


nia ydF Hodge ov yy EOdp. iy) Yah 
eet pe r+y ae ort eee dp oo ae: (aay + dp Zor, 


es B.tde oo - nl . 4|| 
= aus fy dF+ (=<? a.) fy dF +( +a ) [yar ; (18) | 
() ls (ies 


rédy Wed Bes) Y yd F ia 
ee E : 
ae ay a cs (J ans ; as )| 


Die -orsteniiieet dee _erhaltenen ee ales die statischen Nisaieute der slasteHen und | \| 
plastischen Zonen in bezug auf die neutrale Achse dar. Man bezeichnet sie mit S., Su ig Sor 
Auf Grund von (17) ergibt sich aus der letzten Gleichung 


Sut tea ae bd | M+ oa[(Fu—Fy) + Su— Sol an 
, a dg ES.+E,(Su+Sp) : pees | 
_ Die Gleichungen (17) und (19) lassen sich mit Hilfe ¢on (5) folgendermafen darstellen: __ 
: j Fa(Fu—Fo) ; be Gs (r—Y3) ny ) (20) 
dF dF EEG RSE: te 
x E ae i y ie ef é 2 8 
te ( area tala 3 
| oe ae F,,) (Fy) 3 * 
| Mica [(Fu=Fo)t+Su—So) 379) en] 
ES.+E, (Su+5p) - So yy a 


_ Aus den Gleichungen (20), (21) und (16’) kénnen y,, y, und r ermittelt werden. Setzt man | 
in (3), (4) und (15) den Wert von ddq/dg aus (19) ein, so erhalt man aie endgilltigen pC | 
fir die Spannungen im zweiten Spannungszustand : 


Vipeseee M+6.[((F.—Fy))t+Su—Sy)]) sy < ) r Re 
oS sarge S48 S45) Tey ne | 

eh oe et Maa Fle SeeSil fr ys ger | 
eo Meee Oona Ge v SePE(SAE Sp) ere e und 24) 


_ Nimmt man in (9), (10), (20) und (21) an, dab E, 0. 
ist, und ersetzt man 6a durch o;, so ergeben sich die | 
-frither! gefundenen Zusammenhange fir die Biegung — 
des gekriimmten Stabes aus einem. idealplastischen 


Werkstoff: Ue paidess | eae | 
ir 
of r+y ea hae i iL, ae OR tag 
2 yeas 2 ptr enema ee ee | 
: ; oN, cae — =) ere eer, | 
nyt yd F ae Hie , ie . ‘ a 
hege es =, ey Fe A bat (20 ] 
(F.) lacie Leer * “ ety wap chic ey 24 
: Pe + Fu (le—1) 4 Fy het) © (21 
¥ Tess t 721 oe a “i 4 os ‘ es fiak mn “f ) : 
"4 IngArehiv 16 (1948), $.364, 0S 
a fee if ay as on . 
eS “ Ba 
: ¥ 24 Sd h aS “os r 


einem idealplastischen Werkstoff 


XVII. Band 1949 Swida: Die elastisch-plastische Biegung des krummen Stabes. 347 


Setzt man in (22), (23) oder in (24) ca=0 und £=1, so ergibt sich der bekannte Ausdruck 
fir die Spannungen bei der. elastischen Biegung: 


TLL aa 


Sere 

Wir ermitteln jetzt das Grenzmoment, d.h. dasjenige Moment, bei dessen Wirkung die 
maximalen Spannungen im Material eine gewisse Grenze erreichen. Darunter verstehen wir 
entweder eine Bruchspannung oder eine solche Spannung, bei deren Wirkung die relativen 
Dehnungen in den dufersten Fasern einen gewissen gegebenen Grenzwert erreichen. Wir be- 
zeichnen diese Grenzspannung mit Om. 


Setzt man in (23) ou=—om und y=~—h,, so findet man das Grenzmoment 
—h Se 
M, = 2 [S,—Su+r (F\—F.)] + —* (om 4) (= nh) ita oe (25) 
1 
Die letzte Gleichung ist mit (20), (21) und (16’) nach y,, r, y, und M aufzulésen. 


Fur Werkstoffe, die einen kleinen Elastizitatsmodul E, haben, kann in einer Reihe von 
Fallen die Bestimmung des Grenzmomentes vereinfacht werden. 

Tatsachlich mu sich bei der Zunahme des Biegungsmomentes die Hohe der elastischen 
Zone verringern, und diese Verringerung muB (bei gleichen Werten der Momente) fir ein Material 
mit kleinem Elastizitatsmodul E, gréBer sein als fiir ein Material mit groBem E,. 

Bevor nun die Spannung den Wert oy» erreicht, kann sich auf diese Weise der EinfluB der 
elastischen Zone praktisch als derart unbedeutend erweisen (insbesondere, wenn om erheblich 
6; und oa tbersteigt), da® man ihn auBer acht lassen kann. Dann kénnen die Gteichgewichts- 
gleichungen (17) und (18) naherungsweise folgendermaBen dargestellt werden: 


E, 6d Oa (Fyu—Fo) 


dp pydF” a 
r+y ; 
(F) 
“ E,ddo ydF 
M =o, (8,—Su) + (5: Su—r [ on (27) 
| (F) 
wobei F' die ganze Querschnittsflache ist. Setzt man in (4) ou=Om und y=—h,, so erhalt man 
E, ody. rh, 28 
dy h, (Om —Ga) « (28) 
Auf Grund von (28) ergibt sich aus (26) und (27) 
; yd EF se Oa hy (Fu—F) = (26’) 
Teshy (Om— a) (r—hy) 
o | r—h 
M, = 00[S,—Su+r (Fy—Fa)] + 5" (om—a4) (S.+S))- (27") 
> 1 \ : 


Aus (26’) kann der Wert von h, bzw. r ermittelt werden. M, ergibt sich ohne weiteres dann 
aus (27’). Setzt man in (18) und (17) Fy, —0 und o,=o; ein, so erhalt man M, fir einen Stab aus 


Me= 0s(So —Su) ’ 


” wobei S, und S, mit ihren Vorzéichen einzusetzen sind. Dabei ist Fy =Fu. Fur einen Stab mit _ 
- rechteckigem Querschnitt ist bh? 


es Ms 7 Os - 


Nun isecen wir fur den Fall, da8 E>E, ist, die Méglichkeit eines weiteren Anwachsens der 


~ Momente und der Formanderungen bis zum vollstandigen Verschwinden der elastischen Zone zu. 


: 
; 


Aus (5), (16) und (20) ersehen wir,’ daB 
| ‘jim 9 OF es Eert ae alae FENG (29) 
(Fi) 


(in oP Goo \ 2 TE died RUN Gre Urs Rag 
; (y2>9) : o) \ 
weil F,,—F, +0 ist. Dann ist auf Grund von (18): ; 
ej: ‘ : d 
| limM=0, wenn “7? 0... 3 (30) 


fy 


h 
3 
rd 
y 1 
‘4 


Die Gisichate (29) reiie daGB, wenn das vollstandige Verschwinden der elastischen Zone mog- 
lich ware, die Grenze der plastischen Zonen des Querschnittes mit der neutralen Achse des ge- 
; Bone we . 99% 


i . f % : ‘ : ia || 
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krimmten Stabes bei der dastil chen Biegung musmonient len miBte. Aus (30) sieht man, dab 
fur das vollstandige Verschwinden der elastischen Zone ein unendlich groBes Moment. erforderlich 
ware. Man muB auch beriicksichtigen, daB bei groBen Detrmesen die von uns zugrunde 
gelegten Annahmen (Ebenbleiben der Stabquerschnitte, Nichtvorhandensein seitlicher Pres- 
sungen zwischen Langsfasern, Unveranderlichkeit der Querschnittsform bei Biegung usw.) in |f 
gewissem Sinne nicht mehr der Wirklichkeit entsprechen. ; 

Unter Verwendung des von uns friher angegebenen Verfahrens! und der Formeln (8) und_ 
(19) erhalten wir folgenden allgemeinen Ausdruck fir die waagerechte Verschiebung fh eines 
beliebigen Punktes A des krummen Stabes infolge der Formanderung des ganzen Teiles OA 


(Abb. 7): 


‘ie M(yo—y)ds [M+o.(Sut+Fi7)] (vo—y) as 
ms a Based | Bo(S+eS.) 
Sk 
? {M+ 0a [(Fu—Fo)r+Su—Sol} (vo—y) a5 
aif Eo [8:4 #(Sut-So)] 


(30) | 


(s),’) 


wobei s; Cie Lange eines elastisch deformier- 
ten Abschnittes des Stabes, s, die Lange 
eines elastischplastisch deformierten Stab- 
abschnittes mit einseitiger plastischer Zone | 
(an der konkaven Flache) und s, die Ab- 
schnittslange mit beiderseitiger plastischer |} 
Zone sind. In Abb.7 sind x und y die |} 
Koordinaten des Gleitpunktes D. 

Ersetzt man in (30’) y, —y durch x)—x, |] 
so erhalt\man den Ausdruck fir die ver- | 
tikale Verschiebung des Punktes A. End- 
lich, wenn man in 30’) ¥o—y durch | er- 
setzt, so ergibt sich der Drehwinkel des 
Querschnittes A. 

Nun wenden wir die obigen allgemeinen | 
Zusammenhange auf den Fall eines ge- | 


Abb. 7. ; krimmten Stabes mit rechteckigem Quer- jf 


te schnitt an. 


4, Rechteckiger Querschnitt. a) Erster Spannungszustand. Auf Grund der Abb. 8 er- 
halten wir: a 


Bi, = b(r—y1~a)i, C95 
\ os b 6 § 
| i OR) Sle ca)? yl 
pydF (upbda c 
: ou a b(ys 4 pee Say) : 
we DPE age eee eae 
FRE POE ae rr) 


Unter Benut 6 ah ig lain eowonienls : onde al 
cea uae Mela von (6) kann man die Gleichgewichtsgleichungen (9) igs (10) foleondce | 


a. [h—(c—y) In =| —o6]h—yy yo —(¢ yg) In $MM —*] = an (te). (31) 
aes i 2. 4 

ces 3b (c—y,— ee 

M (e— 9, ya)? a] = AEE) [Fot-v)—E, [0 —y_)*—yi]}. 82) | 


Die Lésung der erhaltenen zwei Gleichungen mit zwei Unbekannten kann durch die Lésung ; 
einer Gleichung mit einer Unbekannten ersetzt werden, wenn man die Substitution 


Net jee Rb Bose eg ne jer e 
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‘anwendet.’ Dabei sind v und ¢ die neuen Unbekannten. Dann ergibt sich aus (32) nach einigen ' ( 
Umformungen | re Le Hy : ee Ae 
SS ees _ v(A+Bv+oqbv?)+ Deh : 
ped sete Rao 60) , (34) 
_wobei ee Ea : f 
mel van D =2 bh(os—on) .. B=D—2 obec, pier 
Pe - K=2M-+bh(2 0.c—osh), —  'A=2 M+o0b(c? a2) —D(h+e) on 
Bast. oe / ) ; ‘ : ues 
Auf Grund von (33) und (34) kann man die Gleichung (31) folgendermafen darstellen: “lee eae 
; ¥ c—v c h—v : oh 
—2(n—1)h+(2c—v) (In - nin) +» (2—2=*) _ v(AtBo+oqbe®?) + Deh (35) 
a coe ae In — — la K—v(D+oqbv) : ; 


_wobei n =0s/Oa ea Nach Ermittlung von v aus (35) lassen sich t, y,; und Yq mit Hilfe von (34) ‘S « ie 
und (33) errechnen. Danach kénnen die Spannungen aus (11) und (12) gefunden werden. : Se 
& : i ; ! } ‘a , 


—LWHY 


oe 


Su 0kg 


Sih: 7a Pris Saha Wein eee = sei eon ee ~ . 
ig des Wertes von v ‘kann man In =——-in (35) durch folgende, schnell konver- _ tat 


/ ¥ 


pr si Za 
Iplastischen Wer 


kstoff ist Oa =os, und die Gleichungen 


¥ > é 4 
rs 
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uy 
yd F (u—r)bdu 4 TY 
fae J = b(r y,—a—r In s iE 
(F,) 


[uy jos b(e—r ¥o—t ln aa 


Die Werte von Fy, Gad S. sind odie wie im ersten Spannungszustand. 


Mit Hilfe von (6) und (16’ ) lassen sich die Gleichgewichtsgleichungen (20) und (21) folgender- 


mafen darstellen: 


9 2 | 
Jot [ee (1 29192 1, c\ eye ys YaYar yn, “s] (vo+91) —e-a, (36) | 


i Ye—N1 L\ Ga tp am At Jarek Nn YeoT M1 
P | 
M— 3 (ater) + geet {2% [oa(yetyx) —2(0—e) (0.0%) (0.—02) (c2—a)|=0. (37) | 
2 1 {ne 
a, Wir verwenden die Substitution 


J/ o-1 | v+1 t ; : 
: nat | ’ ya ’ (38) — 
wobei t und v die neuen Unbekannten sind. (Die Substitution von (33) fahrt hier nicht zu er-_ 
winschten Ergebnissen.) Es ergibt sich dann aus (36) nach einigen Umformungen | 
v |e 1 (In wi} ;+4ln-~) ; a 
wobei A=o;/0a—1 ist. Die Cletiong (37) erhilt a4 Grund von (38) und (39) folsom Found: | 
F : | heavy e2 26,b{h(1-+A0)+20P | ee 
a—Bo ae yetel 4 et el 0, (40). 
a fb Aln (v?—1) (in “ “) eat 
g=2M~oab b(a*+c?), A= b(c ® a?) (6s—Oa) Yet bla eMei aa.) 
iste. Diese Gleichung liefert den Wert von v, 
Die Werte von t, y,, yj; r und die Spannungswerte kénnen der, Reihe wack aus (39, (38), 


(16’), (22), (23) und (24) ermittelt werden. — \ , 
Wenn fir das Stabmaterial Oa =Os ist, so vereinfachen sich die Gleichungen (39) und (40) | | 


| 


wobei 


-wesentlich: é 
= i c-+a ; 
: a ¥ , ¥ t —— = * > (39) * | 
oer A | nen ay ed ei yet T1n = Sea 
a sa me v+ 1 2.0,b(c--a)? . hic s Me ae | 
2 ae 4 =1)(In oer + 2 M—o,b(2+-a2) * _ (40 ) | 
Fir den Fall des. Grenzmomentes (bei kleinem F,) ergibt sich nach ead | 
¢ ‘ "i 2 Rak ae, 2a , ‘ a 
: pias ot Ve 8 Cae ts ee Pere re ss | 
eRe Og hy \ ; A Gq are (41) | 
‘> cworln ‘i igs fl ae ae co teh tna a ie wee | 
aS + uf , eal i x os * A : pe L ; a 
x6 1; i i ¥yY => dete In ma sat ry mG) = a(om— Ga) (a1 In = ae eee } 4 
Bee ren NN Se eg | 


eee, Nach Brmitaang von h, ‘Kann der Wert von r=a+h, und M, er) berechnet werden. 


Poe Beuipiele.” a) Erstes Beispiel. Der rechteckige Querschnitt eines gekriimmten Stahes x 
hahe die Abmessungen b=5 cm, h=10 cm. Es wird angenommen M= -2tm,a =10 cm, ¢ =20. cm, 4 
0. =2400 kg/em2, oa —2160 kg/em?, E=2-108 kg/em?, €=0,1. Es sind die Lage der neutralen, | 
Achse der Grenze zwischen der ae und as plastischen Zone und ie este tey * einer _ 
ye ee Faser zu ermitteln: ) “S< 5 

Das Moment, das dem F NeBbceint, entspricht, betragt is (13) 1,55 tm. Bera Be 


Nomen von 2tm tritt der erste, pags bin eelacas aul.” Die Gee hinass (35), ee | 
iefern . 


vs 9,20 em, t= 1,86 ems 3123.67 em, r4=5.58 om, r= 1447 em. 


“xvi ms : 


Die Hohe der. plastischen Zone betragt 0,80 cm. Die ree ergeben sich nach (11) und ie as 
(12) zu a 
_ _ 1044 y ce 04,4 y ae 


Der eS canenesvelent ist in Abb. 10 dargestellt. - 


b) Zweites Beispiel. Die Abmessungen des Querschnittes und das Material des Stabes 
sind dieselben, jedoch wird M =2,8 angenommen. Bei diesem Moment tritt der zweite Spannungs- 
zustand auf. 

Die Gleichungen (40), (39), (38) ‘und (16’) liefern ale. nach v =5,40, t = 2,78 cm, 

41 =2,31 em, y,=3,35 cm, r=14,76 cm. : ees 
Die obere plastische Zone hat eine Héhe von 1,89 cm und die untere von 2,45 cm. Die Span- ele: 
nungen werden auf Grund (22), (23). und (24) 


; 12988 y 1298,8 y 
: Da TRI A Ou= 2160 + 1416? 

Der Spannungsverlauf ist in Abb. 11 dargestellt. = 

- Es ist dabei zu bemerken, da8B das von der elastischen Zone aufgenommene Moment in diesem 
‘Fall nur 25 v. H. des gesamten Momentes betragt. \ ar a 
_-Man ermittelt nun das Grenzmoment M,, wenn om—6160 kg/cm? ist. Bei dieser Spannung 
Pires. die relative Verkiirzung der unteren AuBenfasern 2 v. H. Die Naherungsformeln (41) 
und (27’) liefern h, =4,67 cm, r = 14, a em und M;=5,17 tm. Der Spannungsverlauf i istin Abb. 12 - 
- dargestellt. ae ; i : 


2 


eat 1298,8y . 


4 <4 . ‘ x i 1 f 5 i 


z 
iy tft Baal 
& 


4 


qe AUS ee oo 
HHT 
a oe 

\\ 


“ 


Se 


ikl yy Abb. 12, Spininatiaibets fir das 
ay a5 Grenasterent Waa 


oe 


j 


batrubti ae chet de F LENE | pene werd en 
wir hier nur, daB die Zulassung der plastischen Defor- 
- ergeben kann, wie man aus. ee Hi weiner unten, dar-) Ms 
Stab mit b=5 cm, 
=15. cm ein Moment 
(2475, kg/e cm? betragen 
ng an unteren } Bearer a 


he bhew: amin e€ 
netr scl n Achse v 
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annimmt und die Beziehung h/b=2 und den Krimmungshalbmesser 9 = 15 beibehalt. Der oben: | 
angenommene Wert 6,y,= 2000 kg/cm? ist fiir eine ,.elastische’*’ Berechnung -ziemlich hoch, 
weil o,=2400 kg/cm? ist. Die Materialersparnis betragt infolge der Zulassung des elastisch- 
plastischen Zustandes 39,2 v. H. Die relative Verkiirzung der unteren Fasern ist 0,1 v. H. bei 
einer Spannung von 2000 kg/cm?. 

Fir das Moment M =2,8 tm betragt die gréBte Spannung 2778 kg/cm? bei b=5 cm, h=10 cm 
und ¢=15 cm (Abb. 11). Dementsprechend ist die relative Verkirzung 0,31 v. H. Der Sicher- 
heitsgrad sinkt bis 1,85, was in manchen Fallen noch ge- 
nigen diurfte. 

Auf Grund der ,,elastischen‘‘ Berechnung erhalten wir, 
wenn M=2,8 tm, 6,1 = 2000 kg/em?, h/b=2 und g=15 cm 
ist, h=13,4 cm und b=6,7 cm. Die Materialersparnis be- | 
tragt 79,6 v. H. S . 

Als wesentlicher Faktor zugunsten der Zulassung eines | 
elastisch-plastischen Zustandes in. den vorher betrachteten | 
Fallen erscheint die geringe GroéBe der dabei entstehenden 
Formanderung des Stabes, worauf die Zahlenwerte der rela- 
tiven Verkirzungen der AuBenfasern hinweisen. 

Es ist jedoch zu bemerken, da fiir die Mehrzahl der in | 
der Praxis verwendeten Stahle ohne FlieBteil in den (o, €) | 
Diagrammen der Zahlenwert € kleiner als der in unseren | 
Beispielen angenommene Wert =0,1 ist. Mit der Vermin- | 
derung von & wachst die Dehnung der Faser, die der gleichen | 

Abb. 13 Spannung entspricht. 
Zum SchluB zeigen wir den allgemeinen Weg zur Lésung 
der Aufgabe fir den Fall, daB die Spannungs-Dehnungs-Linie bei Spannungen tber o; einen | 
krummlinigen Verlauf hat (Abb. 13). Die Gleichungen der Geraden pm und der Kurven mn 
und gp sind ; 
| | 


o=E-e o=fi(e); o=fy,(e)- 
Die Dehnungen kénnen durch die Formeln (2) ausgedrickt werden. Es wird angenommen, dab 
das Verhalten des Werkstoffes sowohl beim Zug als auch beim Druck gleich ist. . 
Die Gleichgewichtsgleichungen fiir den zweiten Spannungszustand erhalten folgende all- 
gemeine Form: ° : 


fucar+ [pw@ar+ [fear =o, 


(F.) (F,) _ & 
fEeydF+ [ fileydF+ | fle)ydF=M- : 
(F,) (F.) 148) 


Im ersten Zustand entfallt das dritte Integral in jeder dieser Gleichungen. Auer den Gleich- 
gewichtsgleichungen miissen noch die Bedingungen eines stetigen Uberganges fiir die Span- 
nungen an den Grenzen der elastischen und der plastischen Zonen erfillt werden. 


(Eingegangen am 4. Oktober 1948.) 
Anschrift des Verfassers: Dr.-Ing. Waldemar Swida, (17a) Karlsruhe, Stuttgarter Str. 21. 
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